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ABSTRACT 


Ice crystals were grown in a supersaturated environment pro- 
duced by a dynamic thermal diffusion chamber, which employed two 
horizontal plates separated by a distance of 2.5 cm. Air was circu- 
lated between and along the 1.2 m length of the plates past ice 
crystals which nucleated and grew from a fiber suspended vertically 
between the two plates. Using a zoom stereo microscope with a magni- 
fication which ranged from 3X to 80X and utilizing both 35 mm still 
photographs and 16 mm time lapse cine films taken through the micro- 
scope, the variation of the shape and linear growth rate of ice crystals 
was examined as a function of the ambient temperature, the ambient 
supersaturation and the forced ventilation velocity. The ambient growth 
conditions were varied over the range of temperature 0°C to -40°C, 

over the range of supersaturation 4% to 50% with respect to ice, i.e., 

-3 -3 

over vapor density excesses ranging from 0.07 g m to 0.7 g m , and 
over the range of forced ventilation velocities 0 cm s”^ to 20 cm s - ^ . 

It is shown that the introduction of a ventilation velocity is 
roughly equivalent to increasing the ambient supersaturation. For a 
fixed ambient temperature and ambient supersaturation, the linear 'a 1 - 
axis growth rate is directly proportional to the square root of the 
ventilation velocity, as theory predicts, provided the crystal shape 
does not change significantly. The transitions plate dendrite and 
column needle occur at a lower ambient supersaturation as the venti- 
lation velocity increases. A definite time constant, which is a 
function of the ambient temperature, the ambient supersaturation and 
the magnitude of the change in the ventilation velocity, exists for 
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the transition of both crystal shape and linear growth rate following 
a change in the ventilation velocity. For increasing ventilation 
velocities at a fixed ambient supersaturation the maximum in the linear 
growth rate near -15°C apparently occurs at successively colder tem- 
peratures. Over the temperature range -4°C to -6°C growth occurs 
along a direction up to 25° from the 'c'-axis as the local supersatura- 
tion is increased. Thus, under the proper temperature conditions a 
change in the local supersaturation can induce a change, not only in 
the absolute growth rates, but also in the relative growth rates along 
the 'a* and 'c'-axes. 

In the presence of 10 um mean diameter droplets with concentrations 

o c o 

of 10 to 10 cm droplet accretion accounted for over 90% of the 
growth of both ice crystal columns and dendrites at ventilation veloci- 
ties of 15 cm s -1 . However, even at higher velocities the most exten- 
sively rimed crystals still retained the original orientation of their 
crystalline axes. At velocities less than 1.0 cm s - ^ droplets of 
10 ym diameter or smaller evaporated as they approached a growing ice 
crystal and crystal growth was entirely by vapor diffusion. 

The first vapor grown discoid ice crystals were observed. They 
grew in the temperature regime -5°C to -7°C at low local supersatura- 
tions, i.e., in the regime formerly believed to only support nearly 
equiaxed columns. 

Results from these experiments are interpreted in terms of diffu- 
sion through a local boundary layer, whose thickness is a function of 
ventilation velocity, the diffusivity of water vapor and heat through 
air, and the crystal shape; and interaction with different nucleation 
and growth kinetics on different surfaces. Two-dimensional nucleation 
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and layer growth from corners or edges ts believed to occur at ambient 
supersaturations significantly lower than present theories would pre- 
dict. At low supersaturation and temperature the crystal habit and 
growth may be controlled by the presence of defects sometimes giving 
rise to crystals of habit opposite to that normally observed. 
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I. 

INTRODUCTION 

The growth of ice crystals is of practical interest in meteor- 
ology because of the major role which they play in precipitation 
processes. In fact, nearly all precipitation which falls from clouds 
over large land masses, whether in winter or summer, is at one time 
during the precipitation process composed of the ice phase. On the 
other hand, the growth of ice crystals is also important in providing 
information about crystal growth processes. The fact that ice is 
inexpensive, has a wide range of well-defined forms or habits and can 
be grown without the use of expensive furnaces normally associated with 
crystal growth, makes it an ideal choice for intensive basic crystal 
growth studies. 

Ice can be formed by either of two processes: solidification 

from the liquid phase or condensation from the vapor phase. If 
solidification takes place in a 'pure' substance it is called growth 
from the 'melt 1 . Otherwise it is called growth from 'solution'. 
Crystals grown by condensation from the vapor phase are called vapor 
grown crystals. In either case, crystals will only grow if the solu- 
tion or vapor in contact with them is supersaturated or, in the case 
of a melt, undercooled. 

The actual form, or habit, of ice crystals may vary widely, 
these differences being dependent on the conditions under which the 
growth has taken place. In fact, the growth of ice crystals and their 
resultant habits is determined by a myriad of variables. Nonetheless, 
these variables may be divided into two rather broad and sometimes 
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overlapping categories: those which primarily influence the flux of 

water molecules to, and the transport of heat from, the crystal sur- 
face; and those which primarily govern the distribution of the molecules 
on the crystal surface and their subsequent incorporation into the 
crystal lattice. The former category might be termed "environmental 
or ambient kinetic effects" and the latter "surface kinetic effects". 

The resultant crystal habit will be determined by an integration over 
both time and space of the relative magnitudes of these two terms. Of 
course, it should be realized that the relative magnitudes of these 
terms may vary significantly from one crystallographic face to another 
on the same crystal . 

In this work emphasis is placed primarily on vapor grown ice 
crystals. Of particular interest is their habit and growth rates over 
the range of ambient conditions of temperature, pressure, supersatura- 
tion and ventilation velocity existing in the earth's atmosphere. 

This is important because the habit of the ice crystals not only affects 
their linear growth rates and riming, i.e., accretional growth char- 
acteristics, but may also be important in determining secondary crystal 
production. These characteristics, in turn, determine the rate at 
which precipitation can form in a cloud. 

In the atmosphere, the ambient variables which determine the 
crystal habit and growth rate may fluctuate with both space and time 
in some complicated manner. Therefore, it is profitable to perform 
experiments in the laboratory under as well controlled conditions as 
possible. These laboratory experiments provide insights into the 
fundamental growth characteristics of the crystals; they help verify, 
modify or disprove existing theories; and they complement both surface 
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and airborne field investigations. The extent to which these labora- 
tory studies can be directly applied to the atmosphere depends upon 
the effectiveness with which the laboratory design simulates the 
atmospheric process. In practice, experimentally simulating the 
variations which each of the ambient parameters undergoes in a cloud 
and at the same time providing a method for close observation of the 
growing ice crystal is a virtually impossible task. Thus, in any ex- 
periment designed to study the growth of ice crystals one must first 
determine which parameters he wishes to study; then design and build a 
system in which these parameters are controlled and the remaining 
parameters are held constant at some appropriate value. 

In the past, experiments have been performed in the laboratory 
with temperature, pressure and supersaturation as controlled variables 
but with ventilation dependent upon free convection. Other experiments 
have been performed in which ice crystals fell through a supercooled 
water cloud with temperature as the controlled variable and ventilation 
dependent upon the terminal velocity of the particles. In this case 
growth occurs at, or near, ambient conditions of water saturation and 
an ice supersaturation which depends only on temperature. Still other 
experiments have been performed in which an ice crystal was grown on 
a substrate where temperature, air pressure and supersaturation were 
controlled. However, in this case the substrate dominates the heat 
flow and ventilation is small and uncontrolled. 

The advantage of the chamber used in this study is that it has 
the capability of independently simulating the ambient temperature, 
pressure, supersaturation and ventilation velocity, all under controll- 
ed conditions. To date, no other laboratory instrument has had this 
capability. 
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The purpose of this particular study was to attempt, through the 
use of a controlled environment, to answer some of the questions asso- 
ciated with the growth of ice crystals. I was most interested in deter 
mining how ice crystal growth varies as a function of ambient tem- 
perature, supersaturation and, particularly, ventilation velocity, i.e. 
the fall velocity of the crystal; and how the transition from plates 
to dendrites and from columns to needles is dependent upon the ambient 
supersaturation and the ventilation velocity. I am also interested 
in the growth of ventilated ice crystals in an environment with drop- 
lets, and in the habit of ice crystals at low supersaturation. 
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II. 

REVIEW OF PREVIOUS WORK 

A history of early observations and sketches of ice crystals, as 
detailed by Heilman (1893), is reported in the books of both Nakaya 
(1954) and Hobbs (1974). According to this history, observations and 
sketches of ice crystals were made as early as 1550. Intricate crystal 
structure, however, was first observed following the discovery of the 

i 

microscope in the latter half of the seventeenth century. Even so, 
not until the discovery of photography could the crystal structure be 
accurately and quickly recorded. Following the arrival and widespread 
use of photography, Bentley (1931), using photomicroscopy, produced 
over 6000 spectacular photographs showing the broad variation of forms 
that exist in natural ice crystals. 

Despite these early studies, only recently, with knowledge of 
molecular processes and crystal aspects, has an attempt been made to 
reach a comprehensive understanding of ice crystal growth mechanisms 
and the role which they play in the determination of precipitation. 
Ukichiro Nakaya and his fellow workers of Japan were the first to 
undertake such a project in the laboratory, beginning in 1932. Nakaya 
devised an experiment in which he grew snow crystals in his laboratory 
under controlled temperature and supersaturation conditions. From 
these studies he was able to produce for the first time artificial 
snow crystals which closely resembled those found in the atmosphere. 

His books and various papers contain hundreds of photomicrographs, not 
only of these snow crystals, but also of natural crystals. Nakaya 
made a distinction between a snow crystal and an ice crystal. He 
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defined the ice crystal as the crystal of ice formed in the air by 
condensation of water vapor. This ice crystal then represents what he 
called the 'germ' for the snow crystal. He also called the crystals 
formed from supercooled water droplets by seeding, ice crystals. He 
believed that the form and structure of the snow crystal that is grown 
from this germ of ice or ice crystal are determined by the meteorologic- 
al conditions met after this germ is formed. Although he realized that 
it is not possible with this designation to draw a distinct line be- 
tween snow and ice crystals, he believed that since it was possible 
to grow a snow crystal under different meteorological conditions from 
those under which it was nucleated that this designation would be use- 
ful in eliminating possible confusion between the two processes. This 
distinction is still made by many from the 'Japanese school'; however, 
the 'European school', for the most part, have made no differentiation 
between the two terms and often use them interchangeably. In this 
study no distinction will be made between the two terms. 

Paralleling the early studies of snow crystals by Nakaya and his 
fellow workers were atmospheric observations and both laboratory and 
theoretical studies in various other parts of the world. These obser- 
vations and studies have continued to the present time. The list of 
those participating and their contributions over the years is pro- 
digious. However, detailed histories and summaries of these contri- 
butions may be found in Mason (1971), Lamb (1970) and Hobbs (1974). 

A . Ice Crystal Terminology and Ice C rys tal Habit 

The external shape of an ice crystal is determined by the relative 
rates of growth of the different crystallographic faces. That is, 
those crystallographic faces or planes having the fastest rate of 
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growth also disappear most readily. This fact results in the observed 
most prevalent ice crystal growth form under normal temperature and 
pressure conditions being a prism bounded on the two ends by the 
slowly growing {0001} "basal" faces and on the six sides by the slowly 
growing {1010} "prism" faces. The crystallographic faces are desig- 
nated here by their Mil ler-Bravais indices. Figure 2.1(a) and (b) 
illustrate this most prevalent growth form. Crystallographic directions 



Fig. 2.1 Schematic representation of possible ice crystal growth 
forms. Wolff (1955) 

are indicated by enclosing the Miller-Bravais indices in square 
brackets. For the hexagonal symmetry of ice the crystallographic 
direction [0001] indicates a direction perpendicular to the basal 
face (0001) and is called the 'c'-axis. The basal face contains three 
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'a'-axes oriented 120° from each other. That is, the crystallographic 
direction [1120], for example, which is perpendicular to the crystal- 
lographic face 0120) in Fig. 2.1 Cd), indicates the direction of one 
of these 'a'-axes. Furthermore, due to symmetry considerations the 
'a'-axes are indistinguishable so growth along any of the 'a'-axes can 
be said to be along the <1120> crystallographic direction. Figure 2.2 
shows the unit cell for the hexagonal ice structure. The basic crystal 



Fig. 2.2 The unit cell for the hexagonal ice structure. Gamara (1972) 

habit is determined by the ratio of the lengths along the ' c ' and 'a'- 
axes. If the ratio c/a is large the crystal is prismatic. Fig. 2.1(a). 
Whereas, if the ratio c/a is small the crystal is plate-like, Fig. 2.1 
(b). Ice crystals are commonlyclassified into one of five basic 
crystal habits: hexagonal solid prisms, hexagonal hollow prisms, 

needles, hexagonal plates and dendrites. Figure 2.3 illustrates these 
five basic habits. 
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HOLLOW PRISM 



PLATE NEEOLE DENORITE 


Fig, 2.3 Observed basic ice crystal forms. Gamara (1972) 

Much effort has been expended in delineating the ice crystal habit 
as a function of the ambient atmospheric conditions. It is now general- 
ly conceded that temperature is of primary importance and ambient 
supersaturation is of secondary importance in determining the ice 
crystal habit. These findings are the culmination of several labora- 
tory experiments and a multitude of atmospheric observations. 

0 <e of the first attempts to display the habit of ice crystals as 
a function of both the ambient temperature and supersaturation, or 
vapor density excess, was made by Nakaya. These diagrams are often 
referred to as 'Nakaya diagrams'. Figure 2.4 is such a diagram and 
shows the effect of both temperature and supersaturation on ice crystal 
habit. This diagram is based largely on the laboratory work of Hallett 
and Mason 0958a). If instead of graphing percent supersaturation 
relative to ice, one graphs vapor density excess over ice, one obtains 





10 



Fig. 2.4 Ice crystal habit as a function of temperature and super- 
saturation relative to ice. From Mason (1971) 

a diagram like Fig. 2.5 which shows the results obtained by Kobayashi 
(1958) in his laboratory experiments. Likewise, Magono and Lee (1966) 
have presented this same type 'Nakaya diagram' more pictorially, as 
illustrated in Fig. 2.6, for natural crystals occurring in the 
atmosphere. One of the most striking characteristics common to all 
three of these 'Nakaya diagrams' is the fact that near water satura- 
tion the basic crystal habit makes three transitions between 0°C and 
-30° C. That is, the basic crystal habit is plates from 0°C to -3°C, 
prisms from -3°C to -8°C, plates from -8°C to -25°C and prisms again 
for temperatures colder than -25°C. There is some disagreement about 
the absolute value of the temperatures, for example within 1°C or 2°C, 
at which these transition boundaries occur, but their existence is 
well established. The transition boundaries at -3°C and -8°C are 
sharp, whereas the transition boundary near -25°C is more gradual. 

The effect of ambient supersaturation is to exaggerate the 


GRADE OF VAPOR SUPPLY 

INCREASE =* 

p CLOUD DROPLET REGION 


n 


Direction of preferred growth 



Temperature in C C 

Fig. 2.5 Ice crystal habit as a function of temperature and vapor 
density excess over that of ice. Kobayashi (1958) 



0 - 5 -10 -15 -20 -25 -30 -35 -40 CC ) 

TEMPERATURE 

Fig. 2.6 Pictorial diagram showing the dependence of ice crystal type 
upon both temperature and supersaturation relative to ice. 
Magono and Lee (1966) 
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extremes in habit. In fact, over the temperature region 0°C to -30°C 
and at large values of ambient supersaturation, the ratio of crystal 
axes, c/a, may vary over such a wide range of values as .01 for large 
thin plates to approximately 50 for long narrow needles. 

A comprehensive understanding of ice crystal growth must include 
an explanation of why, near water saturation, and within specific 
temperature ranges, the growth should take place preferentially along 
the 'a'-axis and, at other temperatures, along the 'c'-axis. In any 
attempt to fully explain these different relative growth rates it is 
necessary to consider the growth mechanisms by which the crystal faces 
are propagated. Hence, an account must be made not only of the 
crystal structure of ice and its equilibrium forms, but also of the 
surface kinetic possibilities of adsorption and diffusion of molecules 
across the surface of the individual growing crystal faces, of the 
accommodation of these molecules into the crystal lattice, and of the 
temperature and vapor diffusion fields which surround the crystal. 

In the past, several experiments have been devised and performed 
in an attempt to gain an insight into how the surface kinetics of the 
crystal and other relevant parameters interact under different ambient 
conditions to produce the observed crystal habits. The recurring 
difficulty in these experiments has been in accurately determining the 
actual conditions in which the ice crystal growth took place. 

In general, the fewer the parameters involved in the experiment, 
the simpler it is to separate the effects due to one individual para- 
meter. However, the fewer the parameters involved the less the ex- 
periment simulates actual crystal growth conditions in a cloud. I 
have chosen to separate my discussion of previous experiments into two 
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classes: those experiments utilizing a static environment, that is, 

the carrier fluid flow is a minimum and those experiments utilizing a 
ventilated environment. 

B. Ice Crystal Growth in a Static Environment 

The first experiments dealing with ice crystal growth in a semi- 
static environment were conducted by Nakaya (1954) in a convection 
chamber. He grew the crystals on a fine rabbit's hair suspended in a 
cylindrical cold-chamber. A beaker of water, the temperature of which 
could be controlled and maintained as warm as desired, was located at 
the base of the chamber and served as the moisture source. A fog, 
produced over the warm water, rose due to natural convection, becoming 
supercooled before reaching the ice crystals suspended in its path. 
Therefore, the ice crystals grew in an ambient environment which was 
near water saturation but which also contained small water droplets. 
The ventilation velocity or carrier fluid velocity, although probably 
less than 5 cm s - ^ , was not insignificant. The presence of the water 
droplets and the ventilation velocity made calculation of the "true" 
supersaturation of the environment impossible. So the results were 
not only difficult to interpret, but, because the precise effect of 
neither the ventilation velocity nor the presence of the water drop- 
lets was known, it was also difficult to compare these results with 
those of other experiments where no droplets were present and the con- 
vection velocities were much smaller. Even though it was difficult, 
in these experiments, to separate the independent effects of the 
ambient supersaturation, the ventilation velocity and the presence 
of the water droplets on the crystal habit, Nakaya was able to extract 
enough information to compose the first 'Nakaya diagram'. 
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Later laboratory experiments by Kobayashi (1957), Hallett and 
Mason (1958a), Hallett (1965), and Rottner and Vali (1974) employing 
static thermal diffusion chambers have helped to better establish the 
relationship of ice crystal habit to ambient temperature and super- 
saturation. 

The basic operating principle of the static thermal diffusion 
chamber is the same as that for the dynamic thermal diffusion chamber 
used in this study. The static thermal diffusion chamber consists of 
a sealed chamber with two horizontal plates separated by a vertical 
distance of 1 to 20 cm. The temperature of the plates is controlled 
independently with the top plate being maintained at a warmer tempera- 
ture than the bottom plate to suppress convection due to thermal in- 
stability. If the warmer top plate is coated with ice or water, it 
acts as a water vapor source and the bottom plate as a vapor sink. 

The density of the water vapor at the two ice plate surfaces is equal 
to the equilibrium vapor density of ice at that respective temperature. 
Both the density of the water vapor and the ambient air temperature 
vary almost linearly with height between the two surfaces. However, 
the equilibrium vapor density, at any height, is an exponential 
function of temperature and thus the air between the two plates is 
supersaturated with respect to ice. (For a more detailed account of 
this operating principle see Chapter IV). The maximum supersaturation 
occurs near the center of the chamber and may be changed by varying 
the temperature difference between the top and bottom plates, larger 
temperature differences corresponding to higher supersaturations. 
Therefore, if liquid water is used as the vapor source at the top 
plate, instead of ice, larger temperature differences can be sustained 



15 


and much larger supersaturations may be obtained in the chamber. At 
supersaturations much greater than water saturation, any aerosol's or 
particles present in the chamber nucleate and form water droplets which 
grow and then settle to the chamber floor due to gravity. When such 
nuclei have all settled out, very high supersaturations, several 
hundred percent with respect to ice, can be achieved such that even 
the ions produced along paths of cosmic rays nucleate water droplets. 
Ice crystals are grown on a fine fiber of glass, nylon or hair suspend- 
ed vertically in the center of the chamber. A diagram of the modified 
static thermal diffusion chamber used by Hallett and Mason (1958a) is 
presented in Fig. 2.7. 


DIFFUSION CLOUD CHAMBER 



Fig. 2.7 The static diffusion cloud chamber used by Hallett and 
Mason for growing ice crystals at various temperatures 
and supersaturations. Mason (1963) 


One of the chief advantages of both the convection chamber and 
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the thermal diffusion chamber over many of the other artificial or 
natural environments in which ice crystal growth has been studied is 
the fact that ice crystals can be grown and observed in them over long 

O 

periods of time. The ambient growth conditions can either be main- 
tained relatively constant or they can be altered during the period of 
crystal growth. 

Nakaya (1954), and later Kobayashi (1957), changed the ambient 
conditions of supersaturation and temperature in their convection 
chambers during the period of crystal growth. Likewise, Hallett and 
Mason (1958a), by varying the vertical distance between the horizontal 
plates of their modified diffusion chamber. Fig. 2.7, were able to 
change the ambient supersaturation without changing the temperature in 
the very center. A primary dependence of the crystal habit on tem- 
perature with a secondary dependence on supersaturation was found by 
all. Further, they all found that when a growing crystal was intro- 
duced to a different ambient environment the new growth immediately 
assumed the habit characteristic of that new environment. This strong- 
ly suggested that the ice crystal habit is a product of the ambient 
environment and of the crystal growth processes but not of the nuclea- 
tion of the seed crystal. The secondary dependence of the habit on 
supersaturation was dramatically illustrated by the fact that at a 
constant temperature, near -14°C, dendritic growth was only observed 
when the air was significantly supersaturated with respect to water. 

At somewhat lower supersaturations, but still greater than water 
saturation, crystals which Hallett and Mason (1958a) called "dorites", 
i.e., spear-like, grew. At progressively lower supersaturations 
sector plates, then plates, and finally, nearly equiaxed crystals grew. 
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It should be noted that words such as dendrites, dorites and 
sectors, as applied to ice crystals, are descriptive terms and not 
quantitative assessments of either the crystal habit or linear growth 
rate. In fact, these terms may have quite dissimilar meanings to 
different individuals. Nakaya (1954) separated the crystals which 
were arbitrarily called dendrites into two classes: "dendritic" and 

"sector" crystals. He attempted to classify the crystals according to 
the branch tip, that is, according to the final ambient conditions 
encountered in its formation. If the branch tip was a straight-line 
structure he called the crystals sector. Figure 2.8(a) and (b) illu- 
strate these designations, whereas Fig. 2.8(c) illustrates the widest 
form of the crystal which Hallett called a dorite. Although these 



Fig. 2.8 Distinction between dendritic, sector and doritic crystals. 
Nakaya (1954) 

(a) 'true' dendritic crystal 

(b) dendritic crystal with sector-like structure 

(c) 'dorite' crystal 


definitions are useful, they are somewhat ambiguous. It is sometimes 
helpful, although awkward, to Include another adjective in the 
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descriptive term and thus describe the crystals as "sector plates", 
"sector dendrites", and so forth. 

Still other experiments have been performed in a static environ- 
ment in which ice crystals were grown on a substrate. The experiments 
of Shaw and Mason (1955); Bryant, Hallett and Mason (1959); Hallett 
(1961); Mason, Bryant and Van den Heuvel (1963); Kobayashi (1965); 

Lamb (1970); and Anderson (1974) are representative of this type of 
study. 

In the experiments of Shaw and Mason (1955) and of Lamb (1970) 
the ice crystals grew on a smooth metal surface, whereas in the experi- 
ments of Bryant, Hallett and Mason (1959); Hallett (1961); Mason, 

Bryant and Van den Heuvel (1963); Kobayashi (1965); and Anderson (1974) 
the ice crystals grew on the face of a crystal of dissimilar material. 
In such special growth, where crystals of one material grow in a 
definite orientation on a crystal face of another material they are 
said to exhibit "oriented overgrowth" or "epitaxy". Thus, epitaxial 
growth of ice crystals is said to occur when the direction of the 'c 1 - 
axis of the growing ice crystal bears a simple relationship to the 
crystallographic orientation of the substrate. Bryant, Hallett and 
Mason (1959) studied the epitaxial growth of ice crystals on such 
varied materials as Agl, CuS, Cd^, Pb^, ^ 2 ^ 5 ’ brucite, calcite, mus- 
covite, orthorhombic Hg^ and orthorhombic iodine. However, the 
crystals most commonly used for ice crystal epitaxial growth have been 
Agl, Pb^ and CuS (covellite). 

Bryant, Hallett and Mason (1959) were the first to observe that 
very thin ice crystal plates growing epitaxially on fresh {0001} 
covellite cleavage surfaces exhibited interference colors in reflected 
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white light. These interference colors can be calibrated and hence 
the thickness of the growing ice crystals can be determined. Other 
epitaxial substrates also exhibit this interference phenomena. 

It has been found that at low and moderate supersaturations, 
these epitaxial ice crystal plates will often grow laterally for long 
periods of time without thickening, that is, without a detectable 
change in interference color. However, if the growing crystal comes 
in contact with either another crystal or an obstacle, then successive 
new layers of growth spread from the point of contact. Likewise, if 
the ambient supersaturation is suddenly increased or a large electric 
field is introduced, the crystal begins to thicken but, in these cases, 
much more uniformly. This thickening process has been interpreted as 
the propagation of steps across the crystal surface. Below a critical 
supersaturation these steps are formed by a process of heterogeneous 
nucleation. At high supersaturations, however, the thickening process 
has been interpreted as successive two-dimensional nucleation followed 
by lateral growth. 

Hallett (1961), and later others, found that the speed of propa- 
gation of these steps, a few hundred angstroms in height, is inversely 
proportional to the step height for a given temperature and ambient 
supersaturation. This result is interpreted as meaning that water 
molecules impinging on the exposed ice crystal surface are not incor- 
porated immediately into the crystal lattice but may be adsorbed and 
migrate short distances, i.e., a few pm, over the crystal surface 
before being either incorporated into the steps or evaporating, i.e., 
being desorbed, from the surface. For a constant ambient vapor density 
excess and a given step height, Hallett (1961) found a remarkable 
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temperature variation, shown in Fig. 2.9, in the speed of propagation 



Fig. 2.9 The temperature variation of the rate of growth of a layer 

250 A thick. Excess vapor density 0.25 g m _ 3. Hallett (1961) 

of these steps. Later experiments by Mason, Bryant and Van den Heuvel 
(1963) and Kobayashi (1965) confirmed the general shape of the curve 
found by Hallett (1961) but suggested it should be shifted to colder 
temperatures by about 4°C. However, Ryan and Macklin (1969) pointed 
out that essentially the same equipment was used in all these experi- 
ments and although Hallett (1961) adjusted his experimental data with 
an empirical temperature correction of 4°C, which was apparently neces- 
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sary, the other experimenters did not. Hence, the discrepancy in the 
results is probably not real. 

There are problems, as with other experiments, associated with 
ice crystal growth on a substrate. The substrate dominates the heat 
flow and the ventilation is uncontrolled. Furthermore, when thin ice 
crystals grow on a substrate the local atomic fitting between the ice 
crystal and the substrate is the primary factor initiating the deposit 
orientation. As the ice crystals grow laterally, the degree of lattice 
misfit must affect the stress and stability of the growing ice crystal 
and possibly induces dislocations which could significantly affect the 
growth. Anderson (1974) and Anderson and Hallett (1977), for example, 
interpreted some of their results, particularly for thin crystals, as 
being induced by the presence of this crystal strain. Hence, growth 
on a substrate does not fairly simulate the conditions encountered by 
an ice crystal growing in the free air. 

C. I ce Crystal Growth in a Ventilated Environment 

Marshall and Gunn (1955) were among the first to postulate that 
the carrier fluid velocity or ventilation velocity could play a signi- 
ficant role in determining the 'effective' vapor density excess over 
the crystal and thus its growth. Later, Hallett (1965) suggested that 
a ventilation velocity would cause the supersaturation dependent 
habit transitions, for example, plates to dendrites, and prisms to 
needles, to occur at a lower ambient supersaturation than found in the 
static case. However, the only experiments, prior to studies with 
the dynamic thermal diffusion chamber, designed to study the effect of 
ventilation on ice crystal growth from the vapor have involved ice 
crystals falling through a supercooled water cloud. Representative 
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of this type experiment are the studies of Aufm Kampe, Weickmann and 
Kelley (1951); Reynolds (1952); Mason (1953); Isono, Komabayasi, 
Yamanaka and Fujita (1956); Fukuta (1969); and Ryan, Wishart and Shaw 
(1976). This type of experiment is generally carried out in a cold 
chamber in which a cloud of supercooled water droplets is "seeded" by 
some mechanism to produce ice crystals which then grow and fall out. 

As the crystals grow, they deplete the vapor supply and thus lower the 
ambient supersaturation. Therefore, droplets must be added to main- 
tain an ambient environment near water saturation. Also, as the 
crystals grow their terminal fall velocity and hence their ventilation 
velocity increases. Furthermore, the presence of the droplets and 
ice crystals makes the 'effective' supersaturation near the crystals 
difficult to determine but even more difficult to control. So, although 
these studies closely simulate the crystal growth conditions encounter- 
ed in a cloud, they do not allow close observations of the growing 
crystals, very long growth times (maximum of about 200 seconds), con- 
stant environmental growth conditions, or a separation of the related 
but distinct effects of droplets and the ventilation velocity. None- 
theless, these experiments have proven very valuable in providing 
further evidence that ice crystal habit is a primary function of tem- 
perature with a secondary dependence on supersaturation. Moreover, 
these experiments have not only shown the manner in which ice crystal 
mass growth rate varies as a function of the ambient temperature, but, 
in addition, have provided valuable data which has been used in pre- 
dicting mass growth rates from empirical models. 

Experiments not directly related to ice crystal growth from the 
vapor, but which show, at least qualitatively, the effect of ventila- 
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- tion velocity on crystal growth have been performed in both melts and 
solutions. 

Several investigations, Kumai and Itagaki (1953); Lindenmeyer 
(1959); Hallett (1960); Knight (1962); Hallett (1964); Camp (1965); 
Macklin and Kyan (1965); Knight (1966); Lindenmeyer and Chalmers 
(1966) and Pruppacher (1967), to name a few, have been made of ice 
crystal growth from the melt and the solution under static conditions, 
both growing freely in the fluid and on a substrate. 

Other experiments, Fernandez (1967); Fernandez and Barduhn (1967); 
Poisot (1968); Miksch (1969); Vlahakis (1972) and Vlahakis and Barduhn 
(1974), have been performed in which ice crystals grew in a flowing 
supercooled liquid. 

The results of these experiments relevant to this study may be 
briefly sunmarized. It was found that the growth rate of the ice 
crystals, as expected, increased with increasing supercooling and in- 
creasing velocity of flow. 

In the static case, for crystals growing freely in the fluid, the 
general relation for growth along the 'a'-axis, valid for supercool- 
ings less than 20°C, is given by 

da/dt = C-|(AT) n in cm s"^ . 

AT is the supercooling in °C; C-| is a function of the properties of 
the liquid, being equal to about 0.1 for pure water and decreasing with 
increasing solute concentration; and n is a constant equal to about 
2.0. The growth rate along the 'c'-axis is about 100 times less than 
growth along the 'a'-axis for supercoolings of about 3°C and becomes 
even less at smaller supercoolings. However, at large supercoolings 


24 


of about 15°C the growth rates along both axes are comparable. Hallett 
(1964) suggested that this marked change in growth along the 'c'-axis 
with supercooling might be due to a transition from a growth process 
controlled by concentration of imperfections to one controlled by 
surface nucleation. 

In the ventilated case there is enhanced growth in the upstream 

direction and suppressed growth in the downstream direction. This is 

the result of latent heat being carried downstream from the crystal 

tip and subsequently reducing the supercooling and hence the growth 

rate in the downstream region. Thus in the ventilated case, at vel- 

-1 -1 

ocities as low as 0.042 cm s and as high as 46 cm s , the growth 
rate along the 'a '-axis in cm s”^ is given by 

da/dt = ^ V 1/2 AT 3/2 

where A is a function of the known properties of pure water and ice 
and B is a function of the solution, V is the ventilation velocity in 
cm s' 1 and AT is the supercooling in °C. There is, however, an anomaly. 
That is, the growth rate versus solute concentration has a maximum 
rather than simply decreasing as the solute concentration increases. 

This anomaly is located in the dilute concentration range and the 
effect is not solute specific, but occurs with both electrolytes 
and nonelectrolytes. Furthermore, the maximum growth rate is more 
pronounced with increased supercooling. This anomaly remains to be 
explained. 

Some observations of Hallett (1964) of growth in static super- 
cooled water are also interesting. He found that with increased super- 
cooling both the separation and width of the secondary dendrite arms 
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decreased. He also observed that the angle of the growth front of an 
individual dendrite, determined by the relative rates of growth along 
primary and secondary <1 1 20> branches, increased with decreasing tem- 
perature from an angle of about 60° at -1°C to a maximum of 120° at 
-8°C, so that the limiting shape of the growth front became hexagonal. 
This limiting form occurred when the growth rates of the primary and 
secondary branches became equal and the secondary branches were 
nucleated very near the tip of the primary branch. 

D. Some Other Parameters Affecting Ice Crystal G rowth and Habit 

Several experiments in the past have been performed in an attempt 
to determine the variation in crystal growth due, not only to the 
temperature, ambient supersaturation and ventilation velocity, but 
also to such varied parameters as reduced or increased ambient pres- 
sure, the molecular weight and thus the thermal and vapor diffusivity 
of the carrier gas, contaminants in the carrier gas, an applied 
electric field and even dislocations in the crystals themselves. The 
observed variation of crystal habit associated with each of these 
parameters gives information which is vital not only in distinguishing 
ambient from surface kinetic effects but also in delineating the 
relative magnitudes of their importance under various growth conditions. 

Reduced and Increased Ambient Pressure and Variation of the Molecular 
Weig ht of the Carrier Gas 

Isono, Komabayasi and Ono (1957) and Isono (1958) studied the 
growth and habit of ice crystals produced by seeding a supercooled 
cloud of water droplets with silver iodide, in a cold chamber contain- 
ing different atmospheres. When crystals were grown in either hydro- 
gen at a pressure of one atmosphere or in air at low pressure, 20 to 
80 mm Hg, they observed nearly Isometric growth, c/a = 1, over the 
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entire temperature range -7°C to -16°C. However, at a pressure of one 
atmosphere the growth habits observed in carbon dioxide, nitrogen and 
oxygen were similar to those formed in air at normal pressure. 
Kobayashi (1958) using a convection chamber also investigated the 
effect of reduced air pressure on ice crystal habit and found compara- 
ble results. However, when Van den Heuvel and Mason (1959) performed 
experiments in a static thermal diffusion chamber using reduced pres- 
sures as low as 20 mm of Hg or using gases with differing molecular 
weights, they found variations in the rate of growth in accord with 
the thermal conductivity and the vapor diffusivity of the carrier gas, 
but they observed no real changes in the habits of the ice crystals. 
They suggested that traces of contaminants may have been present in 
the apparatuses of the Japanese workers. Later, it was suggested by 
Lamb (1970) that errors in the determination of the ambient super- 
saturation may have been the real source of conflict giving these 
differing experimental results. As Kobayashi acknowledged, it is very 
difficult to control the ambient supersaturation in an environmental 
gas through which water molecules readily diffuse. 

In an attempt to elucidate the contributions of the vapor dif- 
fusivity and the thermal conductivity of the carrier gas to the over- 
all ice crystal habit, Gonda and Komabayasi (1970), Gonda and 
Komabayasi (1971), Gonda (1976) and Gonda (1977) performed a series 
of experiments in which they grew small ice crystals, 20 to 50 pm 
diameter, in various atmospheres by seeding, with silver iodide smoke, 
a supercooled water droplet cloud contained in a small cold chamber. 

The atmospheres consisted of varying mixing ratios of helium and 
argon gases at both reduced and increased ambient pressures. By 
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changing the mixing ratios of these two inert gases and their total 
pressure it was possible to vary both the diffusion coefficient of 
water vapor and the thermal conductivity of the ambient atmosphere 
independently over a wide range of values. From early experiments 
they drew the following conclusions: 

1) The absolute growth rate of ice crystals increased with 

a decrease in the mean molecular weight of the carrier gas. 

2) The size of the central featureless portion of hexagonal 
plate crystals increased with a decrease in the mean 
molecular weight of the carrier gas. 

3) Hexagonal plates which exceeded a certain critical size 
developed branches at their six corners. This critical 
size increased as the mean molecular weight of the carrier 
gas decreased. 

Gonda's later experiments were conducted to study specifically the 
individual effects of the vapor diffusivity and the thermal conduct- 
ivity. Figure 2.10 shows diagrammatical ly the relative effect as 
determined by Gonda (1976) of the vapor diffusivity and the thermal 
conductivity on the habit of small ice crystals at temperatures of 
-7°C and -15°C respectively. From these results and the results of 
other experiments they concluded the following: 

1) Under conditions of constant vapor diffusivity (D = 0.77 
cm s ), skeletal structures and dendritic structures 
developed more readily as the thermal conductivity of the 
atmosphere was increased. 

2) Under conditions of constant thermal conductivity (K = 

34.4 X 10" 5 cal cm" 1 s" 1 °C _1 ) , skeletal structures and 
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dendritic structures developed less readily as the vapor 
diffusivity was increased. 

3) The dependence of the habit of small ice crystals on tempera- 
ture was most pronounced at high ambient pressures, but had 
a tendency to vanish at low ambient pressures. 
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Fig. 2.10 Schematic diagram of the habit and the growth feature of 
small ice crystals with respect to the vapor diffusivity 
and the thermal conductivity, (a) at -7°C; (b) at -15°C. 
Gonda (1976) 


They also suggested that the habit and growth features of ice crystals 
under identical ambient conditions are dependent upon the size of the 
crystals. Komabayasi (1970) performed a theoretical calculation and 
found that the effect of the diffusion coefficient on the shape of 
ice crystals is much larger when the crystal is small than when it is 
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large. Therefore, he reasoned that considerable discretion should.be 
exercised in comparing the habit and growth features of crystals larger 
than about 1 mm diameter, for example, those of Van den Heuvel and 
Mason (1959), with those of crystals less than 50 ym diameter. 

The experiments of Lamb (1970), which utilized a stainless steel 
substrate in an environment of pure water vapor over the temperature 
range 0°C to -20°C, were particularly valuable in providing informa- 
tion concerning surface kinetic effects. Since no vapor diffusion 
field existed over the crystals grown in an environment of pure water 
vapor and since the natural flow of heat was reversed by growing the 
ice crystals on a substrate of high thermal conductivity, surface 
kinetic effects were much more important than normal in controlling 
the crystal growth. Lamb found that in pure water vapor, at low 
pressure, ice crystals were generally regular hexagonal crystals when 
initially formed. However, as they increased in size both the sharp 
edges and the well defined crystallographic features disappeared. On 
the other hand, when significant amounts of either air or some other 
gas were present, the crystal edges and corners remained sharp even 
when the crystals grew quite large. Whether the gas present, at these 
low pressures, was air or helium made no visually detectable differ- 
ence. In either case, the surface structure of the crystals grown 
under these conditions was more complicated than that of crystals 
grown in pure water vapor. Linder conditions of crystal growth in pure 
water vapor. Lamb ascribed the sharp crystal features, that is, distinct 
crystallographic faces to 'inherent' surface kinetic effects controll- 
ing the growth mechanism. As these crystals grew larger and thicker, 
heat transfer to the substrate through the poorly heat conducting ice 
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crystal was diminished, so the ice crystal surface farthest from the 
substrate warmed up due to continued release of latent heat and hence 
its growth rate decreased to compensate for this reduced heat transfer. 
The thin lateral edges of the crystal, however, maintained good ther- 
mal contact with the substrate and growth in the lateral directions 
continued at an unreduced rate giving crystals a characteristic two- 
dimensional appearance. Lamb attributed the differences in crystal 
features when a carrier gas was present to the introduction of a 
vapor diffusion field rather than to the interactions of the carrier 
gas with the surface kinetics of the crystal . That is, when air was 
present the resistance to mass transfer controlled the growth rate, 
whereas in pure vapor the resistance to heat transfer was the limiting 
factor. 

Contaminants in the Carrier Gas 

Several experiments have shown that the normal ice crystal habits 
can be profoundly modified by the presence of certain organic vapors. 

Vonnegut (1948) observed that the addition of butyl alcohol of about 

_2 

10 mb partial pressure changed the ice crystals growing in a super- 
cooled cloud of droplets at -20°C from hexagonal plates to hexagonal 
columns. Schaefer (1949) found several vapors such as acetone, acetic 
acid, nitric acid, silicones and alcohols effective in modifying the 
crystal growth and habit. Hallett and Mason (1958a) and (1958b) using 
a static thermal diffusion chamber observed that small amounts of 
camphor vapor affected the crystal growth such that needle-like crystals 
appeared at all temperatures between 0°C and -40°C. They also in- 
vestigated in some detail the effect of iso-butyl alcohol on the ice 
crystal habit. Figure 2.11 from Hallett (1968) shows this effect as 



a function of iso-butyl alcohol concentration over a wide range of 
ambient temperatures. As can be noted from the figure, at high con- 
centrations of iso-butyl alcohol the habit at -15°C reverts to a 



Fig. 2.11 The effect of iso-butyl alcohol on ice crystal habit. 

Hallett (1968) 

malformed type of dendritic growth. Anderson, Sutkoff and Hallett 
(1969) showed how the presence of vapors of methyl 2-cyanacrylate 
resin suppressed dendritic growth along the 'a '-axis and induced fine 
fibers a few micrometers in diameter to grow along the 'c'-axis. 

It would appear that since foreign vapor concentrations, too 
small to affect the macroscopic mass and heat transfer processes are, 
nonetheless, effective in upsetting the normal crystal habits that 
they must somehow interfere with the surface kinetics. Furthermore, 
from these and other experiments, it appears that chemicals which are 
somewhat polar are most effective in modifying ice crystal habits. 
Indeed, if surface kinetic processes are, in fact, a growth limiting 
mechanism and the ice crystal surface is polar as it is believed to 
be, then these experimental results are reasonable, although by no 
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means completely understood. 

Effect of an A ppl i ed El ectr i c Fi el d 

Just as with variations of other macroscopic parameters, observa- 
tions of changes in ice crystal growth as a consequence of the appli- 
cation of an external electric field can be used to infer properties 
of the normal ice crystal growth processes. 

Although several studies, each producing essentially similar 
results, have been made of the influence of applied electric fields on 
ice crvstals growing in supercooled water clouds, the investigation of 
Evans (1973) was probably the most thorough. Other experiments, such 
as those by Bartlett, Van den Heuvel and Mason (1963) and Crowther 
(1972), have been performed in a static thermal diffusion chamber 
where the presence of cloud droplets was not necessary to sustain ice 
crystal growth conditions. In all three of the experiments cited, 
individual ice crystals were secured on a fiber in the presence of an 
electric field and could be observed throughout the growth process. 
Maybank and Barthakur (1967) employing a small cold cell also grew 
stationary ice crystals both with and without an applied electric 
field. They used, as fhe vapor source, a single water drop 1-2 mm in 
diameter suspended near the growing crystal. In each case the elec- 
tric field was produced by one of two methods. In the first method 
the electric field resulted from charging, to a high voltage potential, 
either positive or negative, the metal fiber on which the crystals 
grew. In the second method the crystals grew from an insulated fiber 
placed in a uniform electric field. In this case the electric field 
resulted from oppositely charging two parallel plates. 

In the various experiments employing a supercooled water cloud. 
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it was observed that at applied voltages to the metal fiber of less 
than or equal to ± 3000 V, electro-riming occurred. That is, long 
chains of crystals following the electric lines of force outwards 
from the electrodes grew very rapidly. These chains were largely com- 
posed of conglomerates of small crystals and very little could be 
deduced about the effects of the electric field on the individual ice 
crystal habits. However, at applied voltages exceeding ± 6000 V a 
silent discharge occurred creating a space charge in the vicinity of 
the electrode. Most of the cloud droplets entering this region acquired 
a charge of the same sign as the electrode and were consequently re- 
pelled, leaving a virtually droplet- free space around the wire. Sub- 
sequent growth of ice crystals on the wire electrode was from the 
vapor only. The crystals tended to grow as slender "needles" oriented 
in the direction of the field. The tendency for these "needles" to 
develop fragile branches was a strong function of the rate of growth 
which, in turn, was governed by the local vapor density. In fact, a 
critical local vapor density was necessary for the formation of 
branched "needles" at all. 

The experiments conducted in the absence of cloud droplets, for 
example, in the static thermal diffusion chamber, demonstrated that 
below a critical threshold value of the external electric field near 
450 volt cm - ^ , neither electro-riming nor rapid vapor growth called 
'electric growth' occurred. At field strengths somewhat larger than 
this threshold value, a time-lag existed between the application of 
the applied field and commencement of rapid growth. This time-lag 
decreased with increasing electric field strength. This time-lag was 
ascribed to the dependence of the magnitude of the local electric field 
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at the tip of the ice crystal on the shape of the growing tip. This 

local field, being larger the sharper the tip, varied from crystal to 

crystal and even with time for the same crystal. 'Electric growth' in 

-I 

large uniform external fields of about 500 volt cm or larger always 
took the form of "needles" with few if any side branches. When the 
field was removed side branches developed in the usual way. Near -4°C 
and -12°C the growth was generally along the ordinary <0001 > and <1 1 20> 
directions respectively. Sometimes, however, with proper orientation 
of the crystal to the external field, growth occurred along the <1 01 0> 
direction over this whole temperature range. Normal growth has never 
been observed along this direction. The linear rate of 'electric 
growth' was observed to be from 10 to 100 times faster at the same 
ambient temperatures and supersaturations than the linear rate of ordi- 
nary needles or dendrites. Although ice crystals grew equally well 
whether the metal fiber on which they were situated was of either 
polarity, 'electric growth 1 did not occur at all in a 50 cycle s”"* 
alternating field even at peak voltages five times greater than those 
needed to initiate 'electric growth' in a steady potential. Another 
interesting observation was the fact that the 'electric' crystals 
fractured spontaneously with the expelled fragments being highly 
charged. New growth then resumed on the remaining crystal remnant. 
Sometimes, rather than fragmenting, the growing crystal ejected a vor- 
tex ring from its tip which was interpreted as consisting of small 
clouds of ions made visible by their subsequent condensation of water 
vapor. 

It is an understatement to say that the observations of ice 
crystal growth in an electric field are poorly understood. Although 
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the increased growth rate may be exclusively attributable to either 
ambient kinetics, as suggested by Evans (1973), or to surface kinetics, 
it is also possible that 'electric growth* is somehow a complex func- 
tion of both. 

Crystal Dislocations 

The importance of dislocations in facilitating ice, or for that 
matter, any other crystal growth, particularly at low ambient super- 
saturations where two-dimensional nucleation is ineffective, has long 
been recognized. Since it was known that surface irregularities or 
steps on the growing crystal faces make accommodation of free molecules 
easier than on a smooth face, a mechanism for the repeated initiation 
of these steps on each new crystal growth layer was needed to explain 
observed growth at low ambient supersaturations. F. C. Frank (1949) 
proposed the existence of internal spiral or screw dislocations. 

Wherever a screw dislocation emerged on a crystal face a step would be 
initiated and propagated to each new growth layer. In general, this 
particular type of dislocation is only one of several which might be 
present in a given crystal. Depending on both the ambient conditions 
and the history of the crystal , the actual concentration of each type 
of dislocation present may vary over wide ranges. Until quite 
recently, studies of these dislocations have been limited to such in- 
direct methods as etching or observing growth after scratching the 
surface or performing some other such mechanical or thermal disturbance. 
However, with the advent and widespread use of x-ray topographic tech- 
nology in recent years these observations are becoming more direct and 
quantitative. 
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In summary, it is evident that although our knowledge of ice 
crystal growth processes is still lacking, much has been learned in a 
relatively short period of time. 
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III. 

ICE CRYSTAL GROWTH RATE - THEORY AND EXPERIMENT 

It has already been stated that the growth of an ice crystal is 
determined by a large number of variables. This is understandable 
since a water molecule must pass through several steps in going from 
a state of free movement in the carrier gas to a fixed position in the 
crystal. That is, ft must be transported to the crystal where it is 
adsorbed on the crystal surface. In the process of adsorption it 
must release a portion of its latent heat. It must migrate across 
the crystal surface to a growth site, generally a step or a kink in a 
step. Finally, the remainder of the latent heat must be liberated at 
the growth site as the water molecule becomes incorporated into the 
crystal lattice. It should be noted that all these steps in the 
crystal growth process are reversible and if the reverse processes 
take place sublimation is said to occur. 

From a macroscopic viewpoint the growth from the vapor of an ice 
crystal in air primarily consists of the simultaneous transfer of mass 
and heat. Transfer of mass involves the diffusion of water vapor 
through the air to the crystal. Transfer of heat involves both the 
diffusion of heat from the crystal and a radiation balance between 
the crystal and its surroundings. Since both diffusion processes will 
be enhanced by forced and natural convection, the total diffusion 
process is a combination of molecular and convective diffusion. 

A. Molecular Diffusion 

The non-steady state flux equation for molecular diffusion of 
water molecules through air is given by 
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ov 2 c w - ac w /at, 

where C w is the concentration of the water molecules and D is the dif- 
fusion coefficient which is a function of the composition of the 
carrier gas and varies inversely with the carrier gas pressure. The 
term 9C w /9t gives the variation in concentration of the water molecules 
with time. For the steady state 9C w /9t = 0, and a determination of 
the diffusive flux of water molecules reduces to the solution of 
Laplace's equation. Unfortunately, the number concentration of water 
molecules is difficult to measure directly. However, the equilibrium 
partial pressure of the water vapor can be measured directly as a 
function of temperature and expressed, by means of the ideal gas law, 

_3 

as the vapor density, p, in g m . Figure 3.1 shows the vapor densi- 
ties at saturation over plane ice and water surfaces respectively as 
a function of temperature. The difference in these two saturated 
vapor densities is also plotted against the scale on the right hand 
side. Since p is directly proportional to the concentration of the 
water molecules, it is possible to write 

Dv 2 p = 9p/9t. 

If an ice crystal is introduced to a new environment, then after the 
diffusion field surrounding it has attained a steady state, 9p/9t = 0. 

In which case, V 2 p = 0, or the vapor density at any point in space is 
given by the Laplacian equation, subject to p satisfying the boundary 
conditions at the surface of the crystal and at infinity. Thus, an 
analogy, for equivalent boundary conditions, can be made between the 
vapor field around a crystal and the electrostatic potential field 
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Fig. 3.1 Vapor density of water and ice and vapor density difference 
between water and ice versus temperature. 

around a charged conductor with the same shape and size as the crystal. 
As the crystal grows, latent heat is released at the crystal surface 
in an amount directly proportional to the mass of vapor deposited. 

The flux of heat from the surface by conduction is given by 

♦heat ■ -^ T 

where K is the thermal conductivity of the substance through which the 
heat is conducted and is the gradient of temperature surrounding 
the crystal. The quantity which determines the rate of diffusion of 
temperature is called the thermometric conductivity or simply the 
thermal diffusivity and is defined as 

k = -V 

p o G p 
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where K is the thermal conductivity, p Q is the density of the substance 
through which the heat is conducted and Cp is the specific heat at 
constant pressure of the heat conducting substance. Hence, the thermal 
diffusivity determines the rate of cooling due to a given temperature 
distribution according to the equation 

kV 2 T = 3T/3t, 

which implicitly assumes the thermal conductivity is not a function 
of temperature. For a steady state to exist, the latent heat must be 
dissipated at the same rate at which it is released. Otherwise, the 
crystal surface will warm up and the vapor density over the crystal 
surface will increase. In the steady state 3T/3t = 0, in which case 
V T = 0, or the temperature at any point in space is given by the 
Laplacian equation, subject to T satisfying the boundary conditions 
at the surface of the crystal and at infinity. Figure 3.2 demonstrates 
the resulting vapor density and temperature fields surrounding a sta- 
tionary growing crystal. No attempt has been made to illustrate the 
effect of the crystal corners on either the temperature or the vapor 
density field. 

Using the analogy between the vapor field around a crystal and 
the electrostatic potential field around a charged conductor with the 
same shape and size as the crystal and employing the proper boundary 
conditions, i.e., V = V c at the charged conductor surface and V = V ro 

at infinity for the electrostatic potential field and equivalently for 
the vapor field around a growing crystal p = Pj ce ( T c ) at the crystal 
surface and p = p^CTj at infinity, the rate of increase in the mass, 
m, of a crystal with electrostatic capacity, C, is shown by Mason (1971) 
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to be 


dm/dt = 4 ttC D (p^fTj - P ice CT c )) 

where D is the diffusion coefficient of water vapor in air, p^Tj 
the vapor density at a large distance from the ice crystal, and P-j ce (T c ) 
the vapor density at the surface of the growing crystal . Since the 
latent heat released by the molecules deposited on the surface warms 
the crystal surface a few tenths of a degree Celsius above its 
surroundings, the vapor density, Pj ce (T c )> at the surface of a growing 
crystal is larger than the vapor density, P.j ce (Too)> at the surface of 
a non-growing crystal. Hence, the "true" vapor density excess Ap true , 
environmental minus the crystal surface vapor density, with respect 
to the growing ice surface is less than would occur if the crystal 
surface did not warm up. 

The basic equation for diffusional growth of ice in a spherically 
symmetric vapor density field can be written as 

dm/dt = A D dp/dr. (3.1) 

The latent heat released at the surface of a crystal which is growing 
at a rate dm/dt is 

dQ/dt = L dm/dt (3.2) 

where L is the latent heat. For steady state growth the left side of 
eq. (3.2) must be balanced by the rate of heat removal by conduction, 
i.e., radiation is assumed negligible. The thermal conduction equa- 
tion for spherical symmetry can be written 
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dQ/dt = -A K dT/dr. 


(3.3) 


By combining eqs. (3.1), (3.2), and (3.3), the equilibrium relation- 
ship between p and T can be given in the form of the differential 
equation 

dp/dT = - -jjj- (3.4) 


where, in general, K, D and L are functions of T. This differential 
equation defines the variations of p with T between ambient conditions 
and the boundary conditions at the surface of the growing crystal. 

That is, it defines the relative spacings between the surfaces of 
constant vapor density and the surfaces of constant temperature in 
Fig. 3.2. 

Utilizing the vapor diffusion — thermal conduction equilibrium of 
a growing crystal and using the vapor diffusion and thermal conduction 
equations in conjunction with the latent heat, the ideal gas law, and 
water vapor density relations, one can compute a theoretical crystal 
mass growth rate equation in terms of experimentally measurable quan- 
tities. This has been done and is given by Mason (1971) to be 

dm/dt = (3.5) 


where a = [p/p $ (T ) - 1 ] is the supersaturation of the environment rela- 
tive to ice, C is the electrostatic capacity of the crystal and A and 
B are given as 


A a L 

A KT 


LM 

FT 


B = 


RT 

DMP $ (T) • 
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L is the latent heat of sublimation, M is the molecular weight of ice, 
K is the thermal conductivity of air, R is the universal gas constant, 
T is the ambient temperature in °K, D is the diffusion coefficient 
of water vapor in air, and P S (T) is the saturation vapor pressure over 
a plane ice surface at the ambient temperature T. At constant air 
pressure the terms A and B are functions of temperature only. The A 
term is the heat term and denotes the conduction of heat from the 
growing crystal. The B term is the moisture term and gives the rela- 
tion for the transfer of mass to the crystal surface. Figure 3.3 is 
a plot of the magnitude of these two terms as a function of temperature 
at an ambient pressure of 1000 mb. At other pressures the B term would 
be given by P a /1000 of its value where P a is given in mb. Figure 3.3 
indicates that at cold ambient temperatures and high ambient pressures 
the transfer of mass and not the transfer of heat is the primary 
ambient kinetic effect controlling the ice crystal growth. 

B. Convective Diffusion 

If convection is taken into account then both the vapor density 
and temperature fields surrounding the crystal will be altered. That 
is, the gradients dp/dr and dT/dr will be steeper in the direction 
into the flow. A rough model, adapted from Mason (1953), which shows 
the effect of a ventilation velocity on the crystal mass growth rate 
is presented below. 

The solution of the non-steady state equation 

DV 2 p = 3p/3t 

for a sphere of radius r Q with the boundary conditions p = p^CT^) 
everywhere when t = 0 and the vapor density P-f ce ( T c 1 ) at the crystal 
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Fig. 3.3 Comparison of moisture and heat terms from the mass growth 
rate equation as a function of temperature. Uses values 
of K, D and L at -10°C, and at a pressure of 1000 mb. 



surface is constant for all values of t is given by 


pJ t J - P(r) r c 

pJTJ - p ice (T c 'T ' r 

where T 1 is the crystal surface temperature in a static environment. 
So 


1-erf 


r-r. 


2(Dt Y 


(9p/3r) = 

0 


(pJO 


p ice^ T c ^ 



_J 

(ttDT)^ 


Thus, eq. (3.1) becomes 

dm/dt = AD dp/dr = 4 ttCD (pjTj - P ice (T c ')) 



where the capacitance for the sphere is r Q . 


(dm/dt ) ventilated (dm/dt) static 


Thus , 


1 + 


irDt 


where the second term in the square brackets is a measure of the in- 
crease in mass due to ventilation. This term will be large only if 
(•nDt) 1 * « r Q . The quantity (irDt)^ may be defined as the thickness of 
the diffusion boundary layer and has the dimensions of a length. For 
a spherical ice particle of radius r Q within a velocity flow v, the 
time t for which any given parcel of air will remain in contact with 
it before being replaced by fresh air will be given approximately by 
t = 2r Q /v. Hence, the diffusion field will have approximately this 
time in which to spread out from the crystal surface. Therefore, the 
mass growth rate of a ventilated ice sphere will be increased by a 
ventilation factor 
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The condition for the concentration gradient near the crystal to be con- 
siderably augmented due to the presence of the ventilation velocity is 


or 


2irDr 0 /v « r Q 2 
vr Q » 2 ttD. 


Hence, from strictly a mass growth rate standpoint a ventilated crystal 
can be regarded as growing at rest in a steady-state field and eq. 

(3.5) can be applied to the mass growth rate without appreciable error 
provided vr Q /D< 2 tt. However, even in this case, the ventilation 
velocity may affect the way in which the acquired mass is distributed 
over the crystal surface and hence the detailed crystal shape. Since 
the Reynolds number Re is given by Re = — and the Schmidt number Sc 
is given by Sc = v/D where v is the kinematic viscosity of the air, 
their product is given by 


so the ventilation factor is 

f = 1 + -J- Sc 2 Re^. 

2/if 

At 0°C the Schmidt number is given numerically by Sc = .1346/. 226 = .60 
so the ventilation factor becomes 

f = 1 + 0.23 Re h . 


An analogous result is obtained if the temperature field rather than 
the vapor density field is considered. However, in that case, the 
change with time of heat, rather than mass, is considered and one 
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obtains 


(dQ/dt ) vent ^ i atec j - (dQ/d^static 



where the thermal diffusivity k has replaced the analogous vapor 
diffusivity coefficient D. Furthermore, since the Prandtl number Pr 
is given by Pr = v/k, the product of the Prandtl and Reynolds numbers 
is 


and the heat ventilation factor is then given by 


f u = 1 + -L- Pr h Re h 
H 2/n 


where to a first approximation, f = f H . 

Since the ventilation factor f was only "derived" for a sphere 
its application to other shapes is questionable. In particular, one 
might expect the coefficient of the Reynolds number to vary with the 
shape of the crystal and also with the location at which the flow 
separates from the crystal. Nonetheless, no other acceptable theoreti- 
cal method of determining the effect of ventilation is available so 
with the inclusion of the ventilation factor, eq. (3.5) becomes 


dm/dt = 


4TrCaf 
A+B * 


(3.6) 


Equation (3.6) can only be used to predict the mass growth rate if the 
electrostatic capacity C of the crystal is known. Houghton (1950) 
suggested that shapes of conductors of known capacity be used to approxi- 
mate the shapes of actual crystals. McDonald (1963) by experimentally 
measuring the electrostatic capacities of brass models of various ice 



crystal shapes found that complicated shapes gave only small deviations 
from theoretical values of capacitance for symmetrical shapes with a 
similar appearance. 

The capacitance for a circular disk is given by C = 2r/ir and the 

O 

mass of such a disk is given by m = p^irr h where r is the disk radius, 
P d is the bulk density of the disk substance and h is the thickness of 
the disk. In this case eq. (3.6) reduces to 

dm/dt = 2-rrrhpjdr/dt = 4Tr (^/ u ) af 
or 


dr/dt = ( 



of 

A+F 


(3.7) 


which predicts the linear growth rate for a circular disk of constant 
thickness h. Likewise, for other simple shapes of known capacitance 
one can obtain a theoretical prediction of the linear growth rate. 

Applying eq. (3.5) with C = l/4ir to an ice crystal in an ambient 
environment at water saturation gives the calculated mass growth rate 
in g s" 1 as a function of temperature. The result of such calculations 
for two different ambient pressures are shown in Fig. 3.4. To obtain 
the true growth for any other C, the values need only be multiplied 
by 4 ttC. Note that a maximum occurs at -14.25°C at 1000 mb and at 
-16.75°C at 500 mb. The maxima are primarily due to the fact that 
the difference between the saturated vapor densities over water and 
ice. Fig. 3.1, is a maximum near these temperatures. However, the 
pressure effect on the two curves, that is both the increase in 
absolute growth rate and the shift in the maximum growth rate to lower 
temperature as the pressure is decreased, is brought about by the in- 
crease in the diffusivity of water vapor with decreasing pressure. 
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Fig. 3.4 Growth rates as a function of temperature of an ice crystal 
with C = l/4ir in a water-saturated cloud at two pressures. 

The temperature for maximum rate of growth is indicated on 
each curve. Byers (1965) 

If surface kinetic effects were unimportant, one would expect to 
experimentally observe a relationship between mass growth rate and 
temperature similar to the theoretical prediction depicted in Fig. 3.4. 
However, such a simple relationship has not been observed. Hallett 
(1965) found, by growing ice crystals in a static thermal diffusion 
chamber at water saturation, two maxima to exist in the growth rate; 
one at -4°C and the other at -15°C with a minimum growth rate at -8°C. 
These are illustrated in Fig. 3.5. Fukuta (1969) grew ice crystals 
in a cloud of supercooled water droplets and also found a dual maxima 
in the growth rate. Fig. 3.6. However, he observed the maxima to 
occur at slightly colder temperatures than found by Hallett. Ryan, 
Wishart and Shaw (1976) also grew ice crystals in a cloud of super- 
cooled water droplets and found a dual maxima; one at -6°C and the 
other at -15°C. They also presented a diagram. Fig. 3.7, showing the 
individual linear crystal axial growth rates along both the 'a' and 'c 1 - 
axes. Note that the magnitude of the growth rate along the 'a '-axis 
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Fig. 3.5 Relative rate of mass increase of crystals growing at water 
saturation in a static environment. Hallett (1965) 



Fig. 3.6 Mass of ice crystals versus temperature for various times 
after seeding a supercooled water cloud. Fukuta (1969) 
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Fig. 3.7 Variation of estimated crystal axial growth rates with 
temperature. Ryan, Wishart and Shaw (1976) 

at temperatures near -15°C for small crystals, growth time less than 
130 seconds, and for large crystals, growth time greater than 130 
seconds, is quite different. This can be explained by inclusion of 
the ventilation factor. That is, the larger crystals fall faster than 
the smaller ones in accordance with the experimental results of 
Kajikawa (1972), for plane type crystals. Fig. 3.8. Terminal fall 
speeds of crystals of larger dimension are given by Nakaya and Terada, 
Jr. (1935). 

Thorpe and Mason (1966) tested the validity of eq. (3.6) for 
Reynolds numbers ranging from 10 to 200 and ambient temperatures from 
0°C to -20°C by measuring the rates of evaporation of ice spheres. 

They found good agreement with the theory if they used (0.94 + 0.29 Re^) 
as the ventilation factor. Reynolds (1952) and Mason (1953) both made 
experimental measurements of mass growth rates for small crystals less 
than 100 ym diameter <grown in supercooled clouds in small cold cham- 
bers. Their results agreed within experimental error with the mass 
growth rates predicted by theory for a static environment, eq. (3.5). 
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Fig. 3.8 Falling velocity versus dimension of plane type crystals. 

Kajikawa (1972) 

Isono et al . (1956), who seeded a supercooled fog outdoors, found 
mass growth rates at temperatures near -4°C which were nearly twice 
those predicted by theory for a static environment. The crystals 
grown in the outdoor experiment had longer available time for sustained 
growth than those in the cold chambers and were therefore larger being 
typically 300 to 400 pm in length and hence had significant ventila- 
tion velocities. This would suggest that ventilation, as expected, 
may play a significant role in increasing the crystal growth rate. 

Several numerical empirical models have been developed to predict 
the mass or linear growth rate of ice crystals within a cloud. Two 
basic approaches have been used. The first approach fits empirical 
equations to experimental data on the axial growth rates. It assumes 
that the crystal grows only along two axes and different relative 
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growth rates along the two axes result in different crystal habits. 

The models of Todd (1964), Hindman (1968) and Hindman and Johnson 
(1972) are examples of this approach. The other approach uses eq. 

(3.6) with the incorporation of parameters determined from field 
measurements and has been formulated by Jayaweera (1971) and Koenig 
(1971). 

It is important to recognize that although eqs. (3.5) and (3.6) 
give a theoretical prediction of the mass growth rate of crystals in 
a static or ventilated environment, respectively, they do not predict 
the way in which the acquired mass is distributed over the crystal 
surface and hence the detailed crystal shape. They do not explain why, 
within specific temperature ranges, the growth should take place pre- 
ferentially along the 'a'-axis and, at other temperatures, along the 
'c'-axis. A complete understanding of ice crystal growth requires not 
only a knowledge of the rate of mass increase but also the way in 
which the accumulated mass is distributed over the crystal surface. 
Thus, it is necessary to have an understanding of the growth mechanisms 
by which the crystal faces are propagated. 

From a molecular standpoint, the macroscopic growth equation 
represented by eqs. (3.5) and (3.6) has some glaring deficiencies. 
Foremost is the fact that it does not consider surface kinetics. That 
is, it assumes that all molecules arriving at the surface are incor- 
porated directly into the crystal lattice. It also assumes a constant 
vapor density or molecular concentration over the entire crystal sur- 
face. This would imply no surface migration of molecules. 

Sometime ago, Marshall and Langleben (1954) hypothesized that 
vapor density excess was the primary growth mechanism controlling the 
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crystal habit. This would be true if it were not for surface kinetic 
effects. However, surface kinetic effects are very important and their 
relative magnitude of importance is a function of temperature for each 
crystallographic face. The experiments of Lamb (1970) in which ice 
crystals were grown in an environment of pure water vapor clearly 
demonstrated this fact. His results, presented in Fig. 3.9, show the 



Temperature (C) 


Fig. 3.9 Experimental measurements of the linear growth rates of 

the basal and prism faces of ice as a function of tempera- 

ture. The ice crystals were situated in an environment of 
pure water vapor and the excess vapor pressure was constant 
at 10 pm of mercury. Lamb and Hobbs (1971) 

linear growth rate of both the basal and prism face as a function of 

temperature for a constant vapor density excess. These are called 
'inherent' growth rates since they are governed by surface kinetic 
effects and not by the diffusion of water molecules to the crystal 
surface. The exact nature of the surface kinetic effects responsible 
for the shape of the curves in Fig. 3.9 is not known. Bernal (1958) 
speculated that the remarkable variation of ice crystal habit with 
temperature might be dependent on a variation with temperature of the 
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relative surface free energies of the individual crystallographic 
faces leading to different relative rates of nucleation and growth on 
the basal and prism faces as a function of temperature. Mason et al . 
(1963) suggested that the surface kinetic effects may be primarily 
due to the temperature dependence of the effective collection distance 
X c associated with each step. Fletcher (1968) and later Lacmann and 
Stranski (1972) have proposed a quasi-liquid film of varying thickness 
over the crystal surface to account for the temperature dependence of 
the surface kinetic effects. 

The experimental evidence for surface migration is based on the 
finding of Hallett and others that the speed of propagation of small 
steps on a crystal surface is inversely proportional to the step 
height. Direct experimental evidence for a quasi-liquid film on ice 
is sparse. However, recent nuclear magnetic resonance (NMR) experi- 
ments on ice by Kvlividze, Kiselev and Ushakova (1970) and Kvlividze, 
Kiselev, Kurzaev and Ushakova (1974) may indicate such a film. In any 
case, although the precise mechanisms responsible for the surface 
kinetic effects are not fully understood, their effect can be observed 
under the proper ambient conditions as Lamb (1970) demonstrated. 

In general, the crystal growth rate may be limited by insufficient 
transfer of heat, insufficient availability of water molecules, i.e., 
mass transfer, or by surface kinetic effects. Figure 3.3 showed that 
heat transfer is often not the limiting mechanism, at least at high 
pressures and low temperatures. However, if good heat transfer does 
exist, such as for growth on a substrate, the surface temperature of 
the growing crystal is reduced and the true vapor density excess over 
the crystal surface is increased. Hence, a reduction in the crystal 
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surface temperature is roughly equivalent to increasing the ambient 
vapor density. In most cases, the growth rate is limited by both mass 
transfer and surface kinetic effects. Figure 3.10 demonstrates how 
the effective vapor density is determined by the necessity of equality 
between the rates of transfer and deposition. In this figure, p^ is 
the ambient vapor density, p c is the vapor density at the crystal 
surface and p e is an effective vapor density such that (p e -P c ) is the 
true vapor density excess over the crystal surface. Complete surface 
kinetics control corresponds to a nearly vertical mass transfer curve 
and/or a nearly horizontal surface kinetics curve. Conversely, mass 
transfer control corresponds to a nearly horizontal mass transfer 
curve and/or a nearly vertical surface kinetics curve. The introduc- 
tion of a ventilation velocity corresponds to making the mass transfer 
curve more vertical and thus bringing the effective vapor density p fi 
closer to the ambient value p . 
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Fig. 3.10 The growth rate is determined by mass transfer and by inter- 
facial incorporation kinetics, which must take place 
sequentially and at the same rate. Adapted from Wilcox 
0971) 
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In summary although much is known about the general aspects of 
ice crystal growth under various ambient conditions, it is evident 
that a complete understanding of ice crystal growth has not yet been 
attained even under controlled conditions in the laboratory. 
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IV. 

EXPERIMENTAL APPARATUS AND PROCEDURE 


A. Apparatus 

The dynamic thermal diffusion chamber is based on the same opera- 
ting principle as that of the static thermal diffusion chamber. Two 
parallel, horizontal thermally conducting flat plates are separated by 
a small distance, d, as in Fig. 4.1. The top and bottom plate tem- 
peratures are controlled independently and are given by T t and T b , 
respectively. The top plate is maintained warmer than the bottom 
plate to make the system convectively stable. The equilibrium vapor 
density at the top and bottom plate ice surfaces is. given by P e (T t ) 



Fig. 4.1 Parallel conducting plates with associated equilibrium 
temperature and vapor density profiles. Adapted from 
Lamb (1970) 


and P e (T b ). Both the temperature T and the density of water vapor 
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p (x) vary almost linearly with height between the plates. However, 
since the equilibrium vapor density p e (T) over either ice or water 
varies exponentially with temperature as was noted in Fig. 3.1, a 
supersaturation develops in the region between the plates. Unlike the 
static chambers, the dynamic chamber invokes a controlled airflow 
parallel to and along the length of the plates. In this case it is 
necessary to solve the non-steady state heat conduction and vapor 
transport equations 

kV 2 T = 3T/3t and DV 2 p y = 3p y /3t, 

where k is strictly the thermal diffusivity of moist air but may be 

replaced by the thermal diffusivity of dry air since the mole fraction 

-2 

of water vapor within the chamber is only on the order of 10 . T is 
the temperature, t is time, D is the diffusivity of water vapor in 
air and p y is the density of water vapor. The boundary conditions 
for the heat transport equation are given by: 

T = Tjj , for x = 0, all t 

T = T t , for x = d, all t 

T = T(x) , for t = 0, all x. 

Likewise, the boundary conditions for the vapor transport equation are 
given by: 


p v = p vb , for x = 0, all t 
p v = p vt » for x = d, all t 
P v = P v ( x )» for * = °» a11 x - 
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If the incoming sample of air at time t = 0 is assumed to have a con- 


stant initial temperature distribution, T(x) = T Q , and contains water 
vapor saturated at the temperature T 0 and having a density p v (x) = p vq , 
the solutions to the one-dimensional heat conduction and vapor trans- 
port equations are obtained by separation of variables to.be: 


T = 


v T b 


x + T h + l 
b n=l 


•2.(V T o> + 2 < T f T o! 


nr r 


nu 


and 


sin(^) exp ( - (^) 2 kt ) 


(4.1) 


p v = 


p vt“ p vb 


+ P 


vb 


oo 

+ l 


n=l 


" 2 ^ p vb~ p vo^ + 2 ^ p vt _p vo ^ 


n -i 


nfr 


nrr 


sin(^) exp ( -(^ ) Dt ) 


(4.2) 


The time dependent exponential response term in each of these 
equations gives their respective time constants (for n = 1) 

j2 


h 2 

\ = 

T -. 2 , 


and A, 


it 2 D 


Tt"k M v 

After a sufficiently long time, that is several time constants, the 
exponential terms in both these non-steady state solutions go to zero 
and the steady-state linear temperature and vapor profiles are given 

by 

p vt” p vb 


T = 


V T b 


x + T b and p v = 


x + 


p vb' 


These are the same solutions obtained by solving directly the one- 
dimensional Laplacian equations 
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9 2 T/9x 2 = 0 and 9 2 p y /9x 2 = 0. 


The supersaturation at any height x above the bottom plate is 
given by 


Pw(x)-P p [T(x)] 


, where p v (x) = 


P e CT(d)]-p e CT(0)] 


x + P a CT(0)]. 


The supersaturation may be expressed relative to either ice or water 
depending on which equilibrium vapor density is used. In either case 
the supersaturation is often multiplied by 100 and expressed in per- 
cent. For a given set of top and bottom plate temperatures the super- 


saturation relative to ice is always larger than that relative to 


water. In fact, the supersaturation relative to ice in the thermal 
diffusion chamber will never be negative in the steady-state unless 
heat enters the system through the sidewalls. On the other hand, it 
is not unusual for the supersaturation relative to water to be nega- 
tive, that is, undersaturated with respect to water in which case 
water droplets evaporate. This means that ice crystals will grow in 
an ambient environment in which water droplets evaporate. The situation 


of undersaturation with respect to water but supersaturation with re- 
spect to ice is common in natural clouds and is called the Bergeron 


process. 

The air flowing between the horizontal plates of the dynamic ther- 
mal diffusion chamber is required to approach vapor and thermal equili- 


brium with the steady-state moisture and temperature profiles before 
reaching the end of the chamber where the ice crystals grow from a thin 
glass fiber suspended vertically between the two plates as in Fig. 4.2. 
The time needed for the air to approach these equilibrium values to 
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within 98 percent is given by 4A seconds where X is the time constant 
for the process under consideration. The minimum chamber length neces- 
sary is given by L = vt where v is the maximum air velocity and t = 4X. 



I 1.0 cm 1 


Fig. 4.2 Ice crystals grow from a fiber suspended vertically between 
two horizontal parallel plates. 


Thus, for vapor or thermal equilibrium to occur the chamber length 
must be given by 


L = 4vX - 


tFd 


for vapor equilibrium 


and 


L = 4v^ 


4vd 

ir 2 k 


for thermal equilibrium 


where d is the distance over which the diffusion must occur. The dif- 
fusivity of water vapor in air is a little larger than the thermal 
diffusivity of the air. This results in vapor equilibrium occurring 
before temperature equilibrium. Thus, as Fitzgerald (1970) has shown, 
the chamber will experience transient supersaturations higher than 
the steady-state values if the incoming air sample is saturated and 


64 


colder than the top plate temperature. This is of major importance if 

the chamber is employed in the measurement of nuclei spectra but is 

of little consequence in the present application. 

The maximum distance over which the diffusion of moisture and 

heat must operate is the plate separation, d. For air, k = D = .2 
? -1 

cm s and the necessary chamber length is given in terms of the 
plate separation and flow velocity as 

L = 2.03 vd 2 . (4.3) 

Thus the length of the chamber necessary for air with a given velocity 
to reach equilibrium must increase as the square of the plate separa- 
tion. If both the top and bottom plates are coated with ice, one would 
expect vapor diffusion from each plate and the actual diffusion distance 
to be somewhat less than the total plate separation. In this case the 
length calculated from eq. (4.3) would be an overestimate. 

The Reynolds number Re, defined as the ratio of inertial to 
viscous forces, is given for flow within closed tubes as 


o *i 

where r^ is the hydraulic radius and v = .125 cnr s“ is the kinematic 
viscosity of air. The hydraulic radius is used for non-circular pipes 
and channels and is defined as the ratio of the cross-sectional flow 
area to the wetted perimeter. If the width of my rectangular dynamic 
chamber is denoted by w then the hydraulic radius is given by 

wd 

r H = 2Tw+'dy 

and the Reynolds number becomes 
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Re 


_ 2v_ (wd) 
v Cw+d) 


(4.4) 


If the Reynolds number is less than about 2000 the flow in a smooth 
pipe will generally be laminar, Kay (1963). However, if the Reynolds 
number exceeds 2000 or if the surface is very rough the flow will 
become turbulent. 

My dynamic thermal diffusion chamber has a plate separation of 
2.5 cm and a horizontal width of 9.0 cm. This gives a dimensionless 
Reynolds number 


Re = 31.3 v 

where v is the airflow in cm s - ^ . For smooth tubes this implies 
laminar flow for velocities as large as 60 cm s"^ . These values of 
plate separation and horizontal width are both compromises. This 
large a separation is needed to insure that there is sufficient space 
for relatively large crystals, i.e., approximately 1 cm in length, to 
grow with a minimum of inhomogeneities in temperature, supersaturation 
and ventilation velocity. A greater width is preferred over the 
present value to minimize sidewall effects. However, if the width were 
any greater than the present value, the distance from the optical 
viewing window to the crystal would be greater than the focal length 
of our stereo microscope and would result in the sacrifice of a great 
deal of optical resolution. 

Figure 4.3 is a cutaway diagram of my dynamic thermal diffusion 
chamber as viewed from above and Fig. 4.4 is a photograph of the 
chamber as viewed from above, beyond the end. The chamber was generally 
operated as a closed system with the air being continuously recirculated. 
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|-« I METER 


C 

MOISTURE CONDITIONER 



Fig. 4.3 Diagram of the dynamic thermal diffusion chamber. (Top View) 

However, it has the capability of being readily separated at point F, 
Fig. 4.3, and used as an open system. The entire chamber lies in a 
nearly horizontal plane to suppress natural convection around the loop. 
The air drive indicated by point A, Fig. 4.3, consists of a Dayton No. 
2C781 shaded pole blower fan driven directly by a servo motor. The 
28 volt d.c. servo motor system is speed controlled with a 10 turn 
potentiometer and has a maximum speed of 2300 RPM. The regulated d.c. 
current to the servo motor is indicated in microamperes on a meter 
with a 2 3/4" diameter face marked off in 50 divisions. A switch 
varies the full scale deflection of the meter to any one of the values 
50, 100, 150, 200, 500 or 1000 microamperes. Since the RPM of the 
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servo motor and, in turn, the airflow in the chamber is directly pro- 
portional to this regulated current, these meter readings are more 
than adequate to provide a reproducible velocity at each set point. 

For a given meter setting the servo motor system maintains a constant 
ventilation velocity in the working section of the chamber. Further- 
more, the ventilation velocity resulting from this air drive system 
can be varied in infinitesimal increments over the entire range of 
interest, 0.1 cm s’^ < v < 1.0 m s“^. A 1/8 HP Bodine motor capable 
of 5000 RPM and equipped with a Minarik variable speed control is 
connected to the blower fan drive shaft via a V-belt and pulley. This 
motor, when employed, produces ventilation velocities of several m s”^ 
and is used in conjunction with a 2.5 m long scraper stick to blow the 
chamber clear of frost obstructions. 

The components of the chamber between points A and B and points 
C and D, Fig. 4.3, are connected with 3" diameter stainless steel pipe 
jointed with aluminum couplings fitted with Viton 0-rings. 

Point B, Fig. 4.3, indicates a butterfly valve used to close off 
the normal flow through the circular pipe and force the chamber air 
through an air filter. The air filter system, which removes both the 
cloud condensation nuclei (CCN) and the Aitken nuclei, consists of a 
Dyna-Vac Viton diaphragm pump and a circular cannister 22 cm in dia- 
meter and 23 cm in length packed with sterile absorbent cotton. The 
diaphragm pump provides an almost constant volume airflow rate of 6.4 
liter min - ^ even across such large pressure drops as the 10.8 inches 
of water pressure drop across the circular cannister. At this flow 
rate the entire volume of air in the chamber is filtered every 6 or 7 


minutes. 
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The moisture conditioner, indicated by point C, Fig. 4.3, is an ice 
labyrinth which preconditions the air before it enters the working 
section of the chamber. A large stainless steel tray 79 cm wide, 180 
cm long and 3% cm deep, containing cooling coils, has 30 liters of 
distilled water poured into it making the depth of water about 2.5 cm 
and covering all the cooling coils. Another rectangular stainless 
steel tray of just slightly smaller dimensions is inverted and fitted 
inside the first tray. Fins protrude from the top tray into the water 
nearly the length of the trays. When the water is frozen, an air 
tight ice labyrinth is formed with stainless steel sides and top and 
an ice bottom. The air path through the moisture conditioner is 
over 5.0 m in total length and approximately 22 cm wide by 1 cm high. 
This moisture conditioner was designed to precondition the air for a 
maximum ventilation velocity of 50 cm s”^ in the working section of 
the chamber. This corresponds to a volume flow rate of 1.25 liter s”^ 
or a linear velocity of 95 cm s - ^ in the ice labyrinth when Poiseuille 
type flow, that is a maximum velocity 3/2 times the mean, is assumed. 

The time constant for molecular diffusion in this case is given by 

^P v =0.5 seconds. 

The chosen labyrinth dimensions allow over five time constants for the 
air to approach the equilibrium vapor density of the ice surface even 
for ventilation velocities in the working section of 50 cm s~^ . Since 
the ventilation velocities employed in this study are less than half 
this value, the air remains in the ice labyrinth several time constants 
longer than necessary to approach 99% equilibrium. The necessary ice 
labyrinth area could have been reduced to about one- fourth the value 
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used if an ice surface had been employed on both the top and bottom. 
However, it would be very difficult to handle a top tray of the neces- 
sary dimensions if it were full of ice. Also since ice is employed 
only on the bottom, the top can always remain intact except for clean- 
ing and the water can be added directly to the bottom tray and drained 
through a stop cock after each experiment. Likewise, if the separation 
between the ice surface and the top tray were much smaller than 1 cm 
the time constant for diffusion would be greatly reduced, but it 
would be critical how level the tray is, how much warping is present 
in the top tray, how much water is added for each experiment and how 
the water freezes. 

A 36 gauge copper-constantan thermocouple was located just beneath 
the surface of the ice in the labyrinth. This double junction thermo- 
couple was referenced to an ice-water bath mixture and the millivolt 
output was read on a Data Technology Corp. model 350 digital volt meter. 
Counterflow was used in the heat exchanging coils to minimize the tem- 
perature gradients in the ice both across the labyrinth and down its 
length. This resulted in maximum temperature gradients across the ice 
surface of about 0.5°C. For the present study the temperature of the 
ice in the moisture conditioner was not critical but was generally main- 
tained intermediate between the top and bottom plate temperatures of the 
working section. If the ice labyrinth temperature were warmer than the 
top plate temperature of the working section, the air entering the work- 
ing section would be moist and would condense on the chamber walls. If 
the temperature of the ice labyrinth were very cold, the air entering the 
working section would be drier than necessary and the time for vapor 
equilibrium to occur in the working section would be a maximum. 
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The temperature conditioner, indicated by point D, Fig. 4.3, is 
a radiator 25 cm long, 25 cm wide and 10 cm high. Since the airflow 
must make a change in direction of 90° when it enters the temperature 
conditioner, three metal sheets were curved from different points on 
the radiator to the circular pipe where they were spaced to collect 
equal volumes of air. It was necessary, however, to insert a honey- 
comb screen in the circular pipe just ahead of the collection point of 
the three sheets to create turbulence and ensure that the air is always 
divided in the same proportions even at higher velocities when the 
parabolic flow in the circular pipe becomes more marked. There is 0.47 
cm distance between the vertical fins in the radiator and the fins 
are 0.03 cm thick. A honeycomb covers the exit of the radiator. The 
individual honeycombs are 0.67 cm across and the fins are 0.04 cm 
thick. The honeycomb material has a thickness of 1.4 cm. This arrange- 
ment gives a uniform flow at the exit of the temperature conditioner 
across its entire face independent of the air velocity. 

The temperature of the radiator is controlled by regulating both 
the temperature and flow of the coolant passing through it. This 
temperature is monitored with a thermocouple. In this study the tem- 
perature of the radiator was not critical and often the radiator was 
not cooled at all. Nonetheless, its presence gives the chamber an 
additional versatility. It is important to note that the temperature 
conditioner should never be operated colder than the moisture condition- 
er or frost will be deposited on the radiator. 

The effuser, constructed of stainless steel sheet material, is 
25 cm in width and 10 cm in height where it is connected to the tem- 
perature conditioner. Being 61 cm in length, it contracts to a width 
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of 9 cm and a height of 2.5 cm where it joins the working section of 
the chamber. At these low velocities the need for an effuser is 
marginal, so the angle of contraction is not critical. 

The working section consists of three portions: the pre-moisture 

section, the actual working section and the post-working section. 

Each of the three sections is 9 cm in width and 2.5 cm in height. 

Both the upper and lower horizontal plates are constructed of alumi- 
num and the sides of the chamber are of plexiglass which is a poor 
thermal conductor. The pre-moisture section is 64 cm in length. 

Along this length the top and bottom plate temperatures are maintained 
independently of one another and at nearly the same temperature as 
that of the corresponding plate in the actual working section. The 
top plate of the pre-moisture section is not coated with ice. The 
purpose of this section is to allow additional time for the incoming 
air to come to thermal equilibrium with the nearly linear temperature 
profile between the two horizontal plates. 

The actual working section is 122 cm long and has double walled 
plexiglass sidewalls with a double pane optical glass window opposite 
the growing crystals. The plexiglass sidewalls serve to insulate the 
top from the bottom and provide visible access to the entire section. 
The chamber sidewall not containing the optical glass window has nine 
1/2" test holes milled along its length with accompanying plexiglass 
plugs machined to fit. Both temperature and velocity probes are in- 
serted into the working section through holes drilled in the plugs. 

The probes can be moved with respect to the plugs and the plugs can 
be rotated so temperature and velocity profiles can be obtained over 
the entire horizontal and vertical expanse of the working section. A 


1/8" thick by 1" wide neoprene gasket cemented around the bottom of the 
top plates' perimeter acts as a mold for the ice surface. The top 
plate has a thermocouple imbedded just beneath the ice surface. When 
the top plate is operated at temperatures above freezing, a 1/8" thick 
stainless steel substrate 3" in width with .012" thick stainless steel 
wicking material fused to its exposed surface is screwed to the top 
plate fitting snugly into the mold normally reserved for the ice 
surface. The wick coated substrate which serves as the water vapor 
source consists of two lengths which butt together and are gravity 
fed distilled water from each end of the working section. The top 
plate has an oval hole 1.3 cm in width by 3.5 cm in length which is 
located 15 cm from the exit of the working section. This hole is 
fitted with a plug which holds the glass fiber on which the ice crystals 
grow. The plug has an ice or water surface respectively to match the 
rest of the top plate except for the press fitted 0.4 mm diameter tube 
through which the glass fiber is inserted. This plug seals with an 
0-ring and can be easily removed to facilitate the recovery of the 
ice crystals grown on the glass fiber. The bottom plate has two 
thermocouples attached to its unexposed surface. One thermocouple is 
interfaced with the control of the cooling unit, whereas the output 
of the other, along with the output of the thermocouple fixed to the 
top plate, is continuously recorded on a three channel chart recorder 
to ensure that no possible variations in temperature, during the 
period of crystal growth, go undetected. 

The post-working section is 38 cm in length and is not coated 
with ice. Its purpose is to insure that if any thermal gradient 
induced convection cells are present near the exit of the working 
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section they occur well downstream from the growing crystals. 

The cooling for the top plate along the entire working section 
is provided by a Neslab LT-9 circulator which uses isopropyl alcohol 
as the coolant. This circulator maintains the temperature of the 
alcohol stable to within ±0.05°C over the temperature range +10°C to 
-30°C. The cooling for the ice labyrinth, the temperature conditioner, 
and the bottom plate is provided by a built to order circulator from 
Cincinnati Sub-Zero. The coolant is also isopropyl alcohol with a 
temperature stability to within ±0.05°C, this time with an adjustable 
temperature range from 0°C to -50°C. A heat load consisting of a 
large radiator has been added to the system to increase the tempera- 
ture stability. 

B. Calibrations and Procedure 

An ice surface is formed on the top plate of the actual working 
section by lifting the top assembly of this section, inverting it, 
resting it on a special plywood panel, filling the neoprene mold with 
distilled water, freezing it and smoothing it level with the gasket 
height by rolling a heated 4" diameter solid brass cylinder down its 
length. The assembly is then re-inverted, lowered into place and 
sealed by securing the sidewalls to it with a gasket and screws. 
Additional sealing requirements, particularly around the glass fiber 
and sidewall plugs, are met with Apiezon sealing compound, a special 
inert claylike substance with an extremely low vapor pressure. To 
minimize possible temperature gradients both across and along the 
length of the top and bottom plates, the cooling fluid enters alter- 
nate cooling pipes via a manifold arrangement at each end of the 
plates. For both the top and bottom assembly the cooling pipes which 
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are in contact with one another and the aluminum plate are all cemented 
together with a high thermal conductivity epoxy, providing excellent 
heat transfer between the pipes and the aluminum plates, thus inducing 
a more uniform temperature distribution over the exposed plate surfaces. 
The actual working section, salvaged from an earlier prototype dynamic 
thermal diffusion chamber, has been modified but is the same basic unit 
described in more detail by Gamara (1972). 

The uniformity of temperature over both the top and bottom plate 
surfaces has been tested in a series of experiments. Gamara (1972) 
utilizing these same plates made several runs with thermistors imbedded 
in the ice surface on the top plate and found a uniform temperature 
(within ±.05°C) both across and along the length of the plate. Em- 
ploying thermocouples inserted through the plugs in the chamber side- 
walls, I have measured the temperature of the exposed surfaces of both 
the top and bottom plates to be uniform within ±.1°C except within 5 
cm distance of either the entrance or exit of the working section where 
the temperature was .3°C to .4°C warmer. With such temperature uni- 
formity over the two plates, a linear temperature distribution in the 
steady state should exist in the region between the upper and lower 
plates provided no flux of heat occurs through the sidewalls. To 
determine the validity of this postulate both horizontal and vertical 
temperature profiles were taken between the two plates for a variety 
of both top and bottom plate temperatures and ventilation velocities. 
Most of these profiles were taken through the test hole plug located 
15 cm from the exit, i.e., at the position where crystal growth on the 
fiber occurs. The vertical temperature profiles were taken with a 36 
gauge copper-constantan thermocouple shaped in the form of a crank. 
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_! , where it protruded through the sidewall plug. This shape 

facilitated temperature measurements over the entire vertical height 
and also provided a 4 cm length of thermocouple at the same height 
and temperature of the thermocouple tip thus reducing possible errors 
in the temperature measurement. Figure 4.5 shows a horizontal tem- 
perature profile taken 1.00 cm above the bottom plate with the top 



Fig. 4.5 A horizontal temperature profile taken 1.00 cm above the 
bottom plate in the dynamic thermal diffusion chamber 
(T t = -10.0°C; T b = -12.2°C) . 

plate temperature equal to -10.0°C, the bottom plate temperature equal 
to -1 2. 2°C, and the ventilation velocity zero. The dashed line at 
-11. 3°C indicates the temperature at this height assuming a vertical 
linear temperature profile between the two plates. The deviation of 
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the actual temperature from this value near the chamber sidewall oppo- 
site the plug is real. However, the actual deviation near the plug 
sidewall is not really as large as indicated because nearly all but 
the very junction of the thermocouple was outside the chamber and at 
a y/armer temperature resulting in a temperature measurement error. 

The general shape of this horizontal temperature profile and the fact 
that the temperature even at the horizontal center of the chamber was 
0.2°C warmer than it should be indicates that some heat is entering 
the chamber through the sidewalls. Vertical temperature profiles 
taken at zero velocity such as Figs. 4.6 and 4.7 also show that at any 
given height the actual temperature is tenths of a degree Celsius 
warmer than the temperature at that height assuming a linear tempera- 
ture profile. These two figures portray both the actual temperature 
profile, marked with X's, and the theoretical linear temperature pro- 
file, indicated by the solid line, in the region between the upper and 
lower plates. Figure 4.6 shows these vertical temperature profiles 
for a top plate temperature of -10.0°C and a bottom plate temperature 
of -10.3°C. Figure 4.7 gives the vertical temperature profiles for 
an upper plate temperature of -7.9°C and a lower plate temperature of 
-20.7°C. Note from these two figures that the larger deviations from 
the linear temperature profile occur near the middle and upper portions 
of the chamber. 

The top plate temperature is defined as the temperature measured 
by the thermocouple imbedded just beneath the surface of the ice in 
the middle of the top plate and 70 cm from the entrance of the working 
section. The temperature measured by the thermocouple imbedded on 
the backside of the bottom plate and 15 cm from the exit of the working 
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Fig. 4.6 A vertical temperature profile taken at the horizontal 
center of the dynamic thermal diffusion chamber (T^ = 
-10.0°C; T b = -10.3°C). 



Fig. 4.7 A vertical temperature profile taken at the horizontal 
center of the dynamic thermal diffusion chamber (T t = 
-7.9°C; T b = -20.7°C) . 
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section is always colder than the relevant bottom plate temperature. 

The temperature obtained by placing a thermocouple against the upper 
exposed surface of the bottom plate gives the relevant bottom plate 
temperature if no ice crystals grow there. However, at high ambient 
supersaturations, resulting from large temperature differences between 
the top and bottom plates, numerous crystals grow on the bottom plate 
making thermal contact between the thermocouple and the plate poor 
and also raising the local temperature tenths of a degree Celsius by 
the release of latent heat. Hence, under these conditions, as in Fig. 
4.7, the relevant bottom plate temperature lies between the imbedded 
and surface contact temperatures and is obtained by fitting the best 
straight line to the points composing the measured vertical temperature 
profile. 

Since the length of the chamber necessary for steady state transfer 
of heat between the upper and lower plates is a function of ventilation 
velocity and since heat does enter through the sidewalls, it is impor- 
tant to know how the temperature at any given vertical point in the 
chamber varies with the velocity. Figure 4.8 shows this variation of 
temperature with velocity at a height of 1.25 cm above the bottom 
plate and in the horizontal center of the chamber for two different 
distances from the entrance of the working section; 40 cm and 107 cm. 

The position 107 cm from the entrance of the working section is the 
location at which ice crystal growth on the glass fiber takes place. 

The experimental points in Fig. 4.8 were taken with the temperature 
conditioner at room temperature and the moisture conditioner at -4.5°C 
for a top plate temperature of -10.1°C and a bottom plate temperature 
of -10.1°C. Note from Fig. 4.8 that the temperature at the position 
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Fig. 4.8 Variation of temperature with ventilation velocity at a 

fixed height and distance from the entrance of the working 
section. (T t = -1 0 . 1°C; T b = -10. 1 °C ) . 

of the glass fiber, indicated with X's, becomes colder by tenths of a 
degree Celsius as the velocity increases from zero to about 9 or 10 
cm s" 1 and then warms quite rapidly at higher velocities. The lower- 
ing of the temperature at low velocities demonstrates the elimination 
of the wall effects in the horizontal center of the chamber when a 
flow is present. The rise in temperature at higher velocities indi- 
cates the maximum velocity for which diffusive equilibrium of heat 
occurs at that particular distance from the working section entrance. 
At the location of the glass fiber this maximum velocity is about 10 
cm s' 1 , being slightly higher if either the temperature conditioner is 
operated near the top plate temperature or the moisture conditioner is 
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operated at a temperature intermediate between the top and bottom plate 
temperatures. 

The ventilation velocities in my dynamic thermal diffusion chamber 
were in general determined with a Thermo-Systems, Inc. hot-wire ane- 
mometer system employing component models 1050, 1051-2 and 1057 with 
a model 1271-T1.5-6 tungsten hot wire probe. The probe is factory 
calibrated specifically for low velocities at an ambient temperature 
of -10°C and a sensor temperature of 100°C. This sensor temperature 
is a compromise. Higher sensor temperatures require that a smaller 
temperature correction be applied to readings taken at ambient tempera- 
tures other than -10°C. However, larger natural convection velocities 
at the hot wire, leading to bad measurements at low velocities, are 
associated with the higher sensor temperatures. For velocity measure- 
ments, the hot-wire anemometer probe was fixed in the chamber sidewall 
plug at the chamber position of the glass fiber. The hot wire, 1 mm 
in length, was always parallel to both the top and bottom plates of 
the chamber thus lying along an isotherm. The height of the hot wire 
above the bottom plate was adjustable and was determined with a stereo 
microscope equipped with a micrometer eyepiece. This same method was 
used in determining the position of the thermocouple for temperature 
measurements. Hence, the position of either the thermocouple or the 
hot wire probe could be determined in the vertical to within ±0.1 mm 
and in the horizontal to within ±1.0 mm. 

Appendix A shows that the transition from laminar to turbulent 
flow occurred around 40 cm s - ^ with a turbulence level of ±2.5 cm s - ^ . 
Since the turbulent regime occurred only at velocities much larger than 
those employed in this study, a consideration of laminar flow alone is 
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sufficient. 

Figures 4.9 and 4.10 illustrate the vertical and horizontal velo- 
city profiles respectively. The vertical profile was taken 4.5 cm 



Fig. 4.9 Vertical velocity profiles taken at the horizontal center of 
the dynamic thermal diffusion chamber. 

from the plug sidewall, that is, in the horizontal center of the cham- 
ber. The horizontal velocity profile was taken at a height of 1.36 cm 
above the bottom plate. Both profiles were taken with the top plate 
temperature at -10.2°C, the bottom plate temperature at -10.5°C and no 
ice on the bottom plate. As expected from the chamber dimensions, the 
horizontal flow did not become fully developed at as low velocities as 
flow in the vertical. For velocities greater than 10 cm s“^ the maxi- 
mum velocity occurred slightly above the vertical center of the chamber. 

The ventilation velocity for a particular RPM of the fan was a 
marked function of the static pressure. The static pressure was in 
turn a function of how much water was placed in the moisture conditioner 
and how this water froze. This necessitated making a velocity calibra- 
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Fig. 4.10 Horizontal velocity profiles taken 1.36 cm above the bottom 
plate in the dynamic thermal diffusion chamber. 

tion versus RPM of the fan, i.e., motor setting in microamperes, each 
time water was placed in the moisture conditioner and frozen. 

The velocity response of the dynamic chamber, as measured by the 
hot-wire anemometer probe located 107 cm from the working section 
entrance, is presented in Fig. 4.11 as a function of time. Since time 
in Fig. 4.11 increases from right to left, the peak on the right shows 
the velocity overshoot due to a surge of air when the fan was switched 
on. The line with the series of undulations superimposed on it indi- 
cates the magnitude of the steady state velocity, 23 cm s”^ , for this 
motor setting or RPM and the magnitude of the undulations represent 
the turbulence level at this velocity. The dip on the left indicates 
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•*- Time (minutes) 

Fig. 4.11 Velocity response when the fan is turned on and off. 


the response of the hot-wire anemometer when the fan was switched off 
but does not necessarily imply a momentary reverse flow in the chamber. 
In fact, such a reverse flow was never observed. Although the velocity 
response was almost immediate, the steady state velocity occurred about 
10 seconds after the servo motor was switched on at this RPM with an 
overshoot magnitude approximately 10 percent of its steady state value. 
At lower RPM the velocity response was slightly slower but the over- 
shoot magnitude was even a smaller percentage of its steady state value 
and the steady state velocity occurred sooner. 
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When an ice surface was present on the bottom plate, it was ob- 
served that for velocities greater than about 15 cm s“^ small ice 
particles approximately 5 to 10 ym diameter were lifted from the bottom 
plate and carried with the air stream. When the fan motor was switched 
off and on there was a delay of from 15 to 20 seconds before the ice 
particles were observed. The presence of the ice particles makes high 
velocity measurements perilous to the hot wire probe. However, an 
ev$n more serious consequence of their presence is the limitations they 
impose on high velocity applications of the dynamic chamber. 

The accuracy of the hot-wire anemometer velocity measurements was 
verified by two different approaches. In the first approach, the 
shutter speeds of two 35 mm Pentax cameras were electronically calibrated. 
It was determined that for shutter speeds slower than or equal to 1/125 
second the maximum variation in the shutter speed was less than two 
percent of the mean value of twelve experimental measurements at each 
shutter speed. For a given ventilation velocity or motor setting, 
photographs were made of small water droplets, approximately 10 urn dia- 
meter, which followed the airstream and were illuminated with a high 
intensity mercury vapor lamp. The photographs showed the droplet paths 
for a given shutter speed as straight lines. The distance traversed 
by a droplet for a given film exposure time was determined from the 
length of this straight line and divided by the shutter speed giving 
the velocity of the droplet and thus the airstream in cm s“^ . The 
second approach involved using the Dyna-Vac diaphragm pump with its 
nearly constant volume output independent of static pressure as tne air 
driving unit. In this case, the volume flow was measured with a factory 
calibrated airflow meter to be 6.4 ± .1 liter min - ^. The horizontal 
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and vertical velocity profiles in the working section for this flow 
were taken with the hot-wire anemometer. An eighth order curve fit 
was applied to these velocity profile points to give the three dimen- 
sional velocity profile in the working section. This velocity profile 
was then integrated by applying the divergence theorem to give the 
volume flow. Finally, the calculated volume flow was compared with 
the experimentally determined volume flow. Both these approaches gave 
results consistent with the factory calibration of the hot-wire 
anemometer. 

From all available data it was estimated that the hot-wire ane- 
mometer velocity measurements between 3 cm s - ^ and 25 cm s”^ had an 
absolute accuracy of ± .4 cm s - ^ and a relative accuracy of about ± 3 %. 
The percent variation of absolute accuracy was much larger at low 
velocities than at high velocities due to the interaction of natural 
convection from the heated hot wire with the forced ventilation. The 
velocity set-point repeatability was within ± 21 over the entire range 
0 cm s' 1 to 25 cm s - ^ . 

In computing supersaturation profiles in a thermal diffusion 
chamber, it has generally been assumed that the partial density of the 
water vapor is linear with height between the two plates of the chamber. 
Second order diffusion effects such as the mass transfer due to the 
temperature gradient (Soret effect), the energy or heat transfer due 
to the concentration gradient (Dufour effect), and the temperature 
dependence of the diffusion coefficients have been ignored. Only re- 
cently, however, Katz and Mirabel (1975) have solved numerically the 
heat and mass flux equations for a two-component system in which only 
one component, i.e., water, was diffusing without assuming linear 
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temperature, linear partial pressure or linear partial density pro- 
files between the plates. They included effects such as the temperature 
and composition dependence of the binary diffusion coefficient and of 
the mixture (air + water) thermal conductivity, the heat carried by 
the diffusing water molecules and both the Soret and Dufour effects. 
Using a top plate temperature of +20°C for a chamber at a total pres- 
sure of 1 atm, they found that the "best" solution to the complete heat 
and mass flux equations including all effects differed from that ob- 
tained by assuming linear temperature and linear partial pressure pro- 
files by only 1% over a range of temperature difference between the 
two plates of from 0°C to 20°C. They also found that assuming linear 
temperature and linear partial density profiles underestimated the 
value obtained from the "best" solution by about 12%. Hence, by 
fortuitous circumstances the use of a linear partial pressure profile 
between the plates gives a better representation of the ambient super- 
saturation in a thermal diffusion chamber with infinite parallel plates 
than does the use of a linear partial density profile. It will be 
shown that this is not necessarily true for a chamber with finite 
parallel plates. 

I developed two separate computer programs to calculate the 
supersaturation profiles in a thermal diffusion chamber with a plate 
separation, d. One program assumed a linear partial pressure pro- 
file, the other assumed a linear partial density profile and both 
assumed a linear temperature with height between the two plates. 

I used the Goff-Gratch formulation to calculate the equilibrium partial 
pressure of water vapor over both a plane surface of water and a plane 
surface of ice as a function of temperature. The density of water 
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vapor at saturation over an ice surface was obtained from the satura- 
tion vapor pressure according to the equation 


P 



Pei(T) 

W 


(4.5) 


2 

where P e ^(T) is the equilibrium vapor density in g m” over an ice sur- 
face at the temperature T, P el -(T) is the saturation vapor pressure in 
millibars over an ice surface at the temperature T, C y is the "compres- 
sibility factor" for water vapor, introduced to correct for the devia- 
tions of water from ideal gas laws, is the gas constant for water 
vapor, 4.6150 x 10^ erg g‘^ °K - ^ and T is the temperature of the vapor 
in °K. In a similar manner, the density of pure water vapor at satura- 
tion over a plane surface of water, p Q1 ,(T), was found by substituting 
the saturation vapor pressure over a water surface, P ew (T), for the 
saturation vapor pressure over an ice surface in eq. (4.5). It should 
be noted that tables in the Smithsonian Meteorological Tables and else- 
where of saturation vapor pressure over either water or ice were ob- 
tained by computing the Goff-Gratch equations at 0.5°C temperature in- 
crements and interpolating to get the values at 0.1 °C intervals. The 
use of a computer has greatly simplified these computations. There- 
fore, tables of both saturation vapor pressure and density of pure 
water vapor at saturation over both water and ice surfaces obtained by 
computing the Goff-Gratch equations at 0.1 °C temperature increments 
are included in Appendix B. Also included in this appendix is the 
percent supersaturation with respect to ice at which water saturation 
occurs as a function of temperature. 

Figures 4.12 and 4.13 illustrate the computer calculated super- 
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Fig. 4.12 Vertical supersaturation profile in a thermal diffusion 
chamber (Tt = -8.1°C; Tb = -20.6°C). Calculated 
assuming a linear vapor density with height. 



Fig. 4.13 Vertical supersaturation profile in a thermal diffusion 
chamber (Tt = -1.5°C; Tb = -20.5°C). Calculated assuming 
a linear vapor density with height. 



90 


saturation profiles with respect to both an ice and water surface for 
two different pairs of chamber plate temperatures. These supersatura- 
tion profiles were calculated assuming a linear partial density with 
height between the plates. Notice from these two figures that the 
maximum supersaturation with respect to ice is located below the ver- 
tical center of the chamber and is displaced even lower for larger 
temperature differences between the two plates. Also note that the 
maximum supersaturation with respect to water follows this same trend 
but occurs slightly higher in the chamber. 

In most calculations of the supersaturation profile inside a 
thermal diffusion chamber it is assumed that diffusion of water vapor 
and heat takes place between two parallel infinite plates. It is 
therefore very important to know just what aspect ratio, i.e., ratio 
of chamber width to chamber height, is large enough to permit this 
assumption. Twomey (1961) showed that for a static chamber an aspect 
ratio of 4 was the minimum for a saturation ratio of 1.01 with respect 
to water before wall effects became significant. Elliott (1971) calcu- 
lated that the aspect ratio should be no less than 5 and preferably 6 
or 7 or even larger. He further calculated that for a static chamber 
with a 1 cm plate separation and an aspect ratio of 5 one could safely 
use about 2.5 mm on either side of the centerline. For this same 
chamber with an aspect ratio of 7.5 one could use about 1 cm on either 
side of the centerline. These calculations were performed in the in- 
terest of using a thermal diffusion chamber to estimate the population 
of cloud condensation nuclei (CCN). However, this type of analysis is 
applicable to my dynamic thermal diffusion chamber where ice crystal 
growth on a vertical support located in the center of the chamber is 
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the primary interest. 

The most prominent wall effect is the introduction of heat through 
the sidewalls which manifests itself in raising the temperature and hence 
lowering the supersaturation. Furthermore, this introduction of heat 
through the sidewalls leads to the formation of a toroidal circulation 
with rising motion along the walls and sinking motion nearer the middle 
of the chamber. The introduction of heat through the sidewalls has 
its most pronounced effect on the ambient supersaturation at low values 
of supersaturation; that is, when the temperature difference, AT, be- 
tween the top and bottom plates is small. However, the effect becomes 
less significant as the forced ventilation velocity is increased from 
zero to the maximum velocity for which the incoming air continues to 
attain both thermal and vapor equilibrium. 

I performed a series of experiments, in addition to the tempera- 
ture measurements already described, to determine the extent of the 
toroidal circulation in my chamber. I made the chamber an open circuit 
by raising the motors and circulation fan. A 35 mm Pentax camera with 
an attached 200 mm lens and bellows was mounted on a tripod and 
positioned to view along an axis parallel and into the forced ventila- 
tion velocity. With a top plate temperature of -1.5°C and a bottom 
plate temperature of -20.5°C, giving a theoretical supersaturation 
greater than water saturation, i.e.. Fig. 4.13, CCN nucleated to form 
water droplets thus maintaining a supersaturation greater than or 
equal to water saturation. There was, nonetheless, an insufficient 
density of droplets to effectively photograph the toroidal circulation 
as a function of velocity for velocities in excess of 8.0 cm s . For 
higher ventilation velocities it was necessary to introduce droplets 
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generated by an Ultramist III nebulizer and having a 10 ym mean dia- 
meter, into the chamber through the port hole located 5 cm from the 
entrance of the working section. For velocities less than 17 cm s”^ 
the droplets were introduced through a 3/8" diameter polyurethane tube 
13 cm in length which was inserted horizontally into the port hole 
slightly above the vertical center of the chamber. Four holes in the 
tube, each approximately 750 ym diameter and located at distances 2.5, 
4.1, 6.4 and 8.3 cm respectively from the chamber sidewall opposite the 
port hole, allowed the droplets to be introduced parallel to the veloci- 
ty. For ventilation velocities greater than 17 cm s - ^ the holes in 

the polyurethane tube became plugged with ice in a short time. Hence, 

_1 

for velocities greater than 17 cm s the droplets were introduced 
perpendicular to the ventilation velocity without the aid of the poly- 
urethane tube. In both cases the droplets were introduced into the 
dynamic chamber airstream with the aid of a small controlled and mea- 
sured flow of air from a commercially filled compressed air bottle. 

This was effective and good photographs were obtained for velocities 
from 0 to 21 cm s - ^. For ventilation velocities greater than 8 cm s"^ 
the introduction of the droplets created turbulence. This was parti- 
cularly true for ventilation velocities greater than 17 cm s"^ because 
the droplets were introduced perpendicular to the main airstream. 

Figure 4.14(a), (b), (c) and (d) are photographs showing the toroidal 
circulation for ventilation velocities of 0, 6.7, 13.8 ± .9 and 20.5 
± 1.4 cm s" 1 respectively. 

Note from these photographs that a region apparently devoid of 
droplets existed in the upper portions of the chamber. This was a 
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consequence of both the maximum supersaturation occurring near the ver- 
tical center of the chamber and the finite growth times required for 
observable droplets. Also notice from these photographs that the in- 
troduction of a ventilation velocity forced the toroidal circulation 
closer to the chamber sidewall and away from the centerline of the 
chamber. This was the basis for the result obtained in Fig. 4.8 con- 
cerning temperature measurements. That is, at a given point near the 
centerline of the chamber the temperature lowered a couple tenths of 
a degree Celsius when a small velocity was introduced. The direction 
of the toroidal circulation also explains why larger deviations from 
the linear temperature profile occurred near the vertical middle and 
upper portions of the chamber as determined from Figs. 4.6 and 4.7. 

I determined experimentally that for appropriate top and bottom 
plate temperature pairs giving a band of maximum supersaturation slightly 
above water saturation near the vertical center of the chamber, CCN 
nucleated near the upper portions of this band and formed droplets 
which then settled slowly toward the bottom plate. As the droplets 
fell below this band of maximum supersaturation and into a region un- 
dersaturated with respect to water they evaporated and completely dis- 
appeared before reaching the bottom plate. Due to the finite time re- 
quired for either the growth or evaporation process to occur, the band 
of observable droplets was displaced slightly below the band of maxi- 
mum .supersaturation. The vertical width of this band was established 
by the temperature difference between the top and bottom plates, larger 
temperature differences corresponding to higher supersaturations and 
therefore wider bands of droplets. In an experimental attempt to 
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determine how well the computed supersaturation profiles for a thermal 
diffusion chamber with infinite horizontal plates matched the moisture 
profiles in my chamber with an aspect ratio of 3.6, I alternatively set 
the top or bottom plate temperature at a fixed value and incremented 
the temperature of the other plate by intervals of less than 0.1 °C so 
that the band of droplets became successively narrower until freely 
floating water droplets, i.e., less than 5 ym diameter, no longer were 
visible either with the naked eye and the aid of a laser or through our 
stereo microscope. I made the assumption that these pairs of top and 
bottom plate temperatures gave a maximum supersaturation in the chamber 
equal to water saturation. This assumption depends on the postulate 
that only a small number of nuclei activated at relative humidities 
less than water saturation grow to an observable size, i.e., 5 ym dia- 
meter. Table I is a tabulation of such experimentally determined 
temperature pairs. Using these temperature pairs and assuming a linear 

Table I Experimental Temperature Pairs 
at Which Water Saturation Occurs. 


Assumes Linear Assumes Linear 

Vapor Pressure Vapor Density 


T top(°C) 

T bottom (°C) 

Maximum % 
SS wrt Water 

Maximum % 
^wrt * ce 

Maximum % 
SS wrt Water 

Maximum % 

SS a. Ice 
wrt 

-2.6 

-11.6 

+0.14 

7.2 

-0.38 

6.5 

-5.9 

-17.3 

+0.23 

12.6 

-0.68 

11.4 

-9.0 

-22.0 

+0.51 

17.7 

-0.78 

16.1 

-10.9 

-24.4 

+0.11 

20.0 

-1.31 

18.1 


vapor pressure exists between the top and bottom plates, one calculates 
a maximum supersaturation in the chamber greater than water saturation 


by about 0.11 to 0.51 percent with respect to water as the top plate 
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temperature varies over the range -2.5°C to -11.0°C and the correspond- 
ing bottom plate temperature varies over the range -11.5°C to -24.5°C. 
Again using these temperature pairs and assuming a linear vapor density 
exists between the top and bottom plates, one calculates a maximum 
supersaturation in the chamber less than water saturation by about 0.4 
to 1.3 percent with respect to water over the same top and bottom plate 
temperature ranges as above. Hence, at least near water saturation, 
assuming a linear vapor pressure in calculating the supersaturation 
profile in my dynamic thermal diffusion chamber overestimates the actual 
supersaturation. On the other hand, assuming a linear vapor density in 
calculating the supersaturation profile underestimates the actual super- 
saturation. Therefore, near water saturation the ambient supersatura- 
tion at any height in my chamber can be calculated and specified with 
both an upper and lower limit. 

The effect of a temperature variation of ± 0.1°C at either the 
top or bottom plate on the maximum supersaturation is shown in Tables 
II and III for two quite different temperature pairs. Both the maximum 

Table II Effect of Temperature Variation of the Top or Bottom Plate 
on the Maximum Supersaturation, case (a). 


T bottom (°C) 

T top (°c) s \ N J 

-11 

.5 

-11 

.6 

-11 

.7 

-2.5 

+0.23 

(-0.30) 

7.2 

(6.5) 

+0.33 

(-0.21) 

7.3 

(6.6) 

+0.42 

(-0.13) 

7.5 

(6.8) 

-2.6 

+0.05 

(-0.46) 

7.0 

(6.3) 

+0.14 

(-0.38) 

7.2 

(6.5) 

+0.24 

(-0.30) 

7.3 

(6.6) 

-2.7 

-0.13 

(-0.63) 

6.8 

(6.2) 

-0.04 

(-0.55) 

7.0 

(6.4) 

+0.05 

(-0.47) 

7.2 

(6.5) 
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Table III Effect of Temperature Variation of the Top or Bottom Plate 
on the Maximum Supersaturation, case (b) . 


T bottom C°C) 

T top (°c'N x 

-24. 

,3 

-24. 

.4 

-24. 

.5 

-10.8 

+0.18 

(-1.24) 

19.9 

(18.1) 

+0.39 

(-1.06) 

20.3 

(18.4) 

+0.60 

(-0.87) 

20.6 

(18.7) 

-10.9 

-0.09 

(-1.49) 

19.6 

(17.8) 

+0.11 

(-1.31) 

20.0 

(18.1) 

+0.32 

(-1.13) 

20.3 

(18.4) 

-11.0 

-0.37 

(-1.75) 

19.3 

07.5) 

-0.17 

(-1.57) 

19.7 

(17.8) 

+0.04 

(-1.38) 

20.0 

(18.1) 


supersaturation with respect to water, i.e., the left hand figures in 
each column, and the maximum supersaturation with respect to ice, i.e., 
the right hand figures in each column, are presented in these tables. 

The supersaturations in parenthesis were calculated by assuming a linear 
vapor density between the two plates and the unparenthesized super- 
saturations were calculated by assuming a linear vapor pressure. Notice 
that the effect on the supersaturation is more pronounced if a tempera- 
ture variation occurs on the top plate than if a variation occurs on the 
bottom plate. 

At very low supersaturations, resulting from small temperature 
differences between the top and bottom plates, the effect of the side- 
walls on both the temperature and the supersaturation is a maximum and 
in this case assuming either a linear vapor pressure or a linear vapor 
density between the plates probably leads to an overestimation of the 
ambient supersaturation in the chamber. 

Since there is a difference in both the temperature and vapor 
density over an ice crystal which is growing as compared to one that 
is not growing due to the release of latent heat, the relevant growth 
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conditions are not the ambient or environmental conditions but the con- 
ditions which the crystal surface actually experiences. However, since 
these conditions are not readily determined, the ambient moisture, 
temperature and ventilation velocity conditions are generally specified. 
Nonetheless, by using conduction-diffusion equilibrium as discussed in 
Chapter III, eq. (3.4) 

dp/dT = (4.6) 

one can calculate a "true" vapor density excess for a crystal growing 
in a static environment. That is, in a plot of vapor density vs. 
temperature, eq. (4.6) corresponds to a curve such as that shown in 
Fig. 4.15 along with the saturation vapor density curve for ice. In 
Fig. 4.15 point B corresponds to conditions at the crystal surface. On 
the other hand, the ambient conditions correspond to some point A which 
is at a higher vapor density and lower temperature. In general, point 
A is defined by the experimental condition and it is necessary to find 
point B and thus the "true" vapor density excess Ap for growth in a 
static environment. 



Fig. 4.15 Temperature-vapor density conditions at the surface (point B) 
and in the environment (point A) of a growing crystal. From 
Rottner (1971) 
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Marshall and Langleben (1954) using a finite difference form of 
equation (4.6) assumed K, D and L were constant for the interval between 
points A and B. Since their treatment was equivalent to using a straight 
line between points A and B graphical determination of point B was sim- 
plified. 

In the computer program which I developed to calculate the super- 
saturation profile in a thermal diffusion chamber assuming linear 
partial density and linear temperature with height between the two 
plates, I included a section to calculate the "true" vapor density ex- 
cess for ice crystal growth in a static environment and also a calcula- 
tion of the maximum vapor density excess, that is, the difference be- 
tween the vapor densities at point A and point C in Fig. 4.15. In cal- 
culating the "true" vapor density excess in a static environment I in- 
cluded the temperature dependencies of both D and K but I used the value 
of the latent heat of sublimation at 0°C since it only varies 0.1% over 
the temperature range 0°C to -20°C. I also used an environmental 
atmospheric pressure of 850 mb which is the average pressure for Reno, 
Nevada. Figure 4.16 illustrates the resulting computer calculated pro- 
files of both "true" vapor density excess for a static environment and 
maximum vapor density excess for a top plate temperature of -8.9°C and 
a bottom plate temperature of -22.1°C. A sample of the actual computer 
output for this temperature pair is included in Appendix C. Also in- 
cluded in Appendix C is a diagram showing the calculated maximum vapor 
density excess which occurs in the vertical middle of the chamber as a 
function of the top plate and midpoint of the chamber temperatures. 

Note that the profiles of vapor density excess in Fig. 4.16 have 
their maxima near the vertical center of the chamber. These maxima. 
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Fig. 4.16 Calculated profile of maximum and "true" vapor density excess 
for a static environment (T t = -8.9°C; T b = -22.1°C). 

unlike the maxima in the supersaturation profiles. Figs. 4.12 and 4.13, 
occur even higher in the chamber for larger temperature differences be- 
tween the top and bottom plates. Thus, for large temperature differ- 
ences, i.e., greater than 20°C, between the plates the maxima for the 
vapor density excess, the supersaturation with respect to water, and the 
supersaturation with respect to ice all occur at somewhat different 
vertical positions in the chamber. The consequence of this is that 
when any one of these moisture designations is used to specify the 
ambient condition for growth occurring at various vertical positions in 
a thermal diffusion chamber with such large temperature differences be- 
tween the plates, a simple relationship does not exist to determine the 
other two moisture designations accurately. This is not a serious 
problem in this study since temperature differences between the plates 
rarely exceeds 20°C. However, this can be a very serious problem in 
static thermal diffusion chambers where much larger temperature differ- 
ences are often employed. 
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Two quite different definitions implicitly exist in the literature 
as to the vapor density excess for a given temperature at which water 
saturation occurs. One definition simply assumes, that for a given 
temperature, water saturation occurs at a vapor density excess equal to 
the difference in the saturation vapor density over a plane water sur- 
face and the saturation vapor density over a plane ice surface at that 
temperature. In this case, the vapor density excess for which water 
saturation occurs as a function of temperature is given by Fig. 3.1. 

This curve corresponds to the maximum vapor density excess as in Fig. 
4.16. The other definition assumes a conduction-diffusion equilibrium 
as in eq. (4.6) and Fig. 4.15. Thus, this definition takes into account 
the rise in temperature of the crystal surface due to the release of 
latent heat. This was the definition employed by Kobayashi in Fig. 2.5. 
This definition of water saturation corresponds to the "true" vapor 
density excess over a growing crystal in a static environment as in 
Fig. 4.16. Thus, the first definition specifies water saturation in 
terms of the ambient environment and the second definition specifies 
water saturation in terms of the environment with respect to a growing 
unventilated crystal. The first definition is used in this study. 

Since changing the temperature of the ice crystal surface alone 
corresponds to changing the "true" vapor density excess over the crystal 
surface, any process which changes the crystal surface temperature will 
affect its growth rate. In particular, if heat is conducted away from 
the crystal surface, through the crystal itself and into the substrate, 
i.e., the glass fiber, or carried away by a ventilation velocity, the 
crystal will grow faster. Contrariwise, if sufficient heat is supplied 
to the crystal surface from the substrate or through radiation from 
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either an illumination lamp or the top plate of the chamber, the 
crystal will grow more slowly. Since they are very important, the rela- 
tive magnitudes of heat transfer from the crystal surface by conduction 
to the substrate, conduction to the air, convection to the air, and 
radiation were calculated and are presented in Appendix D. Their 
relative magnitudes listed in decreasing order of importance were found 
to be; convection to the air, conduction to the air, conduction to the 
glass fiber and radiation. 

The temperature of the glass fiber was measured by placing a 36 
gauge copper-constantan thermocouple against it. When a ventilation 
velocity was introduced, the temperature at a given height of both the 
air and the glass fiber decreased by about 0.2°C in conformity with 
the results found earlier, i.e., Fig. 4.8. Even so, no detectable 
temperature difference, i.e., less than ± .05°C, was observed between 
the air and the glass fiber at the same height in the chamber regard- 
less of whether or not the glass fiber was in contact with the bottom 
plate or whether or not a ventilation velocity existed. However, when 
the glass fiber was in contact with the bottom plate and the plate 
temperatures were given by one of the temperature pairs in Table I, 
small water droplets condensed, after a long period of time, on the 
glass fiber near the vertical center of the chamber even though freely 
floating water droplets were not observed. This suggested that possibly 
the glass fiber was slightly colder, i.e., less than .05°C colder, than 
the ambient air at a given height. When the glass fiber was replaced 
with a copper wire of comparable diameter and one end of the copper 
wire was terminated one centimeter above the bottom plate but the other 
end extended through the top plate into the room air, the temperature 
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of the copper wire at a given height was tenths of a degree Celsius 
wanner than the ambient air at that height. This temperature difference 
decreased in magnitude when a ventilation velocity was introduced. When 
good thermal contact was made between the copper wire and the bottom 
plate, that is, the copper wire was bent so 1.25 cm of copper wire lay 
along the bottom plate, no detectable temperature difference existed 
between the air and the copper wire at the same height in the chamber 
either with or without the presence of a ventilation velocity. These 
results can be explained in terms of the difference in thermal conduc- 
tivity between copper and glass; the thermal conductivity of copper being 
nearly 500 times that of glass. Therefore, when ice crystal growth 
took place from a glass fiber along a horizontal isotherm, that is at 
an angle of 90° to the vertical fiber, an insignificant amount of heat, 
for all but short crystals at very low ambient supersaturations, was 
conducted from the crystal surface to the fiber. However, when growth 
took place at an angle from the horizontal, the crystal traversed a 
temperature gradient and a small amount of heat was conducted along the 
crystal. If the crystal angled upward, toward the top plate, the base 
of the crystal was colder than the crystal tip and heat was conducted 
away from the crystal tip to the glass fiber enhancing the crystal 
growth rate. On the other hand, if the crystal angled downward, toward 
the bottom plate, the base of the crystal was warmer than the ambient 
air at the crystal tip but, due to the release of latent heat, not 
necessarily warmer than the crystal tip. In this case, heat conduction 
along the crystal was less than for the case when the crystal angled 
upward but its direction was vague. 

The proximity of the ice crystals to one another on the glass 
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fiber sometimes resulted in an overlap of their individual vapor fields 
and a consequent reduction in their growth rates. This competition be- 
tween the growing crystals for the available water vapor was most pro- 
nounced in a static environment at low ambient supersaturation. In 
this case, the vapor density field surrounding a crystal may extend 
several crystal diameters from the crystal. However, at high ambient 
supersaturations the vapor fields typically extended one or two crystal 
diameters from the crystal and even this distance was reduced when a 
ventilation velocity was employed. 

The observed consequence of the superposition of these effects 
was that when two crystals with different base heights, the lower crys- 
tal growing upward and the higher crystal growing downward, approached 
one another, under conditions of low ambient supersaturation, the lower 
crystal with the colder base temperature grew at the expense of the 
upper crystal . 

Since certain organic vapors affect normal ice crystal habits, every 
precaution was taken in this study to exclude from the system possible 
sources of such foreign vapors. As an added precaution, a contamination 
check was made of both the laboratory air and the air from the working 
section of the dynamic chamber. This was accomplished by using a Dyna- 
Vac diaphragm pump to draw an air sample through a cold trap comprised 
of a special glass tube immersed in solid C0 2 (dry ice). The cold trap 
effectively condensed from the air the water vapor and any other gases 
present having a condensation temperature warmer than -78.5°C, the sub- 
limation temperature of solid C0 2 . The resulting solid sample from the 
cold trap was melted, well mixed and then injected into a Varian Aero- 
graph model 1740 gas-liquid chromatograph. The chromatograph column 
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used was a Porapak Q, 7' long by 1/8" O.D. stainless steel column with 
a 150/200 mesh. The carrier gas was nitrogen and a flame ionization 
detector was utilized. With this system the detectable limit was about 
one part unknown per million parts of the water sample. However, know- 
ing the mixing ratio of the sampled air enabled this figure to be 
expressed in terms of the air sample, giving a detectable limit of 
about ten parts per billion of the unknown to the air sample. Using 
this contamination check system only one foreign vapor was detected in 
either the laboratory or the dynamic chamber air. The one detectable 
vapor was identified as isopropyl alcohol, the coolant of the two cir- 
culators. However, even this vapor only had a concentration between 
0.02 and 0.2 parts per million (ppm) in the air and the concentration 
was the same in both the laboratory and dynamic chamber air samples. 
Therefore, possible effects of foreign vapors on the ice crystal growth 
rates were dismissed in this study. 

The normal procedure which was followed in operating the dynamic 
thermal diffusion chamber is given below. For each initial operation of 
the chamber, following a period of non-use, 30 liters of distilled water 
was placed in the moisture conditioner and frozen. Either an ice or 
water surface was employed on the chamber top plate depending on the 
range of ambient temperature and supersaturation desired. The top 
plate was secured in place with a gasket and screws. With both the top 
and bottom plate temperatures near -10°C the air velocity was measured 
with the hot-wire anemometer as a function of the motor setting in 
microamperes. If the ice on the top plate was then replenished once 
or twice a week, depending on how rapidly it sublimated away, the chamber 
could be operated continuously for a period of two to three weeks without 
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a change in the velocity calibration. For day to day operation the 
chamber was made an open circuit by diassembling it at point F in Fig. 
4.3 and raising the motors and circulation fan. Using the 2.5 m long 
scraper stick and employing the high speed Bodine motor, ice crystals 
growing on the bottom plate were removed from the chamber. The end of 
the chamber was re-assembled making it a closed system and the top and 
bottom plates were set at the desired temperatures. When temperature 
equilibrium was achieved for both plates, a vertical temperature pro- 
file was taken. For ambient supersaturations greater than water satura- 
tion the air was then passed through the nuclei filter for a period of 
from 30 minutes to 1 hour to insure that all the CCN were removed. If, 
following the filtration process, no droplets could be detected in the 
working section, a clean glass fiber was inserted into the chamber and 
the servo motor was set for the desired ventilation velocity. If drop- 
lets were detected, an air leak was indicated and steps were taken to 
eliminate it. Following the nucleation of ice crystals on the glass 
fiber, the fiber was rotated so the crystals to be observed and photo- 
graphed grew directly into the airstream. A thermocouple, introduced 
through a small hole in the sidewall plug just opposite the growing 
crystals and sealed with Apiezon compound, could be retracted during 
the run so the thermocouple bead was flush with the chamber inner side- 
wall, thus preventing the growth of ice crystals on its tip, or it could 
be extended to measure the temperature near the ice crystals growing 
from the glass fiber. Hence, accurate ambient temperature measurements 
were obtained during the period of crystal growth. 

Observations of the growing ice crystals were made using an M7- 
Wild Heerbrugg zoom stereo microscope. This microscope had a binocular 
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tube for viewing and a monocular tube for photomicroscopy. With the 
available choice of objectives: IX or 0.5X, and eye pieces; 10X micro- 

meter, 10X, 15X and 25X, this viewing system gave a wide range of mag- 
nification varying from about 3X to 80X. Ice crystal growth, as viewed 
through the stereo microscope, was recorded on film with either a 35 mm 
Pentax camera or a 16 mm Beaulieu cine-camera which was adapted to take 
time lapse photographs over the range one frame per fifteen minutes to 
64 frames per second. This range was more than adequate for all the 
crystals grown. Both the time interval between successive frames and 
its variability were determined by taking time lapse photographs of a 
stop watch. Even at the lowest time lapse rate the variability of the 
time interval between successive frames was negligible. 

In general, background lighting alone, supplied by an American 
Optical microscope lamp positioned on the opposite side of the chamber 
from the microscope, was used. The heat filtered background lighting 
was made uniform by passing it through a 1/8" thick frosted glass. 

This lighting arrangement eliminated possible radiational heating of 
the ice crystals due to the microscope lamp. When colored filters were 
placed between the frosted glass and the ice crystals and the ice 
crystals were illuminated directly with low intensity heat filtered 
light, spectacular photographs of the ice crystals against a colored 
background were obtained, as in Figs. 4.17 and 4.18. Simply changing 
the filter allowed any background color desired. 


* 

i 



108 



Fig. 4.18 Ice crystal dendrites growing against a red background. 
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V. 

RESULTS AND CONCLUSIONS 

A. Crystal Growth in a Droplet Free Environment 

Both the shape and the linear growth rate of ice crystals grown in 
the dynamic thermal diffusion chamber were found to be marked functions 
of the forced ventilation velocity. • This was particularly evident for 
growth under ambient conditions such that the introduction of a forced 
ventilation velocity resulted in a supersaturation dependent secondary 
habit transition, i.e., plate+dendrite or column-meedle. 

Predominant Growth Along the 'a* -Axis 

Figure 5.1 illustrates the difference in the shape of a crystal 
when grown along the 'a'-axis with and without a ventilation velocity. 

In both cases the crystal, indicated with the arrow, grew under ambient 
conditions of temperature near -15°C and at an ambient supersaturation 
near 22.5% with respect to a plane ice surface. Photograph 5.1 (a) 
shows the shape and length of the various crystals 9.0 minutes after 
the glass fiber was inserted into the chamber. During this time a ven- 
tilation velocity, v, flowed from left to right with a constant magni- 
tude of 3.0 cm s' 1 . The velocity was decreased from 3.0 cm s -1 to zero 
as shown in Fig. 5.1 (b) at time t = 21.0 minutes after the fiber was 
inserted. Figure 5.1 (c) shows the shape which the crystals assumed 
in the absence of a velocity. The plates grew at a slightly warmer 
temperature and lower supersaturation than that encountered by the 
spear-like, dorite crystal. The ventilation velocity was increased 
from 0 to 7.0 cm s" 1 at time t = 53.0 minutes, and 6.0 minutes later 
the crystals had resumed a dendritic shape as illustrated in Fig. 5.1 (d). 
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Glass fiber inserted at t = 0 minutes with v = 3.0 cm s - ^ 




Fig. 5.1 


A forced ventilation velocity enhances ice crystal growth 
into the airstream. (T-fc = -4.5°C; T 5 = -20.2°C). 
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These photographs show that a ventilation velocity enhances ice crystal 
growth into the airstream. 

The length of the crystal , indicated with the arrow in Fig. 5.1, 

is plotted as a function of time in Fig. 5.2 for the successive ventila- 

*1 

tion velocities of 3.0, 0, and 7.0 cm s . For a given velocity, the 
crystal length along the 'a'-axis as a function of time can be repre- 
sented with a straight line, whose slope gives the linear growth rate 
of the crystal. The linear 'a'-axis growth rate for this set of ambient 
temperature and supersaturation conditions was 2.2 pm s - ^ for a ventila- 
tion velocity of 3.0 cm s - ^ ; 0.72 ym s“^ for the static case; and 2.8 
ym s" 1 for a ventilation velocity of 7.0 cm s”^ . 

The supersaturation dependent secondary habit transition dendrite-* 
dorite occurred when the ventilation velocity was reduced from 7.0 to 
0 cm s" 1 at an ambient temperature near -14°C and an ambient super- 
saturation of 25% with respect to ice, i.e., a vapor density excess 

O 

equal to 0.39 g m . An example of this transition is given in Fig. 

5.3 (a), (b), (c), and (d). Also note from Fig. 5.3 (e) that the crystal 
resumed a dendritic structure when the ventilation velocity was return- 
ed to its original magnitude, i.e., 7.0 cm s”^ . Notice from these 
photographs that both the crystal shape and the linear growth rate were 
functions of the conditions of ambient temperature, ambient super- 
saturation, and forced ventilation velocity under which growth occurred. 
Secondary branches of dendrites, even though growing under the same 
ambient temperature conditions as a primary branch often grew broader 
or more plate-like because they "saw" a lower local supersaturation due 
to competition with other branches. This competition effect, however, 
was minimized when a forced ventilation velocity was applied. 



Crystal Length (mm) 


Fig. 5 



2 Ice crystal length as a function of time for successive 
ventilation velocities of 3.0, 0 and 7.0 cm s - l , under 
constant ambient conditions of temperature, -15.8°C, and 
vapor density excess, 0.279 g m - 3, at the height of the 
crystal base, 0.70 cm above the bottom plate. Angle of 
crystal from horizontal was +18°. (T^ = -4.5°C; T 5 = 

-20.2°C) . 
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v = 0*7.0 cm s" 1 at t = 20.0 minutes 

e) 1ft 
4 



v = 7.0 cm s 
t = 31 .0 minutes 


Fig. 5.3 (Continued). A forced ventilation velocity affects both 

ice crystal shape and linear growth rate. 


The task of separating the temperature and supersaturation effects 
at a fixed ventilation velocity was complicated by the fact that the 
crystals often nucleated and grew at an angle other than 90° to the 
glass fiber. At mid-chamber temperatures colder than or equal to 
about -10°C, ice crystals nucleated on the solid glass fiber without 
external influence. At warmer temperatures nucleation was induced on 
the fiber by immersing it in liquid nitrogen before inserting it into 
the chamber. In both cases, it was found that employing a forced venti- 
lation velocity immediately after the glass fiber was inserted into 
the chamber helped minimize the angle at which the dominant crystals 
grew from the horizontal, as well as inducing preferred growth in a 
direction parallel to the field of view and into the ventilation 
velocity. Once preferred growth occurred, further nucleation and growth 
of other crystals in the immediate neighborhood of a dominant crystal 
was suppressed due to competition for the available water vapor. Thus 
linear crystal growth rates were rarely assessed in this study if the 
crystal grew at an angle greater than ± 20° to the horizontal. 

It was also important that the crystals grew parallel to the 
field of view of the microscope. If the crystals grew at an angle <j> 
from that field of view, then both the crystal length measurements and 
the subsequent linear crystal growth rates would be reduced by an 
amount (1-cos cj>). For example, for an angle <J> = 20° from this field 
of view the crystal growth rates would be reduced by as much as 6%. 
However, since at high magnification the depth of field of the stereo 
microscope was small, this could easily be detected and then corrected. 
That is, if at high magnification the crystal was partially out of 
focus, the glass fiber was rotated to ensure that the crystal grew 



116 


parallel to the field of view of the microscope and pointed directly 
into the ventilation velocity. 

Fig. 5.3 illustrates a phenomena which was commonly observed 
during the growth of a dendrite under changing conditions of ventilation 
velocity. That is, when a velocity was removed, subsequent growth 
along the primary branch of the dendrite often assumed a slightly dif- 
ferent orientation than before, as in Fig. 5.3 (d). However, when a 
velocity of the same magnitude as before was again employed, the pri- 
mary branch of the dendrite resumed its former orientation, as in Fig. 
5.3 (e). This may be explained in terms of the local vapor field sur- 
rounding the crystal being different in the two cases. The crystal in 
both cases attempted to grow along the maximum vapor field gradient. 

This is similar, although somewhat more pronounced due to the larger 
change in the local supersaturation induced by the introduction of a 
forced ventilation velocity, to an effect sometimes observed on the 
ribs of sector plates. That is, the ribs of sector plates also often 
grow curved in the direction of maximum water vapor. Frank (1974) 
termed this phenomena lacunary growth and Yamashita (1976) described it 
in more detail. In both the case of the dendrite arid the sector plate, 
this curved growth is due to the fact that nucleation and growth occurs 
faster on one of the leading prism faces than it does on the adjacent 
one. 

The photographs in Fig. 5.4 demonstrate how a ventilation velocity 
enhances growth along the leading edges of a crystal, resulting in this 
case in the transition dorite^dendrite. This set of photographs also 
shows in detail both the manner and time scale in which this transition 
proceeds. It fits into the time sequence of Fig. 5.3 between photographs 
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v = CH-7.0 cm s' 1 at t = 20.0 min 
t = 20.0 minutes 


b) 
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v = 7.0 cm s 
t = 21.0 minutes 



Fig. 5.4 A forced ventilation velocity enhances the growth of a 
crystal particularly near its leading edges. 
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5.3 (d) and 5.3 (e). 

The linear growth of the crystal portrayed in Figs. 5.3 and 5.4 
is plotted in Fig. 5.5. Here, as in Fig. 5.2, the crystal length along 
the 'a'-axis as a function of time for a ventilation velocity of 7.0 cm 
s" 1 is well represented by a straight line. However, for no forced 
ventilation velocity, that is, after the velocity was decreased from 
7.0 cm s" 1 to 0 cm s”^ at time t = 0 minutes, the linear growth did not 
progress at a constant rate for about 6.0 minutes, i.e., until the 
doritic crystal tip had grown approximately 0.5 mm in length. The lower 
linear growth rate during this time period can be explained in terms of 
competition for the available water vapor between the primary doritic 
branch and the secondary branches. When the doritic branch had achieved 
a length of approximately 0.5 mm with no new branches this competition 
effect became much less significant and the doritic branch grew at a 
higher linear growth rate as indicated in Fig. 5.5. 

An example of the supersaturation dependent transition dendri tem- 
plate is shown in Fig. 5.6. Within 5.0 minutes after the velocity was 
reduced from 7.0 cm s“^ to 0 cm s”^ the growing dendrite tip had developed 
a rather substantial plate, as illustrated in Fig. 5.6 (d). This crystal 
later resumed a dendritic habit when a velocity of 7.0 cm s"^ was reintro- 
duced. The linear growth response of this crystal to the removal of the 
ventilation velocity is given in Fig. 5.7. This graph shows that the 
linear growth rate did not immediately achieve a constant value upon re- 
moval of the forced ventilation velocity, but rather decreased with 
time, approaching a nearly constant rate after an elapsed time of 1.0 
minute. Since the linear growth rate was higher for a period of 1.0 
minute than the equilibrium rate under the imposed ambient growth 
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Fig. 5.5 Length of the crystal in Figs. 5.3 and 5.4 as a function of 
time for successive ventilation velocities of 0 and 7.0 
cm s~1 , under constant ambient conditions of temperature, 
-13.9°C, and vapor density excess, 0.385 g m~3, at the 
height of the crystal base, 0.78 cm above the bottom plate. 
Angle of crystal from horizontal was +3°. (T+ = -2.3°C: 

T b = -19. 1°C) . 







Crystal Length (mm) 



0 1.0 2.0 3.0 4.0 5.0 6.0 


Time (minutes) 


Fig. 5.7 The linear growth response of an ice crystal under ambient 
conditions of temperature, -13.7°C, and vapor density 
excess, 0.352 g m _ 3, when the ventilation velocity, v, was 
suddenly decreased from v = 7.0 to v = 0 cm s~l at time 
t = 0. (Tt = -4.5°C; Th = -20.2°C). 
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conditions, it can not be argued that there was a vapor competition 
effect. Furthermore, since the time required for the vapor field 
surrounding the crystal tip to come to equilibrium with the new growth 
conditions was approximately 10 sec, an initial supply of water vapor 
near the crystal tip, in excess of the vapor demanded to sustain the 
new equilibrium linear growth rate of the crystal, can not explain this 
observation. Subsequent experiments revealed an explanation. For 
example. Fig. 5.8 illustrates the reverse transition plate-Klorite 
brought about by increasing the forced ventilation velocity from 0 to 
5.0 cm s -1 at time t = 0 minutes. Figures 5.8 (b) and 5.8 (d) show 
this transition in detail and Fig. 5.8 (c) indicates the orientation of 
the crystal both to the glass fiber and to the ventilation velocity. 

This crystal grew at an angle of +30° to the horizontal, somewhat larger 
than normally accepted. However, if the linear growth rates are used 
only over a small increment of length, for example 0.5 mm which is 
greater than necessary in the present case, then this crystal traverses 
a temperature gradient during this growth of less than 0.2°C and a 
moisture gradient of less than 0.3% supersaturation with respect to ice. 
In this case, when the ventilation velocity was increased from 0 to 
5.0 cm s _1 at time t = 0 minutes, it took over 2.0 minutes, as can be 
seen in Fig. 5.9, for the crystal to resume a constant linear growth 
rate. 

The constant linear growth rate along the 'a '-axis for no forced 
ventilation velocity was 0.4 urn s - ^ and for 5.0 cm s"^ was 1.5 pm s“^ . 

In Fig. 5.9, the linear growth rate increased with time in response to 
the introduction of a ventilation velocity at time t = 0 minutes in 
much the same way that the linear growth rate decreased with time in 
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v = 05 cm s' 1 at t=0 min 
t = 0 minutes 


b) 


v = 5.0 cm s 
t = 2.0 minutes 




2.0 mm 


250 ym 


Fig. 5.8 The crystal shape and linear growth rate response to a 

change in the forced ventilation velocity is a function of 
the ambient temperature, the ambient supersaturation and 
the magnitude of the change in the forced ventilation 
velocity. 
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Fig. 5.9 The linear growth response of an ice crystal under ambient 
conditions of temperature, -14.0°C, and vapor density 
excess, 0.243 g m’3, when the ventilation velocity was 
suddenly increased from v = 0 to 5.0 cm s"^ at time 
t = 0. (T t = -7.8°C; T b = -21.5°C). 
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Fig. 5.7 when the ventilation velocity was removed. Further experjr 
ments, documented with both 16 mm time lapse movies and 35 mm time 
lapse photographs, bore out this same symmetrical relationship and 
showed that the crystal shape and linear growth rate response to a 
change in the forced ventilation velocity is a function of the ambient 
temperature, the ambient supersaturation, and the magnitude of the 
change in the forced ventilation velocity. Near -14°C the time constant 
for. the transition of shape and growth velocity change was a maximum 
for the transition plate-Kiendrite. For the most extreme transitions of 
medium thick plates 200 ym thick) to dendrites or of dendrites to 
medium thick plates this time constant was on the order of one or two 
minutes. For the less extreme transitions of dorite-*dendrite or 
dendrite-Klorite the time constant was very small; typically, on the 
order of seconds. 

One of the original objectives of this study was to determine, in 
as much detail as possibie, the effect of a forced ventilation velocity 
on the linear growth rate, as well as on the shape of an ice crystal. 
Although the entire range of ambient temperatures from 0°C to -40°C 
was investigated, the most detailed studies were conducted over the 
temperature range -14°C ± 3°C. Over this range of temperatures, crystals 
were grown under ambient moisture conditions ranging from 4% to 50% 
supersaturation with respect to ice, i.e., at vapor density excesses 

q q 

ranging from 0.07 g m to 0.70 g m , and velocities ranging from 0 
cm s" 1 to 10 cm s - ^ . Plots of crystal length versus time, such as 
Fig. 5.2, were made for each crystal. At least five crystal length 
measurements at each velocity were employed to obtain the linear crystal 
growth rate under that particular set of ambient growth conditions. 



126 


The time period of crystal growth at each velocity was maintained long 
enough to insure that the crystal achieved linear growth rate equilibrium. 
When this condition was purposely not met, the linear crystal growth rate 
at any given ventilation velocity was determined, in part, by the pre- 
ceding velocity. 

The compilation of the linear crystal growth rates over this 
range of ambient temperatures and ambient supersaturations was achieved 
for velocities of 0, 3.0, and 7.0 cm s . Figures 5.10, 5.11, and 5.12, 
which were drawn in an attempt to make use of this compilation of 
individual crystal growth measurements, are plots of the linear 'a' -axis 
growth rates versus ambient temperature along the abscissa and ambient 
vapor density excess along the ordinate for velocities of 0, 3.0, and 
7.0 cm s~^ respectively. Utilizing these individual crystal 
growth measurements, isopleths of linear 'a '-axis growth rates in pm 
s" 1 have been drawn where possible. Error bars for both the ambient 
temperature and the ambient supersaturation include the effect of crystal 
growth occurring along vertical temperature and moisture gradients due 
to the fact that the crystals sometimes grew at an angle other than 0° 
from the horizontal. The fact that the error bars are longer at higher 
ambient vapor density excesses is indicative of the larger vertical 
gradients of both temperature and supersaturation necessary to produce 
these ambient moisture conditions. 

The most obvious effect which can be deduced from these three 
figures of the forced ventilation velocity on ice crystal linear growth 
rates is that the introduction of a velocity reduces the ambient super- 
saturation necessary to sustain the same linear growth rate. This 
result is predicted by theory. In fact, as shown in Chapter I II, theory 
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Fig. 5.10 Isopleths of linear 'a'-axis growth rates (vim s ) for a 
ventilation velocity = 0 cm s~'. Each circle denotes a 
set of individual crystal growth measurements. 
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Ambient Temperature ( °C) 


Fig. 5.12 
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Isopleths of linear 'a' -axis growth rates (ym s - ^) for a 
ventilation velocity = 7.0 cm s~l . Each circle denotes 
a set of individual crystal growth measurements. 
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predicts that the linear crystal growth rate for a given set of ambient 
temperature and supersaturation conditions is directly proportional to 
the square root of the forced ventilation velocity. On the other hand, 
theory does not predict that at higher ventilation velocities the 
maximum linear growth rate at a given ambient supersaturation occurs 
at successively colder temperatures. Although this result is not con- 
clusive from these three figures, due to the relatively large error 
bars, it is interesting to compare this observation with the result of 
other experimenters. Although he did not address this point, the 
experimental results of Fukuta (1969), Fig. 3.6, indicate that the 
total mass of an ice crystal falling through a supercooled water cloud 
has a maximum near -15°C at successively colder temperatures in direct 
proportion to the time after seeding. Of course, the longer the time 
after seeding, the larger the crystal and hence, the higher its termi- 
nal fall speed. If this observation is correct, it would explain why 
Fukuta, using ventilated freely falling crystals, observed the dual 
maxima in the mass growth rate versus temperature plot. Fig. 3.6, to 
occur at colder temperatures than found by Hallett (1965), Fig. 3.5, 
who obtained mass growth rates of crystals in a static environment. 
Further credence is given to this observation by the experimental re- 
sults of Ryan, et. al . , (1976), Fig. 3.7, even though they also did 
not address this point. They grew freely falling crystals in a super- 
cooled water cloud and found that the maxima in the linear 'a '-axis 
growth rate, near -15°C, occurred at an ambient temperature 2°C colder 
when the crystals were collected at times greater than 130 seconds 
after seeding than when they were collected at times less than 130 
seconds after seeding. 



131 


If this result is real, as it appears to be, the explanation may 
be in terms of surface heating due to increased absolute mass growth 
rates when a ventilation velocity is introduced. Figure 5.13, which 
presents isopleths of forced ventilation velocity as a function of the 
ambient vapor density excess along the abscissa and linear 'a' -axis 
growth rate along the ordinate at an ambient temperature of -14.5°C, 
was derived from Figs. 5.10, 5.11 and 5.12. This figure shows that the 
introduction of a ventilation velocity of only 7.0 cm s”^ at this 
ambient temperature is "roughly equivalent" to increasing the ambient 
vapor density excess in the static case by a considerable amount. In 
fact, for dendrites growing in a velocity field of only 7.0 cm s"^ this 
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amount may be greater than 0.3 g m . Conduction-diffusion equilibrium 
as discussed in Chapter IV shows that this would result in a crystal 
surface temperature rise of over 1°C. Larger ventilation velocities 
would correspond to even larger surface temperature rises. Thus, if a 
temperature dependent surface kinetic effect is responsible for the 
growth rate maxima near -15°C in the static case, then since the intro- 
duction of a ventilation velocity may increase the absolute mass growth 
rate sufficiently to induce a rise in the surface temperature of the 
crystal , due to the increase in the release of latent heat, a colder 
ambient temperature is necessary for the surface kinetic effect to 
operate optimally at the same crystal surface temperature as in the 
static case. 

It is interesting to note that reduced pressure, through an in- 
crease in the vapor diffusivity, also encourages the maximum mass growth 
rate to occur at colder temperatures, i.e.. Fig. 3.4. Thus, freely 
falling ice crystals in high clouds, at low pressure, may have their 
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growth rate maxima at an ambient temperature somewhat colder than t15°C\ 

Figure 5.13 demonstrates quantitatively how the introduction of a 
forced ventilation velocity reduces the ambient supersaturation neces- 
sary to sustain the same linear growth rate. It illustrates the 
supersaturation regimes at this temperature in which plates, sector 
plates and dendrites occur. It shows how a small increase in local 
supersaturation, brought about either by an increase in the ambient 
supersaturation or by the introduction of a forced ventilation velocity, 
can have a dramatic effect on both the crystal shape and linear growth 
rate. For an ambient temperature of -14.5°C, this effect is most pro- 
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nounced at an ambient vapor density excess near 0.20 g m , i.e., 
close to water saturation, where the transition plate-*dendrite occurs. 

Previous observations by Hallett (1965) indicated that in a static 
diffusion chamber, where crystals grew in an environment with the car- 
rier gas velocity almost zero, the transition from ice crystal plates 
to dendrites took place somewhat above water saturation. However, 
observations of ice crystals in the atmosphere have indicated that this 
transition takes place at an ice supersaturation somewhat less than 
that equivalent to water saturation. Figures 5.10, 5.11, 5.12 and 5.13 
show quantitatively how the transition plate->dendrite is a function of 
temperature, ambient supersaturation and carrier gas velocity. Further- 
more, these figures demonstrate how it is possible to grow dendrites in 
the atmosphere at ice supersaturations less than that equivalent to 
water saturation. That is, the crystals fall at their terminal veloci- 
ties, thus seeing a relative carrier gas velocity equal to their 
terminal fall speed which may typically be 20 to 30 cm s - ^ , i.e.. Fig. 
3.8. The effect of this relative carrier gas velocity is roughly 
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equivalent to increasing the ambient supersaturation. 

Figure 5.14 shows the relationship of the linear 'a'-axis growth 
rate, at an ambient temperature of -14°C, to the square root of the 
forced ventilation velocity. Each line is indicative of growth at one 
ambient supersaturation, expressed here in terms of the vapor density 
excess. In this figure (derived from Figs. 5.10, 5.11 and 5.12) the 
dashed lines are the extrapolation of the experimentally verified solid 
lines to higher forced ventilation velocities. Notice that the lines 
are quite straight where either plate growth or dendritic growth occurs, 
i.e., at either low or high ambient supersaturations. However, at in- 
termediate supersaturations the crystal shape at this temperature 
changes significantly for small changes of either the supersaturation 
or the ventilation velocity. Thus, the linear 'a'-axis growth rate is 
directly proportional to the square root of the ventilation velocity, 
as theory predicts, if the crystal shape does not change significantly. 

As shown in Chapter III, the present macroscopic theory of ice 
crystal growth rates is based on an electrostatic analogy, with the 
vapor density, p, being analogous to the electrostatic potential, v. 

A basic flaw of this theory, on the molecular scale, is that it assumes 
the crystal surface is at a constant potential or sees a uniform vapor 
density excess. In actuality, the vapor density excess is generally 
different from one crystal face to another on the same crystal , and a 
vapor density gradient often exists across an individual growing 
crystal face. In fact, analogous to the higher electric charge densi- 
ties at the corners of a charged conductor, the crystal corners and 
edges generally have the highest local vapor density excess with lower 
local excess vapor densities near the center of a crystal face. In the 
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Fig. 5.14 Isopleths of ambient vapor density excess (g m ) for 
an ambient temperature = -14 °C. 
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case of two-dimensional nucleation and growth, steps on the crystal sur- 
face are nucleated at the crystal edge and propagate toward the center 
of that crystal face. As the steps propagate, the local vapor density 
field over the crystal surface must continually respond to the move- 
ment of these vapor sinks as illustrated in Fig. 5.15. In this figure 
the dashed lines denote isopleths of constant vapor density and the 
solid lines show the flow of the water molecules. Since the actual 
shape of the steps is not known, no attempt has been made in Fig. 5.15 
to show their minutely exact effect on the vapor density field. Steps 
here are propagating from left to right. The impinging water molecules 
may be adsorbed between adjacent steps on the crystal surface and 
migrate short distances across the surface to a nearby step. As the 
steps move farther from the crystal edge toward the center of the 
crystal face the local vapor density excess decreases and the steps 


Water Vapor 



Fig. 5.15 The vapor field and flow of water molecules over a growing 
ice surface with steps propagating from left to right. 
Adapted from Scott (1977) 


slow down 


As the foremost advancing steps slow down, their collection 
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areas for surface diffused adsorbed molecules overlap and adjacent steps 
merge, further decreasing their advancement velocity. Both the rate at 
which the incoming water molecules can be incorporated into the crystal 
lattice at a given site and the velocity of propagation of the steps 
is a function of temperature and step height. These properties along 
with the magnitude of the flux of water molecules to a given site and 
hence the local vapor density excess are the basic controlling mechan- 
isms for the linear crystal growth rates along both the 'a' and 'c' 
axes and hence of the crystal habit. By changing the relative magni- 
tudes of these controlling mechanisms one can change the crystal habit. 
When a carrier gas velocity is introduced, the absolute and in most, if 
not all, cases even the relative magnitude of the flux of water mole- 
cules over various areas of a crystal face is increased. This is the 
result of vapor rich air flowing over the crystal surface and particu- 
larly enhancing the flux of water molecules to the crystal tip and 
leading edges. With the introduction of a ventilation velocity the 
"true" vapor density excess over the crystal lies between the curves 
in Fig. 4.16 and increases as the velocity increases. This results 
from the ventilation velocity contracting the ambient surfaces of con- 
stant vapor density and temperature as in Fig. 3.2 to very near the 
leading edge of the ice crystal. Diffusion of both water vapor to the 
growing crystal and heat from its surface then occurs through a boun- 
dary layer of approximate thickness 

6 'v (irDt)* 5 , 

as defined in Chapter III. Here D represents either the diffusivity 
of water vapor through air or the thermal diffusivity whichever is 
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appropriate for the boundary layer under consideration. For an ice 
crystal having a radius of curvature, r, at its leading tip, the diffu- 
sion boundary layer thickness there is a function of the velocity accord- 
ing to the relation 


6 'v 


v v > 


5 


where D is the diffusivity of water vapor through air, and v is the 

2 1 

magnitude of the velocity. For r = 10 cm, v = 10 cm s , and D = 0.2 


<5 % 250 pm. 

Since the boundary layer thickness is inversely proportional to the 
square root of the ventilation velocity, it decreases with increasing 
velocity, thus allowing diffusion to occur more rapidly through the 
boundary layer. At high velocity this is equivalent to a nearly verti- 
cal mass transfer curve in Fig. 3.10. Following the introduction of a 
ventilation velocity the local vapor density field over the crystal 
surface is readjusted, almost immediately. However, increased local 
vapor density excesses at constant temperature imply increased nuclea- 
tion rates and hence a greater number of steps. It is the propagation 
and redistribution of these steps on the crystal surface to new equili- 
brium values which leads to the observed time constants for the transi- 
tion of crystal shape and growth velocity change when one or more of 
the ambient growth conditions, such as the forced ventilation velocity, 
is altered. 

For example, the supersaturation dependent secondary habit transi- 
tion plate-Kiendrite occurs near -14°C when a forced ventilation velocity 
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is introduced at a critical ambient supersaturation, because steps 
nucleated on the prism faces at the corners of the plate do not have 
time to propagate to the center of that face before a whole series of 
new steps are nucleated. Thus, if the steps nucleate at the corner of 
the prism face, preferential growth occurs at the corner of the plate 
leading to doritic or dendritic growth. Figure 5.16 represents a situa- 
tion in which the local vapor density excess and hence the nucleation 
rate of steps is higher at the corner BB' than at the other corners. 

For example, this could represent the case in which a ventilation 



Fig. 5.16 A model of layer by layer growth on the prism faces. 

Adapted from Yamashita (1976) 

velocity flows perpendicular to BB' and parallel to CD, i.e., the 
corner BB' grows into the wind. If the crystal continues to grow as 
a plate, even though the vapor density excess and hence the production 
rate of new layers is higher at the corner BB' than at the corners 
AA' or CC 1 , then the growth of the prism faces (ABB'A 1 ) and (BCC ' B * ) 
is controlled by the layers which come from BB'. That is, the steps 
or layers nucleated at the corner BB' would in this case propagate 
beyond the center of the adjacent prism face, for example the prism 
face (BCC'B 1 ), and would merge with the steps or layers propagating 
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from the corner CC'. If, on the other hand, the crystal does not con- 
tinue to grow as a plate but rather grows preferentially from one or 
mor» corners, at least two distinct possibilities exist. In the first 
case, nucleation of steps or layers at the corner BB * may occur at a 
slow enough rate such that the collection areas for surface diffused 
adsorbed molecules of these adjacent steps do not overlap until they 
have propagated some distance from the corner at which they were nucle- 
ated. In this case sector plates result. In the second case nucleation 
of steps or layers at the corner BB' may occur at a very high rate in 
which case the collection areas for surface diffused adsorbed molecules 
of adjacent steps overlap and the steps merge before they are able to 
propagate very far from the corner at which they nucleated. In this 
case dorites or dendrites result. It is important to realize that the 
vapor density excess and thus the nucleation rate of layers or steps at 
adjacent corners of the crystal may be quite different. Also it is 
important to understand that the vapor density excess controls the 
nucleation rate of the steps but the vapor density gradient controls 
the propagation of those steps. That is, if the vapor density gradient 
from a corner is large and negative and the nucleation rate is high, 
the steps repeatedly merge and hence travel a much shorter distance 
along the prism face than if the nucleation rate and the vapor density 
gradient were both much smaller. This was demonstrated by a crystal 
which grew as a dorite at one corner, where the vapor density excess 
was larger, but continued growing as a plate at the other corners. In 
a ventilation velocity field it is even possible for steps nucleated 
at a corner such as CC' in Fig. 5.16 to first move into a negative 
vapor density gradient just away from the corner CC' and slow down. 
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then move into a positive vapor density gradient and speed up as they 
propagate across the prism face (BCC'B 1 ) toward the corner BB ' . 

Depending on the ambient temperature and supersaturation, two- 
dimensional nucleation may occur on either the prism faces, the basal 
faces, or simultaneously on both. Figure 5.17 illustrates a dendrite 
viewed edge-on and shows that the base of the dendrite was thicker 



* 1 .0 mm ^ 

Fig. 5.17 A dendrite viewed edge-on shows the base is thicker than 
the tip. 

than the tip. When the forced ventilation velocity, or carrier gas 
velocity, was increased, the dendrite grew an even thinner crystal tip 
than before. On the other hand, when the velocity was removed, the 
linear crystal growth rate diminished and the crystal tip thickened to 
about the same thickness as away from the edge. These results were 
obtained despite the fact that the ambient supersaturation remained 
constant. This implies that any steps present on the basal face propa- 
gated from the crystal base toward the crystal tip. This must be true, 
because no matter how fast the tip might have grown relative to the 
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microscopic steps, if steps had nucleated on the basal face near the 
tip, the tip would have become thicker than the crystal base. When the 
ventilation velocity was turned off, the velocity of the thin steps 
propagating on the basal face from the crystal base to the crystal tip 
was faster than the linear growth rate of the thicker crystal tip. Hence, 
the thin steps reached the crystal tip sooner than before and produced 
a crystal of almost uniform thickness. The crystal base should always 
be a little thicker than the crystal tip at -14°C unless all the steps 
were eliminated from the basal face and a 'step producer', i.e., the 
glass fiber, no longer existed near the crystal base. This implies 
that near -14°C the prism faces are fairly good nucleation surfaces 
and the basal faces are very poor nucleation surfaces. In fact, near 
-14°C two-dimensional nucleation probably occurs at the edge of the 
prism faces at ambient supersaturations, substantially below water 
saturation, but probably due to competition with the adjacent prism 
face, never occurs at the leading edge of the basal face, regardless 
of the ambient supersaturation. 

Predominant Growt h Along the ' c ' -Axis 

Near -9°C and a little above water saturation nucleation occurred 
at the edges of both the prism and basal faces but was slightly favored 
on the basal face. This resulted in the crystal habit being hollow 
columns, having a c/a ratio slightly greater than 1.0, with macroscopic 
stepped surfaces on both the basal and prism faces. Figure 5.18 illus- 
trates macroscopic steps and their propagation on the basal face of 
such a crystal. These macroscopic steps were a few micrometers in 
height and propagated very slowly. In addition to nucleating and 
propagating steps or layers on the basal face and thus growing along 
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a) 
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v = 3.0 cm s 
t = 0 minutes 


b) 


v = 3.0 cm s" 
t = 17.0 minutes 





v = 3.0 cm s I 1 

t = 63 minutes 250 ym 

Fig. 5.18 Macroscopic steps propagating across the basal face of 
a hollow column. 
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the 'c'-axis, this crystal grew laterally as well, i.e., by nucleation 
and growth on the prism faces. Thus, in the absence of any conflicting 
evidence, the growth of ice crystals at -9°C, near water saturation, 
was interpreted to be by two-dimensional nucleation with subsequent 
propagation of microscopic steps which merged in the way discussed 
earlier to form the observed macroscopic steps. 

The effect of temperature on the relative nucleation rates of the 
prism and basal faces was best illustrated by growing crystals at a 
temperature near -14°C where nucleation is favored on the prism faces 
and then abruptly transferring the crystals to temperatures near -5°C 
where nucleation is favored on the basal faces. Figure 5.19 illustrates 
the results of two such experiments. In the first experiment. Fig. 

5.19 (a), the crystals grew as dendrites into a ventilation velocity 
of 7.0 cm s" 1 at an ambient temperature near -14°C. At time t = 0 
minutes the glass fiber, on which the crystals were growing, was raised 
in the chamber so that the crystals were transferred to a region with 
an ambient temperature near -5°C. With the ventilation velocity main- 
tained at 7.0 cm s“^ the crystals immediately began nucleating and 
growing on the basal faces. This resulted after a period of 22.0 
minutes in columns growing with their 'c'-axes perpendicular to the 
ventilation velocity as in Fig. 5.19 (a). Hence, even though the prism 
faces at the leading edge of the crystal "saw" a larger local super- 
saturation or vapor density excess than the basal faces, due to the 
presence of the forced ventilation velocity, their growth rate was very 
small in comparison to growth along the 'c'-axis. In the second experi- 
ment, a crystal grew as a plate near -10°C, but when transferred to a 
temperature near -4°C at time t = 0 minutes the crystal immediately 







146 


sprouted needles from its corners, the regions of highest local vapor 
density excess, and at time t = 24.0 minutes had assumed the shape 
shown in Fig. 5.19 (b). In this case the velocity was only 3.0 cm s”^, 
and the plate was oriented at an angle from the direction into the air- 
flow of about 60°. Since the needles grew perpendicular to the plates, 
they were oriented at an angle of about 30° from the direction into the 
airfl ow. 

Although the relative rate of nucleation and growth on the basal 
and prism faces of an ice crystal is primarily a function of tempera- 
ture, it is also a function of the local supersaturation or vapor den- 
sity excess. Therefore, when a velocity is either introduced or re- 
moved, altering the local supersaturation over a crystal surface in a 
way already discussed, the resultant nucleation and growth of the crystal 
immediately responds on the microscopic scale to this new local super- 
saturation. However, the crystal does not achieve a new equilibrium 
growth rate and shape on the macroscopic scale until after a measurable 
time period. The modification of the relative rate of nucleation and 
growth on the basal and prism faces of an ice crystal induced by a 
change in the local supersaturation over the crystal was probably best 
demonstrated at an ambient temperature near -6°C and an ambient super- 
saturation of 8.9% with respect to ice, i.e., a vapor density excess of 
0.26 g m . In the absence of a forced ventilation velocity the crystal 
grew as a column, at an angle of 15° from the horizontal, i.e.. Fig. 

5.20 (a), with a growth rate along the 'c'-axis of 0.2 pm s”"*. When a 
forced ventilation velocity of 7.0 cm s" 1 was introduced at time t = 0 
minutes the crystal assumed a new preferred growth direction, as illu- 
strated in Figs. 5.20 (b) and 5.20 (c), which was at an angle of 8° from 
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Fig. 5.20 Ice crystal columns introduced to a forced ventilation 
velocity may assume preferred linear growth along a 
direction other than the 'c'-axis. 
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the 'c'-axis. The linear growth rate along this new axis steadily in- 
creased, after the velocity was introduced, until it assumed an equili- 
brium rate of 1.5 ym s~\ This corresponded to a component of growth 
along the 'c'-axis equal to 1.49 ym s"^ and a growth component along 
the 'a '-axis of only 0.2 ym s“^. When the velocity was removed at 
time t = 45.0 minutes, the crystal resumed growth entirely along the 
'c'-axis but at the higher rate of 0.35 ym s"^. This higher growth 
rate was presumably due to the column tip. Fig. 5.20 (d), being much 
narrower than before the velocity was introduced. Fig. 5.20 (a). When 
crystals were grown at a higher ambient supersaturation but still near 
-6°C, they also assumed a preferred direction of linear growth other 
than along the 'c'-axis, even in the absence of a ventilation velocity. 
Figure 5.21 shows a crystal growing at an ambient temperature of -5.8°C 
and an ambient supersaturation of 13.3% with respect to ice, i.e., a vapor 

_3 

density excess of 0.41 g m . In this case, under conditions of no 
ventilation velocity, the preferred direction of growth was 4° to the 
'c'-axis at a linear growth rate of 1.5 ym s~^ . This corresponded to 
a component of growth along the 'a'-axis of only 0.10 ym s“^ and a 
component of growth along the 'c'-axis of almost 1.5 ym s”^ . At a 
velocity of 7.0 cm s”^ , the preferred direction of growth was 7° to 
the 'c'-axis at a linear growth rate of 2.5 ym s“^. In this case, the 
component of growth along the 'a'-axis was 0.30 ym s”^ and the growth 
component along the 'c'-axis was 2.48 ym s“^ . It can easily be seen 
from these examples that only a small relative change in the growth 
rates along the * c ' and 'a '-axes is necessary to induce growth in a 
new preferred direction. 


149 


Fiber inserted at t = 0 minutes with v = 7.0 cm s 
v = 0 cm s~l at t = 58.0 minutes 


a) 


b) 



v = 0 cm s 
t = 1 hr 33 minutes 
v = Cb-7.0 cm s _1 at 
t = 1 hr 33 min 



v = 7.0 cm s" 
t = 1 hr 43 minutes 


2.0 mm 


Fig. 5.21 At higher ambient supersaturations preferred linear growth 
may occur along a direction other than the 'c'-axis, even 
in the absence of a forced ventilation velocity. 

Hallett (1965) grew crystals in a static vapor diffusion chamber, 
at temperatures near -4°C and at an ambient supersaturation near 250% 
with, respect to ice, which showed preferred growth along a direction 
25° ± 3° to the 'c'-axis direction. These crystals grew hollow prisms 
behind their tips in shielded regions of lower supersaturation much as 
did the crystal in Fig. 5.21. Hallett confirmed the direction of the 
'c'-axis by subsequently lowering the crystal to -15°C and growing 
dendrites on the ends of the prisms. With the realization that the 
preferred growth direction did not lie in a simple crystallographic 
direction, he called these crystals "22° needles". The observations 



150 


in this study indicate that there is nothing inherently special about 
the angle of 22° at -4°C, but rather, growth along a preferred direc- 
tion over the temperature range -4°C to -6°C is a function of the local 
supersaturation. That is, a change in the local supersaturation can 
induce a change, not only in the absolute growth rates, but also in the 
relative growth rates along the 'a' and 'c'-axes. 

B . Crystal Growth in an Environment with Dropl ets 

The riming efficiency of a crystal is determined by the crystal 
size, habit, and orientation, the relative carrier gas velocity, and 
the droplet size distribution and concentration. Since the crystal 
habit is itself a function of the carrier gas velocity and ambient 
supersaturation, riming of crystals in the atmosphere is a very complex 
process. In an attempt to gain an insight into the riming process and 
to demonstrate the versatility of the dynamic thermal diffusion chamber, 
a series of riming experiments were performed. Droplets of 10 ym mean 
diameter were generated and introduced into the chamber as described in 
Chapter IV. The crystals to be rimed grew from a stationary solid glass 
fiber suspended vertically between the two horizontal chamber plates. 
Figure 5.22 is a sequence of photographs showing the extent of riming 
on columnar crystals, growing into the wind, for four different venti- 
lation velocities. In this case, the top plate temperature of the 
chamber was -0.8°C and the bottom plate temperature was -11.5°C, giving 
a theoretical maximum ambient supersaturation in the chamber slightly 
greater than water saturation. Of course, the presence of high con- 
centrations of water droplets ensured that the supersaturation was 
close to water saturation. On the other hand, at higher ventilation 
velocities, i.e., greater than 15 cm s" 1 , the incoming undersaturated 
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air no longer had sufficient time to approach equilibrium with either 
the linear temperature or vapor density profiles between the chamber 
plates, before it came into contact with the crystals. Hence, in the 
absence of droplets at higher velocities the crystals sublimated rather 
than grew. Figure 5.23 is a sequence of photographs showing the extent 
of riming on dendritic crystals for four different ventilation velocities. 
In this case, the chamber top plate temperature was -9.0°C and the 
bottom plate temperature was -22.5°C. This also resulted in a theoreti- 
cal maximum ambient supersaturation in the chamber slightly greater than 
water saturation. Since, unfortunately, the droplet size distribution 
and concentration was a complicated function of the ambient super- 
saturation, the ventilation velocity and the horizontal and vertical 
position in the dynamic thermal diffusion chamber, it was impossible 
to obtain meaningful quantitative results from these riming experiments. 
However, some qualitative observations were made. As the ventilation 
velocity was increased from 0 to 20 cm s"^ , in the presence of high 
droplet concentrations of mean diameter 10 ym, the crystal growth was 
progressively controlled more by droplet accretion and less by vapor 
diffusion. In fact, in the presence of a ventilation velocity of 15 
cm s" 1 droplet accretion accounted for over 90% of the crystal growth 
in both these experiments. This is a much lower velocity for the onset 
of riming than is typically observed in the atmosphere. However, the 
droplet concentrations in these experiments were very high, estimated 
to be typically several thousand per cubic centimeter, whereas, even 
in heavily polluted continental clouds, droplet concentrations with 
this mean diameter rarely exceed 1000 cm" . The rimed columnar crystals 


in Fig. 5.22 appear much more delicate than do the rimed dendritic 
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Fig. 5.23 Rimed ice crystal dendrites. 
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crystals in Fig. 5.23. Nonetheless, in both these experiments even the 
most extensively rimed crystals retained their original crystalline axes 
orientation. This was demonstrated by allowing vapor diffusional growth 
to take place on the rimed crystals as, for example, is shown in Fig. 
5.23 (d). 

For ventilation velocities less than 1.0 cm s' 1 , it was observed 
through the stereo microscope that water droplets of 10 pm diameter or 
smaller evaporated and completely disappeared as they approached a 
growing ice crystal. The distance from the crystal at which the drop- 
lets disappeared was a function of the initial droplet size, the venti- 
lation velocity and the crystal shape. For a given crystal shape this 
distance decreased with increasing droplet size or increasing velocity. 
At higher velocities the droplets flowed around the crystal unless 
they encountered a stagnation point, in which case the droplet would 
decelerate and impact upon the crystal surface. Following impaction 
the evaporating droplet would skitter along the surface until it com- 
pletely evaporated, adhered to the surface and froze, or moved off the 
leeward side of the crystal. 

C. Growth of Discoid Crystals from the Vapor 

While growing columnar crystals from the vapor at ambient tempera- 
tures over the range -5°C to -7°C, I observed that adjacent to the solid 
glass fiber, in regions of low local supersaturation resulting from 
competition between crystals, thin discoid crystals grew, as for exam- 
ple, on the leeward side of the fiber in Fig. 5.24. These were the 
first vapor grown discoid ice crystals ever to be observed. They had a 
thickness of only about 10 to 20 ym and grew at a linear rate of 0.03 
ym s' 1 . When they obtained a diameter of approximately 300 ym they 
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proceeded to nucleate and grow columns perpendicular to their tips as 
illustrated in Fig. 5.24 (c). Note that Fig. 5.24 (c) is an enlarge- 
ment of the discoids and the crystal indicated with the arrow in Fig. 
5.24 (a), (b), and (d). Thus, it was established that the discoids grew 
in the basal plane. This, in itself, was interesting since the c/a 
axial ratio for crystals grown in this temperature regime was formerly 
believed to never be smaller than 1.0, much less 0.03 as was the case 
for these discoids. When the liquid nitrogen cooled, 250 pm diameter 
solid glass fiber was repeatedly inserted into the chamber at an ambient 
supersaturation equal to water saturation, that is, at an ice super- 

o 

saturation of 6% or a vapor density excess of 0.19 g m , the small 
crystals which nucleated and grew on the fiber were always columns when 
there was a large distance between them and always discoids when their 
population density was high. The time needed for these crystals to 
nucleate and grow to an observable size was of the order of minutes, 
i.e., much longer than the time needed for the liquid nitrogen cooled 
fiber to come to thermal equilibrium with the flowing airstream. Hence, 
the fact that the fiber was initially very cold did not determine the 
observable growth of the crystals. The repeated observations of either 
columns or discoids, depending on their population density is consistent 
with the argument that discoids grow in the temperature regime -5°C to 
-7°C at low local supersaturations. Since the tip of a discoid crystal 
was always thinner than its base, nucleation did not occur on the basal 
face near the tip. When, however, the local supersaturation was in- 
creased above a critical threshold value, by the introduction of a 
ventilation velocity, by reduced vapor competition between crystals, 
or by increasing the ambient supersaturation, nucleation occurred on 
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the basal face and the crystal habit became columnar. At still higher 
local supersaturations, the relative nucleation and rate of growth on 
the basal and prism faces, respectively, became even more pronounced 
and the columnar crystals became sheath or needle-like as did the upper 
column in Fig. 5.24 (d). 

Following up these initial investigations, McKnight (to be pub- 
lished) has also grown both discoids* and plates at an ambient tempera- 
ture of -5°C and an ambient supersaturation of 2.5% with respect to 
ice, i.e., a vapor density excess near 0.1 g nf^. McKnight employed a 
sealed vertical ampoule about 5 cm in diameter near the base and 15 cm 
in length containing, as the vapor source, a sugar-water solution of 
known concentration and vapor pressure. The ampoule was immersed in a 
well -control led temperature bath of glycol. The ice crystals grew from 
a solid glass fiber which was suspended vertically inside the ampoule 
and terminated about 1.0 cm above the sugar-water solution. The dis- 
coids grown in the ampoule had comparable thicknesses and growth rates 
to those grown in the dynamic thermal diffusion chamber. However, in 
this case, the local supersaturation was much nearer the ambient super- 
saturation and the only crystals which became columnar were those near 
the tip of the fiber directly over the sugar-water solution, and then, 
only when the fiber tip was positioned within a few millimeters of the 
solution surface. Even in the ampoule, the discoids rarely exceeded 
500 ym in diameter before they grew scalloped edges and subsequently 
faceted prism faces. As the discoids grew in diameter they thickened 
from the base to the tip. When they became plates, typically, 0.5 mm 
in thickness and 2 to 3 mm in diameter their linear growth rates were 
reduced to about 0.005 ym s~\ 
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Although neither the mechanism by which the vapor grown discoids 
grow, nor the reason why they grow scalloped edges and then faceted 
faces when they reach a critical diameter; is as yet understood, these 
experimental observations should contribute significantly to our total 
understanding of ice crystal growth. 

D. Growth of Crystals at Interm ediate to Low Ambient S uper saturation 
Another area of experimental study which is very interesting and 
also provides much needed information about basic crystal growth pro- 
cesses, but plays only a minor role under atmospheric conditions, is 
the growth of ice crystals at intermediate to low ambient supersatura- 
tion. Current ideas of crystal growth mechanisms differentiate between 
growth by continuous surface or two-dimensional nucleation on a molecular- 
ly smooth surface, and growth which takes place by propagation of dis- 
crete steps on an otherwise molecularly smooth surface. The discrete 
steps may originate from surface accidents - particulate or molecular 
impurity, from regions elsewhere at higher supersaturation where surface 
nucleation has already taken place, such as from corners as illustrated 
in Fig. 5.16, or from dislocations. If the crystal dislocations are 
screw- type, the steps perpetuate from growth layer to growth layer. 
Continuous two-dimensional nucleation over the crystal surface is be- 
lieved to be limited to higher ambient supersaturations, being under 
these conditions the dominant growth mechanism. Under very low ambient 
supersaturation growth conditions, dislocation initiated steps are 
believed to be the dominant growth mechanism. At intermediate to low 
ambient supersaturations, there has been considerable question as to 
what growth mechanism prevails. 

McKnight (to be published) using vapor grown ice crystals from my 


159 


dynamic thermal diffusion chamber, vapor grown ice crystals from an iso- 
thermal chamber containing a little supercooled water or supersaturated 
sugar solution as the vapor source, natural frost crystals, and ice 
crystals grown from a slightly super-cooled melt, has studied ice 
crystal dislocation densities under various ambient conditions. One of 
his observations, Anderson, Keller, McKnight and Hallett (1976) is that 
crystals sometimes grow, apparently dislocation free, at ambient super- 
saturations significantly lower than present theories would allow. 

At an ambient temperature near -30°C and at ambient supersatura- 
tions near 20% with respect to ice, i.e., at a vapor density excess 
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near 0.07 g m , simultaneous growth of ice crystal plates and columns 
under identical ambient conditions occurred in my dynamic thermal dif- 
fusion chamber as illustrated both in Fig. 5.25 and in the paper by 
Anderson, et al . (1976). At much lower ambient supersaturations the 
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Fig. 5.25 Simultaneous growth of columns and plates sometimes occurs 
at low supersaturation. 

crystals grew more or less equiaxed and at much higher ambient super- 
saturations the crystals grew as hollow columns. Table IV delineates 
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the habit of the ice crystals grown near -30°C as a function of the 
"effective" ice supersaturation. 

TABLE IV Habit of ice crystals growing at -30°C. 

"Effective" ice 


supersaturation 


Habit 

c/a 

0 




5 


equiaxed 

1.0 

10 

f 

simultaneous 

1.0 

15 


sol id 

to 

20 


columns and 

0.3 

25 


plates 




hollow 

5 

30 water saturation 



to 

-Kzero velocity) 



10 

>35 


col umns 

>20 


Anderson (1974) and Anderson, et. al . (1976) observed that ice 
crystals growing epitaxially on either covellite or silver iodide some- 
times have basal faces which do not grow at all in the presence of a 
supersaturation between water and ice but can be induced to grow by 
increasing the ambient supersaturation to a critical value which is 
variable from crystal to crystal due, at least in part, to strain 
originating between the crystal and substrate. 

The three independent observations of McKnight, Keller, and 
Anderson are consistent with the following hypotheses for growth within 
the range of conditions covered by these experiments: 

(a) growth at and somewhat below water saturation is initiated in 
the absence of defects by a two-dimensional surface nucleation process. 

(b) some crystal surfaces fail to grow at all at small ice super- 
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saturation; finite growth requires the presence of defects. 

(c) growth habit under low supersaturation will be determined by 
the presence of defects; presence of defects in only the 'a' or ' c ' 
direction will lead to growth only on that face and could give rise to 
crystals of habit opposite to that usually observed for a crystal with 
defects in both directions. 

(d) habit at higher supersaturation and the details of the skele- 
tal growth will be controlled by two-dimensional nucleation processes 

at the crystal corners and edges where the crystal protrudes into higher 
supersaturation regions. 

For crystal growth habit under atmospheric conditions, the control- 
ling factor will almost always be the surface nucleation process since 
most cloud processes take place at or a little below water saturation. 

In those situations in the atmosphere for which small ice supersatura- 
tions do exist, probably confined to regions of slow cooling by radia- 
tion or slow ascent situations with an abundant concentration of ice 
crystals, the habit, which is likely controlled by the initial defect 
distribution being retained during growth, may be opposite to that 
usually observed. 

The surface nucleation process at the crystal corners and edges 
and subsequent propagation of steps to regions of lower local super- 
saturation, as described earlier, was the controlling factor for the 
majority of the crystals grown in the dynamic thermal diffusion chamber. 
In fact, all the crystals taken from my chamber and subsequently x-ray 
topographed were apparently dislocation free. However, topographs of 
crystals grown in my chamber at low supersaturation remain to be taken. 
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E. Growth of 'Pecu liar 1 or Prism-Faced Crystals 

At temperatures near -30°C, some of the crystals which grew from 
the solid glass fiber in the dynamic thermal diffusion chamber had very 
peculiar shapes. Several of these crystals were similar to the 'peculiar' 
natural crystals reported by Nakaya (1954, pg. 65), Kikuchi (1970), and 
Kikuchi (1971) and also to the artificially grown crystals reported 
by Yamashita (1971). Kikuchi observed this type of crystal in both 
Antarctica and Japan and noted that it was frequently mixed with a com- 
bination of bullets, columns and crossed plates, comprising, in general, 
less than 1% of the total number of crystals falling from clouds at 
temperatures colder than -25°C but warmer than -35°C. The actual con- 
ditions of growth, however, were difficult to ascertain since the 
crystals grew in natural clouds but were normally collected at the 
ground. Yamashita, on the other hand, observed this type of crystal 
when he seeded a supercooled water cloud near -25°C with a metal rod 
previously immersed in either dry ice or liquid nitrogen. Yamashita 
believed the temperature of the seed rod might be an important factor 
in the subsequent formation of these peculiar crystals. Kobayashi (1976) 
suggested that some of these peculiar-shaped crystals may be explained 
as twins and he proposed a possible mechanism for their formation em- 
ploying a generalized Coincidence-Site Lattice (CSL) concept. Kobayashi 
suggested that the twinned seed may originate in the atmosphere from 
the freezing of water droplets to form polycrystals and, in particular, 
to form twins at temperatures lower than -20°C. He emphasized that 
although these crystals seem to have a peculiar shape, far from the 
normal hexagonal symmetry as a whole, they are indeed hexagonal and do 
not form a new classification of crystals. 
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Figure 5.26 shows two stages in the growth of a prism- faced or a 
particular 'peculiar' ice crystal which grew in the dynamic thermal dif- 
fusion chamber. This crystal grew at an ambient temperature of -28.6°C 
and an ambient supersaturation of 15.5% with respect to ice, i.e., a 

O 

vapor density excess of 0.06 g m . At all stages of growth of this 
crystal there was no forced ventilation velocity. However, in the 
latter stages of growth, illustrated by Fig. 5.26 (b), the local super- 
saturation was lower due to competition with other crystals and this 
crystal grew a column on the tip of the rib, parallel to the column 
already growing at the opposite end of the crystal. The direction of 
the 'c'-axis is indicated in the figure. The surface structure of this 
crystal was investigated using several different lighting arrangements, 
i.e., background lighting, direct illumination, and specular reflection, 
and was found to be very complex. This crystal, particularly in the 
stage of growth shown in Fig. 5.26 (b), was stepped from its center to 
its outer edges. The upper edge of the crystal containing the exposed 
basal faces of the two columns was much thicker than the lower edge. 

The basal face of the small column was by no means smooth. In fact, 
with the exception of the lower left-hand corner, this crystal had few 
smooth surfaces of significant area. Notice in Fig. 5.26 (a) that the 
rib grew at an angle of about 50° from the 'c'-axis and extended into 
the vapor field slightly beyond the rest of the crystal thus providing 
a possible nucleation site for continued growth. 

The effect of a forced ventilation velocity on the growth of this 
type of crystal is demonstrated in Fig. 5.27. This crystal grew at an 
ambient temperature of -33.7°C and an ambient supersaturation of 57.7% 
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with respect to ice, i.e., a vapor density excess of 0.135 g m . When 
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a ventilation velocity was introduced, it was only a few minutes before 
it was obvious that nucleation and growth was predominant at the tip of 
the rib as illustrated in Fig. 5.27 (b) and (c). Since the prism plane 
of this crystal did not lie quite perpendicular to the viewing field, 
it is difficult to accurately determine the angle between the rib and 
the 'c'-axis. However, this angle would certainly appear much smaller 
than the angle of 50° measured for the crystal in Fig. 5.26. 

The variety of unusual shapes which the 'peculiar' ice crystals may 
assume is demonstrated in Fig. 5.28. The number of 'peculiar' crystals 
grown was insufficient to accurately assess how the shape of these 
crystals is a function of both the ambient temperature and supersatura- 
tion, or if, in fact, an initial "twinned seed" nucleation process con- 
tributes significantly to later stages of growth of the crystal. How- 
ever, enough crystals were grown to provide the following qualitative 
observations. The crystals responded to an increase in the local super- 
saturation, i.e., through the introduction of a ventilation velocity, 
in much the same way as crystals at other temperatures. Nucleation and 
growth of 'peculiar' crystals occurred on the solid glass fiber even 
though the fiber was not cooled to either liquid nitrogen or dry ice 
temperatures, but, rather, was inserted into the chamber at room tem- 
perature. Nucleation and growth of 'peculiar' crystals occurred on the 
solid glass fiber even when the ambient supersaturation was maintained 
below water saturation. Over this range of ambient temperature and 
supersaturation the number of 'peculiar' crystals comprised a large 
percentage of the total number of crystals which grew. The angle be- 
tween the rib and the 'c'-axis varied significantly from crystal to 
crystal . 
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VI. 

APPLICATION OF RESULTS AND SUGGESTIONS 
FOR FUTURE STUDY 

A. Application to the Atmosphere 

Caution should be exercised in applying either the quantitative 
linear growth rates or the onset of riming criteria obtained from these 
experiments directly to the atmosphere. Since ice crystals in these 
experiments grew from a stationary support their orientation in the 
moving airstream was fixed. On the other hand, the flow pattern of 
air around the freely floating ice crystals in the atmosphere may be 
quite variable. In the atmosphere, ice crystal plates, for example, 
often fall in an oscillatory manner with their basal plane perpendicular 
to the air flow. In these experiments, the basal plane of plates was 
usually oriented parallel to the air flow, resulting in preferential 
growth of the crystal corner oriented into the airstream. The ventila- 
tion velocities employed in this study in the absence of water droplets 
were less than or equal to 10 cm s"^. This corresponds to terminal 
fall speeds of only very small crystals, typically less than 0.5 mm dia- 
meter. The presence of water droplets in the atmosphere both compli- 
cates the computation of the local supersaturation surrounding a grow- 
ing ice crystal and introduces the possibility of ice crystal growth 
occurring by a combination of vapor growth and riming, both these growth 
characteristics being functions of the droplet size distribution and 
the droplet concentration. Considering these differences in growth 
conditions, the rates of growth determined in these experiments when 
employing a ventilation velocity of 7.0 cm s - ^ , i.e.. Fig. 5.12, com- 
pare favorably with the linear growth rates determined by Ryan, et. al . 
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(1976), i.e.. Fig. 3.7, for small crystals 100 ym to 300 ym in diameter 
falling at velocities between 5 cm s - ^ and 10 cm s - ^ through an environ- 
ment containing a cloud of water droplets at or near water saturation. 

It should be noted in comparing linear crystal growth rates obtain- 
ed in the present experiments with those of Kobayashi (1958), i.e.. 

Fig. 2.5, Kobayashi (1961), Rottner (1971), or Rottner and Vali (1974) 
that they used a different definition of vapor density excess than 
normally employed in this study. They implicitly defined the vapor 
density excess over a growing crystal to be 


AP ■ PJU-P 1ce (T c ), (6.1) 

where p^T^) is the ambient vapor density and Pj ce (T c ) is the saturated 
vapor density over a plane ice surface at the temperature of that 
growing crystal surface. They calculated the ice crystal surface tem- 
perature, T , from a conduction-diffusion equilibrium for a static 
environment in the way described in Chapter IV, i.e.. Fig. 4.15. That 
is, in their experiments T was always given by T ' , i.e., the ice 
crystal surface temperature in a static environment. Their definition 
of the vapor density excess was defined as Ap stat ^ c in Fig. 4.16. In 
actuality, eq. (6.1) gives the definition of the "true" vapor density 
excess for either a static or ventilated crystal and it could be used 
successfully if the crystal surface temperature, T c , were known. 

However, since the surface temperature of a growing ice crystal is a 
function of the crystal shape, the ambient supersaturation and the 
ventilation velocity its accurate determination is immensely difficult . 
Hence, in this study the environmental vapor density excess was employed. 
This quantity is defined to be 
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Ap = Poo( T oo)-Pice^ T oo)’ ( 6 - 2 ) 

where p.-„(T_) is the saturated vapor density over a plane ice surface 
at the ambient temperature . That is, this definition gives the vapor 
density excess over a non-growing crystal. This definition of the 
vapor density excess was defined as Ap max in Fig. 4.16 and is more nearly 
equivalent to the normal definition of ambient supersaturation as used 
in the 'Nakaya diagrams' of Hallett and Mason (1958a) and Mason (1971), 
i.e., Fig. 2.4. 

Exercising the aforementioned cautions, the significance of these 
experimental results as applied to the atmosphere may be briefly sum- 
marized. I have demonstrated how it is possible to grow dendrites and 
needles in the atmosphere at ice supersaturations less than that equi- 
valent to water saturation. That is, as the crystals fall they "see" 
a relative carrier gas velocity equal to their terminal fall speed. 

The effect of this carrier gas velocity is roughly equivalent to in- 
creasing the ambient supersaturation. This is an important result 
since accretion of droplets is a function of the crystal habit, with 
skeletal or dendritic structures offering larger surface areas and 
having better riming characteristics. A very important process in the 
atmosphere may be the initial growth of a small ice crystal plate, for 
example, near water saturation and at low terminal fall speed. As the 
crystal grows and reaches a higher terminal velocity its linear growth 
rate may increase substantially. As a result of the increased relative 
carrier gas velocity the crystal may go through the transition to a 
dendrite which offers a much larger collection area for both water 
droplets and pollutants. Hence, both the initiation of precipitation 
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and scavenging of pollutants may occur faster in these cases than has 
previously been predicted. Since the crystal growth rate is a function 
of temperature, having maxima near -6°C and -15°C which are a function 
of the crystal fall speed, optimum injection of ice nuclei in cloud 
seeding experiments utilizing dynamic cloud growth is a function of 
the temperature at the level of injection and the stage of cloud 
maturity. 

Finally, I have shown that when an ice crystal undergoes a plate 
to dendrite or vice-versa transition there is a time constant for a new 
equilibrium linear crystal growth rate to be established. This means 
that the history of a crystal having fallen through changing ambient 
conditions of temperature, supersaturation and ventilation velocity 
cannot be uniquely specified. 

B . Ap plication to Crystal Growth 

Diffusion of both mass and heat takes place through the boundary 
layer surrounding a growing crystal . To facilitate good transfer, this 
layer should be as thin as possible. Since the boundary layer thickness 
is inversely proportional to the square root of the ventilation velocity, 
as shown in Chapter V, an increase in the velocity decreases the boun- 
dary layer thickness. (In practice the boundary layer thickness is 
slightly different for the two processes). A decrease in the boundary 
layer thickness increases the flux of water vapor to and heat from a 
growing ice crystal and also induces partial pressure gradients of water 
vapor from the leading edge across the crystal surface. Hence, enhanced 
growth in a preferred direction results, and needles or dendrites, for 
example, grow depending on the ambient temperature. 

The introduction of a ventilation velocity is not equivalent to 
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increasing the diffusion coefficient of water vapor through air, as a 
comparison of the results from these experiments with those shown in 
Fig. 2.10 demonstrates. That is. Fig. 2.10 shows that at a constant 
thermal conductivity of the air an increase in the diffusion coefficient 
leads to increased growth rates but to a decrease in the development of 
skeletal or dendritic structures. The experiments described in this 
study show that the introduction of a ventilation velocity leads to in- 
creased growth rates and to an increase in the development of skeletal 
or dendritic structures. This is consistent with the fact that an in- 
crease in the diffusion coefficient not only increases the flux of 
water molecules to a growing ice crystal but also increases the boun- 
dary layer thickness in direct proportion to the square root of the 
diffusion coefficient as shown in Chapter V. This results in enhanced 
growth, but since sustained partial pressure gradients of water vapor 
across the crystal surface are in this case smaller, the growth is also 
more uniform. This is also consistent with the observations of Lamb 
(1970), as described in Chapter II, that the surface structure of ice 
crystals grown in the presence of either air or helium at low pressure 
was more complicated than that of crystals grown in pure water vapor. 

Figure 2.10 shows that under conditions of constant vapor dif- 
fusivity skeletal and dendritic structures developed more readily as the 
thermal conductivity of the atmosphere was increased. This is not un- 
reasonable, since an increase in the thermal conductivity of the air 
or any other process which lowers the growing crystal surface tempera- 
ture corresponds to increasing the "true" vapor density excess, i.e.. 


supersaturation, over the crystal surface as described in Chapter IV. 
Thus, an increase in the thermal conductivity is nearly equivalent to 
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increasing the ambient supersaturation. 

It is important to appreciate the fact that a ventilation velocity 
is not equivalent to increasing the diffusion coefficient. On the other 
hand, the ventilation velocity does not enhance growth by just lowering 
the crystal surface temperature. Rather, in addition to lowering the 
crystal surface temperature by conducting heat away, the ventilation 
velocity decreases the thickness of the boundary layer preferentially 
over the crystal surface. This increases the flux of water vapor to 
the crystal while maintaining partial pressure gradients of water vapor 
across the crystal surface. 

In many crystal growth applications, such as in industry, it is 
desirable to grow uniform crystals in the least possible time. For 
these applications, the results obtained for ice crystals indicate that 
at a fixed supersaturation a large thermal conductivity of the carrier 
fluid is desired to promote higher growth rates, while maintaining a 
large ratio of the diffusion coefficient to the thermal conductivity 
to induce uniform growth. It is also important in these applications 
to maintain as uniform fields of temperature and vapor surrounding the 
crystal as possible. 

In practice it is very difficult to maintain uniform fields of 
either vapor or temperature around a growing crystal. In addition to 
'inherent' inhomogeneities in the vapor and temperature fields, the 
release of latent heat at the growing crystal surface often leads to 
natural convection. In the experiments described in this study the 
natural convection resulting from a growing dendrite was estimated to 
be 1 mm s"^ . This measurement was obtained by observing, with a stereo 


microscope equipped with a micrometer eyepiece, entrained water droplets. 
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approximately 5 pm in diameter, convecting upward from growing dendrites 
whose basal plane was parallel to the horizontal. Smaller values of 
natural convection result from more slowly growing crystals. 

Natural convection can be greatly reduced in 'Zero Gravity' experi- 
ments conducted in space, resulting in more uniform fields of both 
vapor and temperature. Consequently, crystal growth under these con- 
ditions should be more uniform. 

C. Suggestions for Future St udy 

Several ideas for possible future studies exist which employ the 
dynamic thermal diffusion chamber. If the working section of the cham- 
ber were doubled or tripled in length to obtain equilibrium ventilation 
velocities of 20 to 30 cm s - ^ , while retaining the nearly linear ver- 
tical profiles of both temperature and vapor density, the present ex- 
perimental results could be extended to crystal fall speeds approaching 
those encountered in the atmosphere. If the working section of the 
chamber were modified, consideration should be given to increasing the 
chamber aspect ratio. A greater number of quantitative growth rates 
of columns and needles as a function of temperature, supersaturation and 
ventilation velocity are needed to enlarge upon the present results. 
Quantitative growth rates should be obtained for crystals at various 
orientations to the air stream to better assess the direct applicability 
of the present results to growing ice crystals falling freely in the 
atmosphere. Both the study of discoid crystal growth at low super- 
saturation and the study of prism- faced or 'peculiar' crystals at low 
temperature should be continued in conjunction with x-ray topographic 
studies of these crystals. An intensive study involving the growth of 
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scalloped plates near 0°C would be interesting. A comparison might be 
made between their growth and the growth of scalloped plates from dis- 
coids near -6°C. It would be interesting and beneficial to study the 
simultaneous effect, on the growth rate and shape of ice crystals, of 
a direct current electric field and a ventilation velocity, particularly 
if they were applied perpendicularly to one another. Perhaps such a 
study would answer conclusively the question of whether or not the 
application of a d.c. electric field is equivalent to the introduction 
of a ventilation velocity. 
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APPENDIX A 

AIR DRIVE MEASUREMENTS 

Figure A-l is a plot of the motor setting in microamperes, i.e., 
a measure of the regulated d.c. current to the servo motor which drives 
the shaded pole blower fan, versus the shaft speed of the motor, i.e., 
motor speed in RPM. The shaft speed of the motor was determined with 
a General Radio type 1538-A Strobotac electronic stroboscope. At low 
RPM the stroboscope calibration was not very accurate so there was some 
scatter in the experimental points at low motor settings. 

Figure A-2 is a diagram showing the ventilation velocity in the 
working section of the chamber as a function of the motor setting in 
microamperes for one particular freezing of the water in the moisture 
conditioner, i.e., one static pressure. The bars on the experimental 
points in Fig. A-2 are not error bars. Rather, they denote the 
variability of the measured ventilation velocity, i.e., the magnitude 
of the turbulence, for that particular motor setting. Figure A-2 shows 
that the transition from laminar to turbulent flow occurred around 40 
cm s' 1 with a turbulence level of ± 2.5 cm s"^. 
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APPENDIX B 
WATER VAPOR TABLES 

Tables B-I and B-II give, respectively, the saturated vapor pres- 
sure in millibars over plane ice and water surfaces as a function of 
temperature in degrees Celsius. Tables B-I I I and B-IV give, respectively, 

_3 

the saturated vapor density in g m over plane ice and water surfaces 
as a function of temperature in degrees Celsius. The values in these 
four tables were obtained by evaluating the Goff-Gratch equations at 

0. 1 °C temperature increments, taking into account through a compres- 
sibility factor the deviation of water vapor from the ideal gas law. 

The value in making computations at 0.1 °C temperature increments, rather 
than at 0.5°C increments and interpolating, lies mainly in internal 
consistency resulting in smoother calculated supersaturation curves, 

1 . e. , Figs. 4.12 and 4.13. 

No completely satisfactory formula exists for the vapor pressure 
over liquid water at temperatures below 0°C. However, direct extra- 
polation of the Goff-Gratch formulation gives values of the saturated 
vapor pressure with respect to water in the middle of the range suggest- 
ed by other experimental investigators. Therefore, this formulation 
was adopted in Resolution 164 of the Twelfth Conference of Directors of 
the International Meteorological Organization (Washington, 1947) for 
the range 0°C to -50°C pending further research. 

Table B-V gives the percent supersaturation with respect to a 
plane ice surface at which water saturation occurs as a function of 
temperature in degrees Celsius. 
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.1343 

.1327 

.1313 

.1298 

-40.0 

.1253 

.1269 

.1254 

.1240 

.1226 

.1212 

.1199 

.1105 

.1172 

.1159 

-41.0 

.1145 

. 1132 

.1120 

.1107 

.1094 

.1002 

.1069 

.1357 

.1045 

.10 3$ 

-42.0 

.10213 

•'13 0 9 6 

.09951 

.09066 

.09753 

. 09641 

.09530 

• u9 420 

.09312 

.09204 

-43.0 

.0909ft 

.08 99 3 

.08089 

.00786 

.08685 

.08584 

.00484 

.08306 

.00208 

.08192 

-4 4.0 

.05097 

. 08002 

•07909 

.07817 

.07725 

.07635 

.07546 

.07457 

.07370 

. 0 7 2 ft 4 

-45.0 

• 0719ft 

.07 113 

.07030 

.06947 

.06865 

.06784 

.06734 

. 06625 

.06547 

• G f » 4 6 9 

—46 .0 

. 06393 

. 0631 7 

.06242 

.06160 

.06094 

. 06022 

• 0 6950 

.05079 

. 0 5809 

.05740 

-47.0 

• 0 56 7 t 

. 05603 

.05536 

.05470 

.05405 

.05340 

.05276 

.05212 

.05149 

.05007 

-4 5.0 

.05026 

.'04 965 

*04905 

• 0 4846 

.04788 

.04730 

.04672 

.04616 

.04560 

.04504 

-49.0 

.0 4449 

. 04395 

•04342 

.04289 

.04237 

. 04185 

.04134 

• 04 383 

• 04033 

.03934 

-50.0 

.03935 

. 03886 

.03839 

.03791 

.03745 

. 03699 

.03653 

. 03 60 0 

.0 3563 

.03619 

-51.0 

.03476 

.03433 

•03390 

•0334ft 

.03306 

. 03265 

.03225 

. 031HS 

.03145 

.03106 

-52.0 

.03067 

.03029 

.02991 

.02953 

.02916 

.02880 

.02843 

.02800 

.0 2 772 

.02737 

—53.0 

• 027 3 3 

. C2669 

• 0 2635 

.02602 

.02569 

.02537 

.32604 

.02475 

.02441 

. 0 ?4 10 

-54.0 

.02350 

. 02 34 9 

•02319 

• uZ29G 

.02261 

.02232 

.02203 

.02175 

.02147 

. j212u 

-55.0 

.02392 

. 02066 

.02039 

. C2U13 

.01987 

.01961 

.01936 

.91911 

.0183b 

.01062 

— 5 6.0 

.0133ft 

.'01814 

. 91/90 

.01767 

.01744 

.01722 

.01699 

.01677 

.C 1655 

.0 16 54 

-57.0 

.01612 

. 01591 

.01670 

.01550 

.01629 

.01509 

.01489 

.01470 

.01451 

.41431 

-5 5.0 

.01413 

. 01 394 

.01376 

.01357 

•01339 

.01322 

.01304 

.012B7 

. u 1270 

. 0 1253 

-59.0 

.0 12 3 6 

.01220 

.01203 

.01187 

.01172 

.01156 

.01140 

.01125 

•0111U 

.01 095 



Table B-II 


SATURATION VAPOR PRFSSURE OVER MATER 


TEMPCC) 

.0 

• 1 

-2 

.3 

UNIT 

MB 

MB 

MB 

MB 

-0.0 

6. 1 >J7 8 

6.0636 

6.0196 

5.9759 

-1.0 

5.6780 

5.6366 

5.5954 

5.5544 

-2.0 

5.2753 

5.2364 

5.1979 

5.1595 

-3*0 

4.8931 

4.3617 

4.8256 

4.7397 

-<♦.0 

4.5451 

4.51 11 

4.4773 

4.4437 

-5*0 

4.2148 

4. 1830 

4.1514 

4.1200 

-6.0 

3.9061 

3. 8764 

3.3469 

3.8175 

-7.0 

3.6177 

3.5399 

3.5623 

3.5349 

-6.0 

3.3484 

3. 3224 

3.2967 

3.2711 

-q.o 

3.0971 

3.0729 

3.0439 

3.0250 

-10*0 

2.862 7 

2.3402 

2.3173 

2.7955 

-11. 0 

2.6443 

2.6233 

2.6024 

2.5317 

-12.0 

2.4409 

2.4213 

2.4019 

2.3826 

— 13.0 

2.2515 

2.2334 

2.2153 

2.1974 

-14.0 

2.0755 

2.0585 

2.0417 

2. (251 

-15.0 

1.9118 

1.3961 

1.8805 

1.8650 

-16.0 

1.7597 

1.7451 

1.7307 

1.7163 

-17.0 

1.6136 

1.6051 

1.5916 

1.5783 

-18.0 

1.4877 

1.4752 

1.4627 

1.4503 

-19. 0 

1.3664 

1.3543 

1.343? 

1.3318 

-go*o 

1.254 0 

1.24 32 

1.2326 

1.2219 

-21.0 

1.1500 

1.1400 

1.1301 

1.1203 

-22.0 

1.0538 

1.0446 

1.0354 

1.0264 

-23.0 

.9649 

.9564 

.9479 

.9395 

-24.0 

. 382 7 

.8749 

.8671 

.8594 

-25.0 

. 8070 

.79 97 

.7925 

.7854 

-26.0 

.7371 

. 7304 

.7238 

.7172 

-27.0 

.6727 

.66 66 

.6605 

.6 544 

-28.0 

.6134 

• 69 7 8 

.6022 

.5966 

-29.0 

.5539 

. 55 37 

.5485 

.5434 

-30.0 

.5 03.3 

.5040 

.4993 

• 4946 

-31.0 

.4628 

.4584 

.4540 

.4497 

-32.0 

.4205 

.4165 

.4125 

.4086 

-33.0 

.3818 

.3781 

.3745 

.3701 

-34.0 

. 3463 

. 3430 

.3396 

.3363 

-35.0 

.3139 

.3108 

.3077 

.3047 

-36.0 

.284 2 

.2814 

,2736 

.2758 

-37.0 

.2571 

.2*45 

.2519 

.2494 

-38.0 

.2323 

.22 99 

.2276 

.2253 

-39.0 

.2097 

.20 76 

.2054 

.2033 

-40.0 

.1891 

.13 72 

.1153 

.1833 

-41.0 

.1704 

• 1636 

.1669 

.1651 

-42.0 

.1534 

.1518 

.1502 

• 1486 

-43.0 

.1379 

.1364 

.1350 

.1336 

-44.0 

.1239 

.12 25 

.1212 

.1199 

-45.0 

.1111 

. 10'99 

.1087 

• 1076 

-46.0 

. 0 9961 

.09 852 

.09744 

.09637 

-47.0 

.08918 

.03819 

* 08722 

.08625 

-48.0 

.07975 

. 07886 

.07793 

.07710 

-49.0 

.07124 

. 07 043 

.06964 

.06885 

-50.0 

.06356 

. 06233 

.06212 

•06141 

-51.0 

.05664 

.0559° 

.05534 

• 0 547 J 

-52.0 

. 05042 

. 04 933 

.04925 

.04367 

-53.0 

.04482 

.04429 

.04377 

•04326 

-54.0 

.03980 

. 03 932 

.03836 

.03839 

-55.0 

.03529 

. 03487 

.03445 

• C 3 4 0 3 

-56.0 

.0,5126 

. 03033 

.03050 

•03113 

-57.0 

.02765 

. 02731 

.02697 

•02664 

-58.0 

• 0 24 4? 

.02412 

.02332 

.02352 

-59.0 

.02154 

.02127 

•Q21C0 

.02074 


.4 

.5 

' • 6 

.7 

• 8 

• 9 

MB 

Mb 

MB 

Mt> 

M3 

Mu 

5.9325 

5. 8894 

5.8466 

6. 80 4C 

0. 7618 

9.7198 

5.5138 

5. 4734 

5.431? 

5.3934 

5.3537 

5.3144 

5.1214 

5.0836 

5.0460 

5.00 87 

4.9716 

4.9347 

4,7541 

4.71 H7 

4.6835 

4.64 85 

4 . b 1 3 .8 

4.5/93 

4.4103 

4.3772 

4.3443 

4.3116 

4.2791 

4.2469 

4.0388 

4.0579 

4. 3271 

3.9966 

3.9662 

i .9351 

3.7884 

3. 7594 

3.7307 

3. 7Q22 

3.6733 

3.6457 

3. 5077 

3.4807 

3.4539 

3.42 72 

3.4008 

3.3745 

3.2457 

3. 2205 

3.1955 

3.1706 

3.1 459 

3.1214 

3.0013 

2.9771 

2.9545 

2.9313 

2.9183 

2.8854 

2.7735 

2.7516 

2.7298 

2.7012 

2,6867 

2.6654 

2.5612 

2.5403 

2.5205 

2.5004 

2.4804 

2.4606 

2.3635 

2.3445 

2.3256 

2.3069 

2.2883 

2.2699 

2.1796 

2.1619 

2.1444 

2.1269 

2.1097 

2 • G 925 

2.0085 

1.9921 

1.9758 

1.9596 

1.9436 

1.9276 

1.8496 

1.8343 

1.8192 

1.8042 

1.7893 

1 .7744 

1.7020 

1.6879 

1*67 38 

1.6599 

1.6460 

1.6323 

1 .5651 

1.5519 

1.5389 

1.526U 

1.5131 

1.5004 

1.4381 

1.4259 

1.4138 

i.4019 

1.3899 

1.3761 

1.3204 

1. 3091 

1.2979 

1.2868 

1.2758 

1.2649 

1.2114 

1. 20 10 

1.1906 

1.1803 

1.1 701 

1.1600 

1.1106 

1. 1009 

1.0914 

1.0819 

1.0724 

1.0631 

1.0174 

1. QC 84 

.9996 

.9908 

.9821 

.9734 

.9312 

.9230 

.9141 

.9067 

.8 986 

.8907 

.8517 

.8441 

• 8365 

. 8291 

.3216 

.8143 

.7783 

.7713 

• 7 644 

.75 75 

.7506 

.7438 

.7107 

.70 42 

.6978 

.6915 

.6852 

.6789 

.6484 

.6425 

.6366 

• 6307 

.6249 

.6191 

.5911 

.5856 

.5802 

.5748 

.5695 

• 5642 

.5384 

. 5333 

.5283 

.5234 

.5185 

.5136 

♦ 4899 

.4853 

.4807 

.476 2 

.4717 

.4672 

.4454 

.4412 

.4370 

.4328 

.4287 

.4246 

.4046 

.4007 

.3969 

.3931 

.3 893 

.3855 

.3672 

.3637 

.3602 

.3567 

.3532 

.3497 

.3330 

.3297 

.3265 

.3233 

.3201 

.3170 

.3017 

.2987 

.2957 

.2926 

.2899 

.2870 

.2730 

.2703 

.2676 

.2649 

.2625 

.2597 

.2459 

.2444 

.2419 

.2395 

.2371 

.2347 

• 2230 

• 2207 

.2185 

.2163 

• 2141 

.2119 

.2013 

.1992 

.1971 

.1951 

.1931 

.1911 

.1814 

.1796 

.1777 

.1758 

.1740 

.1722 

.1634 

• 1617 

.1600 

.1583 

.1567 

• 1550 

.1470 

.1455 

.1439 

.1424 

.1409 

.1394 

.1321 

.1307 

.1293 

.1279 

.1266 

.1252 

.1186 

.1174 

.1161 

.1148 

• 1136 

.1124 

• 1064 

.1052 

.1041 

• 1030 

.1018 

. 1 Du7 

.09531 

.09426 

.09323 

.09220 

•09116 

.09018 

.08529 

.08435 

.08341 

.08248 

•06156 

•08065 

.07624 

.07531 

.07454 

.07370 

• U 7287 

.07205 

.06807 

.06730 

.06653 

.06573 

.06503 

.06429 

.06070 

•Q6GU1 

.05932 

.05864 

.05797 

• 0 573G 

.05407 

.05345 

.05213 

.05221 

.05161 

.05101 

• 048 11 

.04754 

.04699 

.04644 

• U 4589 

.04535 

.04275 

.04224 

.04174 

.04125 

.0 4076 

.04G27 

.03793 

.03748 

.03703 

.03659 

.03615 

.03572 

.03362 

.03322 

.03282 

.03242 

. 0 320 3 

• u 31 64 

.02976 

.02943 

.02904 

.02869 

•02834 

.02799 

.02631 

.02599 

.02567 

.32535 

•02504 

.02473 

.02323 

.02294 

.02265 

.02237 

.02209 

• C 2 1 12 

•02048 

.02022 

.01997 

.01972 

.01947 

.01922 


Table B-III 


00 

ro 


DFNStTY OF PURE WATER vA^OR AT SATURATION OVER ICF 


TEHPICI 

• 0 

.1 

• 2 

.3 

.4 

.5 

.6 

. 7 

• 8 

. 9 

UNIT 

G/M 3 

G/M 3 

G/M 3 

G/M 3 

G/M 3 

G/M 3 

G/M 3 

6/M 5 

G/M 3 

o/\\ 3 

-0.0 

4.847 

4.309 

4.771 

4.734 

4.696 

4.659 

4.623 

4.586 

4.550 

4.514 

-1.0 

4.4 79 

4.44 3 

4.408 

4.373 

4.339 

4. 3C4 

4.273 

4.23 6 

4.213 

4.16 9 

- ?• Q 

4.136 

4. 10 3 

4.070 

4.038 

4.006 

3.974 

3.942 

5.910 

3.8 79 

3.P48 

-3,0 

3.817 

3.737 

3.756 

3.726 

5.696 

3.666 

3.637 

3.608 

3.578 

3.55U 

-4.0 

3.521 

3.492 

3.464 

3,436 

3 .408 

3. 331 

3.353 

3. 526 

3.299 

5.2 72 

-5.0 

3.246 

3.219 

3.193 

3.167 

3.141 

3.115 

3.090 

3.065 

3 .039 

3 * v 15 

-6.0 

2.990 

2.965 

2.941 

2.917 

2.893 

2.869 

2.345 

2.822 

2.799 

2. 775 

-7.0 

2.763 

2.730 

2.707 

2.685 

2.662 

2.640 

2.613 

2 . ‘>9 7 

2.575 

2,954 

-8.0 

2.532 

2.5U 

2.490 

2.470 

2.449 

2.428 

2.408 

2.338 

2.368 

2.348 

-q.o 

2.323 

2.30 9 

2.289 

2.270 

2.251 

2.232 

2.213 

2.195 

2.1 76 

2.15*8 

-10.0 

2.119 

2.121 

2.103 

2. 085 

2.067 

2.050 

2.032 

2.U5 

1.998 

i. m 

-11.0 

1.964 

1.94 7 

1.931 

1.914 

1.898 

1.832 

1.865 

1.349 

1.8 34 

1.818 

- ! 2.0 

1.8 02 

1.737 

1.771 

1.756 

1.741 

1.726 

1.711 

1.496 

1.682 

1.667 

-13.0 

1.653 

1.633 

1.624 

1.610 

1.596 

1.58? 

1.563 

t.555 

1.541 

1.523 

-14.0 

1.514 

1.501 

1.438 

1.475 

1.462 

1.449 

1.437 

1.42 4 

1.411 

1 . 3 99 

-15.0 

1.3 87 

1. 374 

1.362 

1.350 

1.338 

1.327 

1.315 

1 • !G 3 

1.292 

1 , Z 8 U 

-16.0 

1.269 

1.258 

1.246 

1.235 

1.224 

1.214 

1.203 

1.192 

1.181 

1.171 

-17.0 

1.160 

1.15 0 

1.140 

1.129 

1.119 

1.109 

1.099 

1.039 

1.080 

1 . 0 7 J 

-18.0 

1.0603 

l . 05 0 7 

1.041? 

1.0318 

1.0225 

1.0132 

1.0041 

. 1950 

.9859 

.9770 

-15.0 

.9631 

.9593 

.9506 

.9419 

.9333 

. 9243 

.9164 

• 9 0 8 C 

.8997 

.8915 

-20.0 

.8833 

.37 52 

.8671 

. 8591 

.851 2 

.8434 

. 3356 

.3279 

.8203 

.8127 

-21.0 

.3 15? 

. 7973 

.7904 

.7831 

.7759 

. 7687 

.7615 

• 7545 

.7475 

. 74 J5 

-22.0 

.7136 

.7263 

.7200 

.7133 

« 7 C66 

. 7000 

.6935 

.6370 

.6806 

.6742 

-23,0 

.6679 

.6616 

.6554 

.6492 

.6431 

.6370 

.6310 

.6251 

.6192 

.6133 

-24.0 

.6075 

.6013 

.5961 

.5904 

.5848 

.5793 

. 5738 

. 5683 

,5629 

.5575 

-25.0 

.5522 

.5469 

.5417 

.5365 

.5314 

.5263 

.5213 

.5163 

,5113 

. 5 Li 6 4 

-26.0 

.5015 

.49 67 

.4919 

.4872 

.4825 

.4773 

.4732 

• 46 86 

,4641 

,4 59b 

-27.0 

.455 1 

.45 07 

.4463 

.4420 

.4377 

.4334 

.4292 

.4250 

.4209 

.4168 

-28.0 

.4127 

.4087 

.4047 

. 4007 

.3968 

.3929 

• 3 390 

.3352 

.3814 

. 3776 

-29.0 

. 3 73 9 

.370? 

.3666 

.3629 

,3594 

.3553 

.3523 

. 3483 

.3453 

. 3419 

-30.0 

.3335 

.33 51 

.3318 

.3285 

V325? 

.3219 

.3187 

.3155 

.3124 

. 3 093 

-31.0 

.3062 

. 3031 

.3000 

.2970 

.2940 

. 2911 

.2881 

• 28 52 

.2824 

.2715 

-32.0 

.2767 

. 2739 

.2711 

.2634 

.2656 

.2630 

.2603 

• 25 76 

.2553 

• 2524 

-33.0 

.2493 

.2473 

• 2448 

.2423 

.2398 

.2373 

.2349 

.2325 

.2301 

.2277 

-34.0 

.2254 

• 2231 

• 2203 

.2185 

.2163 

.2140 

.2118 

.2996 

.2 075 

• 2 0 5 3 

-35.0 

.2 332 

.2011 

.1990 

.1969 

.1948 

.1928 

.1908 

.1888 

, 1 868 

.1 349 

-36.0 

.1830 

• 18 10 

.1791 

.1773 

.1754 

.1736 

. 1717 

,1699 

.1681 

. 1 6t»4 

-37.0 

.1646 

.1629 

.1612 

.1595 

.1578 

.1561 

.1544 

. 1528 

.1512 

.1 496 

-38.0 

.1430 

• 14 64 

.1448 

.1433 

.1418 

• 1403 

.1383 

.1373 

.1358 

.1343 

- 39.0 

.1329 

.1315 

.1301 

.1237 

.1273 

.1259 

.1246 

.123? 

,1219 

. 12t»b 

-40.0 

.1193 

.1130 

.1167 

.1154 

.1142 

. 1129 

.1117 

.1105 

. 1 093 

.10 81 

-41.0 

. 1 0690 

.10573 

.10458 

.10343 

.11) ?30 

.10118 

.10007 

. U 9397 

.09788 

• osr 3u 

-42.0 

.09574 

• 09 46 8 

.09364 

.09260 

.09158 

.09057 

,08956 

.09857 

.0 8759 

,1,8462 

-43.0 

.0 35 65 

. 03470 

.03376 

.08283 

.08190 

.06099 

.08008 

,07919 

.07830 

.'7743 

-44.0 

.07656 

. 07 570 

.07485 

.07401 

.07318 

.07235 

.07154 

• u7 0 7 3 

• 0 6993 

• U b 4 1 4 

-45.0 

.06336 

.06759 

.06682 

• 0 6606 

.06531 

• 0645 7 

.06384 

,06 311 

,06239 

,06168 

-46.0 

.06098 

.0602? 

.05959 

.05391 

.05824 

.05757 

.05691 

.05625 

.05561 

.05497 

-47.0 

.05434 

. 05 3 71 

.05309 

.05248 

.05187 

.05127 

.05068 

• C 5009 

,04951 

.04894 

-48.0 

.3 43 37 

.04781 

.04725 

•04670 

.04615 

.04562 

.04508 

.04456 

.04434 

• u 4 35? 

-49.0 

.04301 

. 04251 

.04201 

.04151 

.04102 

.04054 

.04006 

.03959 

.03912 

.0 ?866 

-50.0 

. 033 20 

.03 775 

.03731 

.03636 

.03643 

.03599 

.03557 

. 03514 

,03472 

.03431 

-51.0 

. 0 3390 

. 03 350 

.03309 

.03270 

.03231 

.03192 

.03154 

• 03 11 b 

• u 3 vi 7 B 

• U 3d 41 

-52.0 

.03005 

• C2 969 

.02933 

.02397 

.02862 

.02823 

.02793 

• G2 760 

,02726 

• U 2693 

- 53*0 

.02660 

.02628 

.02596 

.02564 

.02533 

. 02502 

.02472 

.02441 

,02411 

.02332 

-54.0 

.02353 

. 02 324 

.02295 

.02267 

.02239 

.02212 

•Q21R4 

.02157 

.02131 

.02104 

-55.0 

.02073 

. 02053 

.02027 

.02002 

.01977 

.01952 

.01923 

. J1904 

•01880 

• U 1 8 57 

-56.0 

.01834 

.01311 

.01788 

.01766 

.01744 

•01722 

.01703 

.01679 

,01 658 

.01637 

-5 7.0 

.01616 

• 01 596 

.31576 

.01556 

.01536 

.01516 

.01497 

*01478 

.0 1459 

.01441 

-58.0 

.01423 

.01404 

.01387 

.01369 

.01351 

.01334 

.01317 

.01300 

.01283 

.01267 

-59.0 

.01251 

.01235 

.01219 

•01203 

.01188 

.01172 

.01157 

.01142 

.01127 

.11113 



Table B-IV 


DENSITY OF PURE WATER VAPOR AT SATURATION OVER WATER 


TEMPI Cl 

.0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

UNIT 

6/M 5 

G/H 3 

G/M 3 

G/H 3 

G/M 3 

G/M 3 

G/M 3 

G/M 3 

G/M 3 

o/M 3 

-0.0 

4.847 

4.814 

4.781 

4.748 

4.715 

4.683 

4.650 

4,618 

4.586 

4.554 

-1.0 

4.523 

4.491 

4.460 

4.429 

4.398 

4.368 

4.337 

4. 3j 7 

4.277 

4.24 7 

-?.o 

4.218 

4.188 

4.159 

4.133 

4.101 

4.072 

4. 0 43 

4.015 

3.986 

3.956 

— 3.0 

3.930 

3.933 

3.875 

3.848 

3.521 

3.793 

3.767 

3.740 

3.713 

3.607 

-4.0 

3.661 

3.63 5 

3.609 

3.583 

3.557 

3.532 

3.507 

3.482 

3.457 

3.4 32 

-5.0 

3.407 

3.38 3 

3.359 

3.334 

3,310 

3.207 

3.263 

3.239 

3.216 

3.193 

— 6.0 

3*170 

3.14 7 

3.124 

3.131 

3.079 

3.056 

3.034 

3.012 

2.990 

2.96 8 

-7.0 

2.947 

2.92 5 

2.904 

2.882 

2.861 

2.840 

2.020 

2.799 

2.778 

2.750 

-8.0 

2.738 

2.717 

2.697 

2.677 

2.658 

2.638 

2.618 

2.399 

2.580 

2.561 

-q.o 

2.542 

2.523 

2.504 

2.485 

2.467 

2.448 

2.430 

2.412 

2.394 

2.376 

-10.0 

2.358 

2.341 

2.323 

2.305 

2.288 

2. 271 

2.253 

2.236 

2.22Q 

2.2G3 

-11.0 

2.186 

2.17 0 

2.15? 

2.137 

2.121 

2.105 

2.089 

2.073 

2.057 

2.041 

-12.0 

2.026 

2. 010 

1.995 

1.930 

1 .965 

1.950 

1.935 

1.920 

1.905 

1.890 

— 13.0 

1.876 

1.861 

1.847 

1.833 

1.819 

1.805 

1.791 

1.777 

1.763 

1.7 49 

-14.0 

1.736 

1.722 

1.709 

1.696 

1.682 

1.669 

1.656 

1.643 

1.630 

J , 6 1 0 

-15.0 

1.605 

1.593 

1.580 

1.565 

1.555 

1.543 

1.531 

1-519 

1.507 

1.4 95 

-16.0 

1.483 

1.471 

1.460 

1.448 

1 • 437 

1.425 

1.414 

1.-4Q3 

1.392 

1.381 

-17.0 

1.3 70 

1.359 

1.345 

1.337 

1.326 

1.316 

1.305 

1.295 

1.2 84 

1.274 

-18.0 

1.264 

1.25 4 

1.243 

1.233 

1.223 

1.214 

1.204 

1.194 

1.184 

1.175 

-19.0 

1 . 1 65 

1.156 

1.146 

1.137 

1,128 

1.119 

1.109 

1.10 0 

1.091 

1.082 

-20.0 

1.0736 

1.0648 

1.3559 

1.0473 

1.0387 

1.0301 

1.0216 

1.0132 

1.0049 

.9966 

-21.0 

.9883 

.9802 

.9720 

.9640 

.9560 

.9401 

.9402 

,9324 

.9246 

.9169 

-?2.0 

• 9 092 

.9017 

.8941 

.8866 

.8792 

.0719 

.8646 

. 8573 

• 8501 

. 0429 

-23.0 

. 8359 

.’8288 

.8218 

.8149 

• 8080 

.8012 

.7944 

.7877 

.7810 

. 7744 

-24.0 

. 767 8 

.7613 

.7543 

.7484 

.7420 

.7356 

.7294 

.7231 

.7169 

.7 lu 1 

-25.0 

.7047 

.6987 

.6927 

.6067 

.6800 

.6749 

.6691 

.6633 

.6576 

.6519 

-26.0 

.6463 

.6407 

.6351 

.6296 

• 6242 

.6187 

.6134 

• 60 00 

.6027 

• 5 975 

-27.0 

.5922 

.5871 

.5519 

.5768 

.5718 

.5668 

.5618 

.5568 

.5519 

.5471 

-28.0 

.5423 

.5375 

.5327 

.5280 

.5233 

.5187 

.5141 

.5096 

• 5 G50 

• 5 L05 

-29.0 

.4961 

.4917 

.4573 

.4829 

.4786 

.4743 

.4701 

.4659 

.4617 

.4576 

— 30.0 

.4535 

.44 94 

.4453 

.4413 

.4373 

.4334 

* 4295 

.4256 

.4217 

.4179 

-31.0 

.4141 

• 4104 

.4066 

• 4029 

.3993 

.3956 

.3920 

.3884 

.3849 

.3014 

-32.0 

.3779 

.37 44 

.3710 

• 3676 

.3642 

. 3609 

.3575 

. 3542 

.3510 

.3477 

-33.0 

.3445 

• 3413 

.3382 

.3350 

.3319 

.3288 

. 3258 

. 3221 

• 3 1 98 

.3160 

-34.0 

.3138 

.'3109 

.3080 

.3051 

.3022 

. 2994 

.2966 

• 2938 

.2910 

. 280 1 

-35.0 

.2856 

.28 29 

.2802 

.2776 

.2749 

.2723 

.2697 

.26 72 

.2646 

,2621 

-36.0 

.2596 

.2572 

.2547 

.2523 

.2499 

. 2475 

.2451 

.2428 

.2404 

.2301 

-37.0 

.2359 

. 23 36 

.2313 

.2291 

,2269 

.2247 

.2225 

.2204 

.2183 

.2161 

— 38.0 

.2140 

.2120 

.2099 

.2079 

.2058 

.20 38 

.2018 

.1999 

.1979 

• 1 9 4 u 

-39.0 

.1941 

.1922 

.1903 

.1884 

.1865 

.1847 

.1829 

.1811 

.1793 

.1775 

-40.0 

.1750 

.1740 

.1723 

.1736 

.1689 

.1672 

.1656 

.1639 

.1623 

.1607 

-41.0 

.1591 

.15 75 

.1559 

.1543 

.1528 

.1512 

.1497 

.1482- 

.1467 

.1452 

-42.0 

.1438 

.14 23 

.1409 

.1395 

.1380 

.1366 

.1353 

.1339 

.1325 

.1312 

-4 3.0 

.1298 

.12 85 

.1272 

.1259 

.1246 

.1233 

.1221 

.1238 

.1196 

.1103 

-44.0 

.1171 

.1159 

.1147 

.1135 

.1124 

.1112 

.1101 

.1089 

.1078 

.1067 

-45.0 

.10555 

. 10 445 

.10 336 

.10228 

.10122 

•10016 

.09911 

.09007 

.09704 

,19603 

-46.0 

.09502 

. 09402 

.09303 

.09205 

.09108 

.09011 

.00916 

.08822 

,00728 

.01636 

-47.0 

.0 1544 

.08453 

.08363 

.08274 

•08186 

.08099 

,00012 

.07927 

.07842 

• U 775 8 

—4 8.0 

.0 76 75 

.07592 

.07511 

.07430 

.07350 

.07271 

.07192 

.07115 

.07038 

.06961 


.06886 

• 06811 

.06737 

• 06664 

.06592 

.06520 

,06449 

.06371 

« 063C9 

.06240 

-50.0 

.361 71 

• 06104 

.06037 

.05973 

.05905 

• 05840 

.06775 

. 05712 

•05649 

.05536 

—5 1.0 

.05524 

.05463 

.35403 

.05343 

.05283 

.05224 

.05166 

.05109 

.05052 

.04995 

-52.0 

.04939 

.04884 

.04829 

.04775 

.04722 

.04669 

.04616 

,04564 

.04513 

. 044 62 

-53.0 

•04411 

•04 361 

.34312 

.04263 

.04215 

.04167 

.04119 

.04072 

.04026 

.0 398 0 

-54.0 

.03935 

.03890 

.03845 

•03801 

.03757 

.03714 

.03672 

. 1 3629 

• G 3588 

.03546 

-55.0 

.03505 

.03465 

.03425 

» C 3385 

.03346 

• 03307 

.03268 

.03230 

.03193 

.03166 

—56.0 

.03119 

.03082 

.03046 

.03011 

.02975 

.02940 

.02906 

.02872 

.02038 

.12804 

-5 7.0 

.02771 

• 32 7 19 

.02706 

.02674 

.02642 

.02611 

.02583 

.02549 

.02519 

• 0 2409 

—5 8.0 

.02459 

.02430 

•02401 

.02372 

•02344 

.02316 

.02280 

.02263 

•02233 

• fj 22 G 6 

-59.0 

.02180 

.02153 

.02127 

•02101 

.02076 

•02051 

,32026 

•02001 

.01977 

.01953 


00 

00 


Table B-V 

PERCENT S . S . H . R . T . ICE AT WHICH WATER SATURATION OCCURS 


TEMPtO ) 

• 0 

.1 

.2 

.3 

.4 

.5 

• 6 

.7 

.0 

.9 

- 0.0 

.011 

.10 8 

.205 

• 303 

.400 

.497 

.595 

.692 

.7 90 

. 080 

* 1*2 

.99 

1.00 

1.10 

1.20 

1.30 

1.48 

1.57 

1.67 

1.77 

1.87 

- 2.0 

1.9 7 

2.07 

2.17 

2.27 

2.37 

2.47 

2.57 

2.67 

2.77 

2.87 

- 3.0 

2.96 

3.07 

3.17 

3.27 

3.37 

3.47 

3.57 

3.67 

3.7 7 

3.87 

- 4.0 

3.97 

4.07 

4.17 

4.27 

4.38 

4.48 

4.51 

4.68 

4.78 

4.88 

- 5.0 

4.99 

5.09 

5.19 

E . 29 

5. 40 

5.50 

5.60 

5.70 

5.81 

5.91 

- 6.0 

6.01 

6.12 

6 . 2 ? 

6.32 

6.43 

6.53 

6. 64 

6.74 

6 . H 4 

6.95 

- 7.0 

7.06 

7 . 16 

7.26 

7.37 

7.47 

7.57 

7.60 

7.78 

7.89 

7.99 

- 0.0 

8.10 

8.21 

0.31 

*•42 

1.52 

8.63 

8.73 

0.34 

0.95 

9.05 

- 9.0 

9.16 

9.27 

9.37 

9.48 

9.59 

9.69 

9.80 

9.91 

10.02 

10.12 

- 10.0 

10.2 

10.3 

10.4 

10.6 

10.7 

10.8 

10.9 

11.0 

11.1 

11.2 

- 11. 0 

11.3 

11.4 

11.5 

11.6 

11.7 

11.9 

12.0 

12.1 

12.2 

12.3 

- 12.0 

12.4 

12.5 

12.6 

12.7 

12.8 

13.0 

13.1 

13.2 

13.3 

13.4 

-i?.o 

13.5 

13.6 

13.7 

13.0 

14.0 

14.1 

14.2 

14.3 

14.4 

14.5 

- 14.0 

14.6 

14.7 

14.8 

15.0 

15.1 

15.2 

15.3 

IS. 4 

15.5 

15.6 

— 1 9 . Q 

15.7 

15.9 

16.0 

16.1 

16.2 

16.3 

16.4 

16.5 

16.7 

1 b . 8 

- 16.0 

16.9 

17 . 0 

17.1 

17.2 

17.3 

17.5 

17*6 

17.7 

17.0 

17.9 

- 17.0 

J 0.0 

18.1 

18.3 

18.4 

10.5 

18.6 

18.7 

10.0 

19.0 

19 . L 

- 18.0 

19.2 

19.3 

19.4 

19. 5 

19.7 

19.8 

19.9 

2 U . 0 

20.1 

20.2 

- 19.0 

20.4 

20.5 

20.6 

20.7 

20.0 

21.0 

21.1 

21.2 

21.3 

21.4 

- 20.0 


21.7 

21.8 

21.9 

22.0 

22.1 

22.3 

22.4 

22.5 

22.6 

- 21.0 

22.7 

22.9 

23.0 

23.1 

23.2 

23.3 

23.5 

23.6 

23.7 

23.8 

- 22.0 

23.9 

24.1 

24.2 

24.3 

24.4 

24.5 

24.7 

24.0 

24.9 

25.0 

- 23.0 

25.2 

25.3 

25.4 

25.5 

25.6 

25.8 

25.9 

2 6.0 

26.1 

26.3 

- 24.0 

26.4 

26.5 

26.6 

26.8 

25.9 

27.3 

27.1 

27.2 

27.4 

27.5 

- 25.0 

27.6 

27.7 

27.9 

28.0 

28.1 

28.2 

28.4 

28.5 

28.6 

28.7 

— 26.0 

28.9 

29.0 

29.1 

29.2 

29.4 

29.5 

29.6 

29.7 

29.9 

30.0 

- 27.0 

30.1 

30.2 

30.4 

30.5 

30.6 

30.8 

30.9 

31.0 

31.1 

31.3 

— 28.0 

31.4 

31.5 

31.6 

31.8 

31.9 

32.0 

32.2 

32.3 

32.4 

32.5 

- 29.0 

32.7 

32.8 

32.9 

33.1 

33.2 

33.3 

33.4 

33.6 

33.7 

33.8 

- 30.0 

34.0 

34.1 

34.2 

34.4 

34.5 

34.6 

34.7 

34.9 

35 .0 

35.1 

— 31.0 

35.3 

35.4 

35.5 

35.7 

35.8 

35.9 

36.0 

36.2 

36.3 

36.4 

— 32.0 

36.6 

36.7 

36.0 

37.0 

37.1 

37.2 

37.4 

37.5 

37.6 

37.8 

— 33.0 

37.9 

38.0 

38.2 

38.3 

31.4 

31.6 

38.7 

38.8 

39.0 

39.1 

- 34.0 

39.2 

39.4 

39.5 

39.6 

39.0 

39.9 

40.0 

40.2 

40 .3 

40.4 

- 35.0 

40.6 

40.7 

40.8 

41.0 

41.1 

41.2 

41.4 

41.5 

41.6 

41.8 

— 36.0 

41.9 

42.0 

42.2 

42.3 

42.5 

42.6 

42.7 

42.9 

43.0 

43.1 

- 37.0 

43.3 

43.4 

43.5 

43.7 

43.6 

44.0 

44.1 

44.2 

44.4 

44.5 

- 38.0 

44.6 

44.8 

44.9 

45.0 

45.2 

45.3 

45.5 

45.6 

45.7 

45.9 

— 39.0 

46 • 0 

46 . 1 

46.3 

46.4 

46.6 

46 . 7 

46.8 

47.0 

47.1 

47.3 

- 40.0 

47.4 

47.5 

47.7 

47.8 

47.9 

48.1 

48.2 

48.4 

48.5 

48.6 

- 41.0 

40.8 

48.9 

49.1 

49.2 

49.3 

49.5 

49.6 

49 . d 

49.9 

50.0 

- 42.0 

50.2 

50.3 

50.5 

50.6 

50.7 

50.9 

51.0 

51.2 

51.3 

51.4 

— 43.0 

51.6 

51.7 

51.9 

52.0 

52.1 

52.3 

52.4 

52.6 

52.7 

52.9 

- 44.0 

63.0 

53.1 

53.3 

53.4 

53.6 

53.7 

53.8 

54.0 

54.1 

54.3 

- 45.0 

54.4 

54.5 

54.7 

54.0 

55.0 

55.1 

55.3 

55.4 

55.5 

55.7 

- 46.0 

55.8 

5 6 .’ 0 

56.1 

56.3 

56.4 

56.5 

56.7 

56.1 

57.0 

57.1 

- 47.0 

57.2 

57.4 

57.5 

57.7 

57.6 

50.0 

58.1 

58.2 

58.4 

58.5 

- 48.0 

58.7 

56.1 

59.0 

59.1 

59.2 

59.4 

59.5 

59.7 

59.8 

bu . 0 

- 49.0 

60.1 

6 0.2 

60.4 

60.5 

60.7 

60.8 

61.0 

61.1 

61.2 

61.4 

-5 0.0 

61.5 

61.7 

61.8 

62.0 

62.1 

62.2 

62.4 

62.5 

62.7 

62.8 

- 51.0 

63.0 

6 3.1 

63.2 

63.4 

63.5 

63.7 

63.0 

64.0 

64.1 

64.2 

- 52.0 

64.4 

64.5 

64.7 

64.0 

65.0 

65.1 

65.2 

65.4 

65.5 

65.7 

- 53.0 

65.0 

66.0 

66*1 

66.2 

66.4 

66.5 

66.7 

66.8 

67.0 

67. 1 

- 54.0 

67.2 

67.4 

67.5 

67.7 

67.0 

60.0 

66.1 

68.2 

68.4 

6 R .5 

- 55.0 

68.7 

6 0.0 

68.9 

69.1 

69.2 

69.4 

69.5 

6 9.7 

69.8 

69. 9 

- 56.0 

70.1 

70.2 

70.4 

70.5 

70.6 

70.0 

70.9 

71.1 

71.2 

71.3 

- 57.0 

71.5 

71.6 

71.1 

71.9 

72.0 

72.2 

72.3 

72.5 

72.6 

72.7 

— 50.0 

72.9 

73.0 

73.2 

73.3 

73.4 

73.6 

73.7 

73.8 

74.0 

74.1 

- 59.0 

74.3 

74.4 

74.5 

74.7 

74.0 

75.0 

75.1 

75.2 

75.4 

75.5 
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APPENDIX C 

CALCULATION OF AMBIENT CONDITIONS 

Table C-I gives specimen calculations of the specifications for 
the crystal growth environment for a chamber top plate temperature of 
-8.9°C and a chamber bottom plate temperature of -22.1°C using a linear 
vapor density profile between the plates. The column headings from left 
to right denote the following: TEMP (C) is the temperature in °C, 

X (CM) is the height above the bottom plate in centimeters, TC (C) is 
the surface temperature in °C of a growing ice crystal calculated from 
the conduction-diffusion equilibrium relationship for a static environ- 
ment, DLT (C) is the difference in temperature in °C between the ambient 
air and the calculated crystal surface temperature, TC (C), at a given 
height in the chamber, MVDX is the maximum vapor density excess in g 

O 

m , TVDXS is the 'true* vapor density excess over a growing crystal 

with a surface temperature, TC (C), VDICE and VDWAT are, respectively, 

-3 

the saturated vapor densities in g m over a plane ice or water sur- 
face at the ambient temperature, SSICE and SSWAT are, respectively, the 
supersaturations in percent of the ambient environment with respect to 
a plane ice or water surface, WATSAT is the supersaturation with 
respect to ice at which water saturation occurs at that ambient tempera- 
ture, and finally, XVDW is the vapor density difference in g m be- 
tween VDWAT and VDICE. 

Figure C-l is a diagram showing the calculated maximum vapor den- 
sity excess, Ap max , in g nf^ which occurs in the vertical middle of the 
chamber, i.e., at x = 1.25 cm above the bottom plate, as a function of 
the top plate and midpoint of the chamber temperatures. The chamber 



186 


bottom plate temperature, Tj J0t , can be calculated from the equation 
T bot = 2 T mid “ T top 

where is the midpoint of the chamber temperature and T top is the 
chamber top plate temperature. Thus, from Fig. C-l the chamber top 
and bottom plate temperatures needed to produce any desired ambient 
conditions of temperature and vapor density excess at the vertical 
center of the chamber can be determined. 


Table c-I Specimen Calculations of the Ambient Chamber Conditions. 
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Fig. C-l Calculated ambient crystal growth conditions at 
the chamber center can be determined for any 
pair of top and bottom plate temperatures. 


Top Plate Temperature ( C) 
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APPENDIX D 

RELATIVE MAGNITUDES OF HEAT TRANSFER 

Heat transfer from a growing ice crystal surface takes place by 
conduction to the substrate, conduction to the air, convection to the 
air and radiation to the surroundings. The relative magnitudes of 
these terms were computed for this experiment. 

I. Heat Transfer by Conduction 

t 

Heat transfer by conduction is calculated from the equation 
dQ/dt = -KA dT/dx where 

dQ/dt = Heat transferred per unit time 

K = Thermal conductivity of the heat transferring media 

A = Cross-sectional area through which heat transfer occurs 

dT/dx = Temperature gradient in the heat transferring media. 

A. Heat Transfer by Conduction to the Substrate 

The substrate in this case consisted of the vertical solid 
glass fiber from which the ice crystals grew. Thus, conduction to the 
substrate means conduction along the length of the ice crystal to the 
glass fiber. Since in this calculation the maximum heat transfer to 
the substrate is desired, the assumption is made that heat transferred 
to the fiber at the position of the base of the crystal is efficiently 
dissipated away. In actuality, the solid glass fiber was a poor heat 
conductor and was coated along its length with growing ice crystals, so 
the actual heat transfer by conduction to the substrate was probably 
somewhat less than that calculated. 

The thermal conductivity of ice, K- ce , is given approximately 


by 
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K tce “ 5 x 10 ~ 3 cal cm_1 S_1 ° C_1 * 

For an ice crystal growing at an ambient supersaturation well above water 
saturation, the crystal tip may be a maximum of approximately 1.0°C 
warmer than either the ambient air or the glass fiber at that height in 
the chamber. Hence for an ice crystal 0.5 cm in length, dT/dx = 2°C 
cm”^. If the contact area of the crystal base is assumed to be as large 
as 10“ 4 cm 2 then 


dQ/dt = -1.0 x 10' 6 cal. s" 1 . 


B. Heat Transfer by Conduction to the Air 

The thermal conductivity of air, K air , is given approximately 
by 


K • = 4.0 x 10 -5 cal cm - ^ s“^ °C - ^ . 

a 1 1 


If one assumes a temperature gradient near the crystal of dT/dx = 

I Op Q I 

qj — = 10°C cm" and assumes conduction takes place from both sides 

2 2 2 

of a flat crystal of total surface area A = 2 (1mm) = 2 x 10 cm then 

dQ/dt = -8.0 x 10" 6 cal s" 1 . 


II . Heat Transfer by Convection 

The equation employed in this study to calculate heat transfer by 
convection was the one used by Hall and Pruppacher (1976). That is. 


dQ/dt = 4„CKF a (VT s )f Q , 


where 


J 


111 I I I II II III I III I 


in nr 


t — i H r 



191 


!*, Kf Q 

Q ^ a V aC p 
r* = S/(4 ttC) 



1.00 + 0.14( Np y 3 N R V L 2 * ) 2 , N p y 3 N R ’{ L 2 . <10 
0.86+0.28 N p y 3 N R y 2 *. .0 


K = K d n-(1.17-1.02 KyK/bpyPa" 1 ] 

K d = (5.69+0.168T) x 10" 5 cal cm" 1 s' 1 °C" 1 
K v = (3.73+0.020T) x 10" 5 cal cm" 1 s" 1 °C" 1 . 


The individual parameters are defined in the following way: 

C = capacitance = r for a spherical ice particle 
K = thermal conductivity of moist air 
= environmental temperature 
T s = temperature at ice particles' surface 
fq = heat ventilation factor 
S □ total surface area of ice particle 
P a = density of air 

v r = average velocity of the air molecules striking the ice 
particles' surface = (SR^/tt) 1 ^ 

R d □ specific gas constant for dry air 
a = thermal accommodation coefficient = 1.0 
Cp = specific heat of moist air at constant pressure 
N pp E Prandtl number [=nC p /K] = 0.72 

n 2 dynamic viscosity of air [= 1.8325 x 1 0“ 4 (T/296 . 1 6 ) 1 * 5 


(T+120)/(T+296.16) , with T in °K] 
N Re L * = Reynolds number based on L 
v = kinetic viscosity of air = .2 cm 2 



s 


-1 
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★ 

L = characteristic length defined as the total surface area 
divided by the perimeter normal to the flow for the ice 
particle 

P v = ambient water vapor density 

For chamber plate temperatures of T t0 p = -1.8°C and T bot = -19.9°C at 
an ambient pressure of p = 850 mb, i.e., the average pressure at Reno, 
Nevada, the ambient temperature, T^, at a height of x = 1.36 cm above 
the bottom plate is = -10.2°C. Other values there are p y = 2.664 g 
m~ 3 , T $ = -9.1°C, Ap max = 0.563 g m" 3 , Ap statl - C = 0.359 g m“ 3 , 

P a = 1 .1 kg m" 3 , 

R d = m = 2,871 x 1C)6 erg ° K_1 » 


Cp = c p Cdry air) + wCp (water vapor) 


w = mixing ratio = — = 10' 2 to 10' 3 so C = 0.24 cal g -1 °C -1 , at 


T = -10°C 


v r " 


7T 


1/2 


K d = 4.01 x IQ" 5 cal cm' 1 s' 1 °C" 1 
K y = 3.53 x 10' 5 cal cm' 1 s" 1 °C' 1 
K = 4.01 x 10" 5 cal cm' 1 s' 1 °C' 1 


= 4.39 x 10 4 cm s' 1 . 


For a sphere of radius 1 mm 




* 


r = 0.1 cm and L 


2 

_ 4irr 
2irr 


2r = 0.2 cm. 


N Re,L* ■ = V for v = 5.0 cm s' 1 , N Re _ L * = 5.0 

and N p J /3 N R 172 * = 2.004 so f Q = 1.421, 1* = 1.967 x 10' 5 cm 

so F„ = 0.9998 and dQ/dt = -7.87 x 10' 5 cal s” 1 for v = 5.0 cm s' 1 . 

(X 
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III. Hea t Transfer by Radiation 

The net radiative heat flux $ R is equal to the difference between 
the amount of impinging radiation upon the ice particles' surface 
which is absorbed by the ice particle, and the radiation that is emitted 
by the particles' surface, i.e., 

* R = (1-TOD TO7-effT s 4 , 


where R(X) is the reflectivity, A(xJ the absorptivity and e the emis- 
sivity of ice. 


TOT = 



B(T,X)R(X)dX 


0.034 


TOT = 


5F 


B(T,X)A(X)dX = 1.0 


where B(T,X) is the Planck spectrum function. 


e = 0.95 

a = Stefan constant = 1.355 x 10”^ 3 cal cm" 3 °K” 4 s"^ 


T 4 + T 4 

bot + top 


$ R = 0.966 - 0.95 a T $ 4 but = 1/2 ct 

For T top = “ 1 * 8 ° c = 271.36°K and T bot = -19.9°C = 253.26°K 

= 6.4609 x 10” 3 cal cm" 3 s"^ . At a height, x, in the chamber of 
x = 1.36 cm the ambient temperature was T ro = -10.2°C = 262.96°K and 
the ice crystal surface temperature was computed from conduction- 
diffusion equilibrium to be T = -9.1°C = 264.06°K. Thus, 


$ R = [6.2412 x 10” 3 - 6.2585 x 10" 3 ] cal cm” 3 s”^ 
or 

$ R = -1.7 x 10" 3 cal cm" 3 s"^ . 
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The term, $ R , is very sensitive to the ice crystal surface temperature, 
T $ , being much larger or even positive for a surface temperature only 
0.2°C colder. 

The resultant heat transfer by radiation is given by dQ/dt = S $ R 

O 

where S is the radiating surface area of the crystal. For S = 2 x 10 

p 

cm , i.e., the same as assumed for conduction to the air, 
dQ/dt = -3.4 x 10“ ^ cal s“^ . 
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